WO2021193296A1 - 日射推定システム、空調制御システム、空調装置、車両、建造物、日射推定方法、空調制御方法、及び日射推定プログラム - Google Patents

日射推定システム、空調制御システム、空調装置、車両、建造物、日射推定方法、空調制御方法、及び日射推定プログラム Download PDF

Info

Publication number
WO2021193296A1
WO2021193296A1 PCT/JP2021/010883 JP2021010883W WO2021193296A1 WO 2021193296 A1 WO2021193296 A1 WO 2021193296A1 JP 2021010883 W JP2021010883 W JP 2021010883W WO 2021193296 A1 WO2021193296 A1 WO 2021193296A1
Authority
WO
WIPO (PCT)
Prior art keywords
solar radiation
amount
total
solar
sky
Prior art date
Application number
PCT/JP2021/010883
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 JP2022510009A priority Critical patent/JPWO2021193296A1/ja
Publication of WO2021193296A1 publication Critical patent/WO2021193296A1/ja

Links

Images

Classifications

    • B—PERFORMING OPERATIONS; TRANSPORTING
    • B60—VEHICLES IN GENERAL
    • B60H—ARRANGEMENTS OF HEATING, COOLING, VENTILATING OR OTHER AIR-TREATING DEVICES SPECIALLY ADAPTED FOR PASSENGER OR GOODS SPACES OF VEHICLES
    • B60H1/00—Heating, cooling or ventilating [HVAC] devices
    • 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/30—Control or safety arrangements for purposes related to the operation of the system, e.g. for safety or monitoring
    • 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/70—Control systems characterised by their outputs; Constructional details thereof
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
    • G01J1/00—Photometry, e.g. photographic exposure meter
    • G01J1/42—Photometry, e.g. photographic exposure meter using electric radiation detectors
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01W—METEOROLOGY
    • G01W1/00—Meteorology
    • G01W1/12—Sunshine duration recorders

Definitions

  • an air conditioner control device that is installed in a vehicle equipped with a navigation device and controls a vehicle air conditioner.
  • the air conditioner control device determines the presence or absence of solar radiation on the vehicle and the direction of intrusion based on the solar radiation sensor arranged in front of the vehicle and the road information including the solar position information obtained from the sun position search map and the information on the buildings around the road. It is provided with a solar radiation direction calculation means for calculating the amount of solar radiation, and a solar radiation amount calculation means for obtaining the amount of solar radiation based on the calculation result of the solar radiation direction calculation means and the detection value of the solar radiation sensor.
  • the solar radiation amount calculation means corrects the solar radiation amount based on the information unique to the vehicle (see, for example, Patent Document 1).
  • the XYZ coordinate system is defined as the coordinate system of the vehicle 10.
  • the X direction is the front-rear direction (length direction) of the vehicle 10, and the rearward direction is the + X direction.
  • the Y direction is the lateral direction (width direction) of the vehicle 10, and the right direction is the + Y direction.
  • the Z direction is the vertical direction of the vehicle 10, and the upward direction is the + Z direction.
  • FIG. 4 is a diagram showing a solar azimuth angle A and a solar zenith angle h.
  • the solar azimuth angle A (degree) represents the azimuth angle of the sun S at a certain time (year, month, day, hour, minute, second) on a certain date, and is east of the south direction. Is represented by a positive value and the west side is represented by a negative value.
  • the solar zenith angle h (degree) is an altitude angle in the direction in which the sun S exists.
  • FIG. 5 is a diagram showing the incident angle ⁇ of direct solar radiation on the surface 20.
  • FIG. 5 shows the solar azimuth angle A and the solar zenith angle h as in FIG.
  • the inclination angle of the surface 20 is ⁇ (degree), and the azimuth angle is A'.
  • the incident angle ⁇ (degrees) of direct sunlight on the surface 20 is an angle formed by the direction of the normal vector n of the surface 20 and the direction specified by the solar azimuth angle A and the solar altitude angle h.
  • FIG. 7 is a diagram showing the solar radiation (total of direct solar radiation, sky solar radiation, and ground reflection solar radiation) incident on the window glass 30.
  • the solar radiation incident on the window glass 30 is divided into the amount transmitted through the window glass 30 (transmission amount), the amount reflected on the surface of the window glass 30 (reflection amount), and the amount absorbed by the window glass 30 (absorption amount). Be done.
  • the absorbed amount raises the temperature of the window glass 30 and re-radiates heat from the window glass 30 to the outside and the inside.
  • the roof 11 is used as a reference plane.
  • the vehicle 10 is an example of an object
  • the roof 11 is an example of a reference surface on the outer surface of the object.
  • each of the windowpanes 12WS to 12RG is an example of a predetermined outer surface of the outer surface of the object.
  • the total solar radiation amount on the roof 11 which is the reference plane is obtained before calculating the solar radiation heat acquisition amount in the windowpanes 12WS to 12RG.
  • the amount of solar radiation in the sky is taken into consideration in addition to the amount of direct solar radiation and the amount of ground reflection solar radiation.
  • the value obtained by multiplying the amount of sky solar radiation by the sky solar radiation transmittance SRT is added.
  • the convergence calculation is performed so that the total amount of solar radiation obtained for the roof 11 substantially matches the total amount of solar radiation obtained by the pyranometer.
  • the sky solar transmittance SRT and the atmospheric transmittance P are corrected to give an accurate total solar radiation amount that is substantially the same as the total solar radiation amount obtained by the pyranometer.
  • the corrected sky solar radiation transmittance CSRT and the corrected atmosphere. Obtain the transmittance CP.
  • the atmospheric transmittance P represents the ratio of the solar radiation before entering the atmosphere that passes through the atmosphere.
  • the corrected sky solar transmittance CSRT and the corrected atmospheric transmittance CP obtained for the roof 11 which is the reference plane are used for the calculation of the total solar radiation amount in each of the windowpanes 12WS to 12RG.
  • the roof 11 and the windowpanes 12WS to 12RG, which are the measurement target surfaces, are very close to each other, so that it is considered that a common corrected sky solar transmittance CSRT and a corrected atmospheric transmittance CP may be used.
  • I 0 is the solar constant
  • R is the radius of the earth
  • P is the atmospheric transmittance of direct solar radiation
  • h is the solar zenith angle.
  • ⁇ is the incident angle of direct solar radiation on the roof 11, and the surface 20 shown in FIG. 5 may be replaced with the roof 11.
  • the solar constant I 0 the value described in the science chronological table or the like may be used as an example. Since the value of the solar constant I 0 may be updated, the updated value may be used when the value is updated.
  • the calculation method described in Document (1) "Technical Explanation of Meteorological Data System Co., Ltd., 2016. Calculation of Sun Position, by Shinichi Matsumoto” is used. Just do it.
  • the atmospheric transmittance P may be obtained by the method described in Document (2) "Direct scattering separation and calculation of slope solar radiation, Meteorological Data System Technorogy Explanation Co., Ltd., 2017. by Hiroshi Akasaka" as an example.
  • an appropriate initial value is set.
  • An appropriate initial value of the atmospheric transmittance P is 0.7 as an example.
  • cos ⁇ may be obtained by the formula (5) described later.
  • the calculation unit 120 uses the value for the roof 11 at the time of acquisition of the total solar radiation amount by the pyranometer 101 for cos ⁇ . This is because cos ⁇ changes depending on the date and time.
  • the position data (including the time data) used when obtaining cos ⁇ the time data included in the position data acquired from the navigation ECU 300B or the time data used by the HVAC-ECU 200 may be used.
  • the main control unit 110 reads the total solar heat acquisition amount CQ TS from the memory 160 and transmits it to the control unit 210.
  • the control unit 210 stores the total solar heat acquisition amount CQ TS in the memory 220 and uses it for controlling the HVAC unit 470.
  • Window glass 12FL is the left front of the side glass, per solar radiation (Asahi) is in the morning at the time of the experiment, since the solar radiation will not hit almost become a shadow In the afternoon, in the morning and the calculated value J TS of the global solar radiation The measured value is high, gradually decreasing toward noon, and low in the afternoon. Since the window glass 12FL is not exposed to solar radiation in the afternoon, it is considered that the total amount of solar radiation obtained in the afternoon is mainly the amount of solar radiation from the sky. Thus, it can be seen that it is possible to obtain the global solar radiation J TS (calculated) for the window glass 12FL.
  • the direct solar radiation amount JD (calculated value) obtained for the window glass 12WS by the solar radiation estimation system 100 on October 31, 2019, and the solar radiation meter 101 are arranged on the window glass 12WS separately. It shows the measured value of the total amount of solar radiation obtained by the pyranometer. The measured value shown in FIG. 15 is the same as the measured value shown in FIG. In FIG. 15, the direct solar radiation amount JD (calculated value) and the measured value of the total solar radiation amount are compared.
  • the total solar radiation amount J TS (calculated value) and the measured value of the total solar radiation amount obtained by the pyranometer are almost the same, and the total solar radiation amount J TS (calculated value) can be obtained for the roof 11. You can see that.
  • Figure 19 shows the global solar radiation calculations (J TS) and the correction air permeability CP and correction sky solar transmittance CSRT when the measured values are acquired as shown in FIG. 18. From the start of measurement to, for example, around 12:27:06, the corrected atmospheric transmittance CP is zero, and the corrected sky solar transmittance CSRT changes between about 0.5 and about 1.0. There is. At this time, there was no direct solar radiation. In addition, after 12:27:06, the value of the corrected atmospheric solar transmittance CP may partially increase because there is a time zone when there is direct solar radiation, and the value of the corrected sky solar transmittance CSRT is It was higher than the value before 12:27:06. In FIG. 19, it can be seen that the corrected atmospheric transmittance CP and the corrected sky solar transmittance CSRT required by the solar radiation estimation system 100 are changing from moment to moment.
  • solar radiation estimation system 100 can determine the global solar radiation J TS entering the interior of the vehicle 10 with high accuracy.
  • Global solar radiation J TS is an estimated value obtained using the corrected atmospheric transmittance CP and correction sky solar radiation transmittance CSRT.
  • the HVAC unit 470 can be controlled very precisely.
  • the HVAC unit 470 When controlling the HVAC unit 470 using the outside air temperature and the room temperature of the vehicle 10, the HVAC unit 470 is operated after detecting an increase or decrease in the outside air temperature and the room temperature. At this time, due to the influence of the heat capacity of the window glass 12WS to 12RG and the air in the room, there is a delay between the time when the sunlight hits the room of the vehicle 10 and the time when the room temperature rises. Therefore, when the HVAC unit 470 is controlled by using the outside air temperature and the room temperature of the vehicle 10, there is a delay of about 5 minutes after the sunlight hits.
  • the solar radiation estimation system 100 the total amount of solar radiation heat gained in the room of the vehicle 10 due to solar radiation can be estimated instantly and accurately, so that the HVAC unit 470 can be operated at an early stage and the delay can be suppressed.
  • Room temperature can be controlled quickly and appropriately.
  • the HVAC unit 470 can be operated at an early stage and the room temperature can be controlled while the degree of increase in the room temperature is small, the operating amount of the HVAC unit 470 can be reduced, which can contribute to the improvement of the fuel consumption and the cruising range of the vehicle 10.
  • the vehicle 10 is an EV, the influence of the ground reflection solar radiation on the battery life can be grasped more accurately.
  • the direct solar radiation amount J D , the sky solar radiation amount J S, and the ground reflection solar radiation amount J G are added to obtain the total solar radiation amount J TS. Since the ratio of the ground reflection solar radiation amount J G is about several% to about 10%, the direct solar radiation amount J D and the sky solar radiation amount J S are added without obtaining the ground reflection solar radiation amount J G. it may be obtained global solar radiation J TS Te.
  • the reference plane is the roof 11, a solar radiation meter 101 arranged on the roof 11, the calculator 120 has been described embodiment to calculate the global solar radiation J ST in the roof 11.
  • the reference surface may be any surface other than the roof 11 as long as it is the outer surface of the vehicle 10.
  • Calculator 120 may calculate the global solar radiation J ST in the reference plane of the vehicle 10.
  • the conversion rate data may be any data representing the conversion rate given by the reciprocal of the transmittance of the window glass 12WS.
  • the transmittance of the window glass 12WS represents the ratio of the permeation amount to the amount of solar radiation on the outer surface of the window glass 12WS, and is the same as the ratio of the permeation amount to the amount of solar radiation on the outer surface of the window glass 30 shown in FIG.
  • a pyranometer 101 is arranged on the roof of the ship 10C shown in FIG. 25 (C).
  • the electronic control device having the main control unit 110, the calculation unit 120, the correction unit 130, the estimation unit 140, the calculation unit 150, and the memory 160 may be arranged inside the ship 10C as an example.
  • the total solar heat gain amount CQ representing the amount of heat entering the interior of a plurality of window glass 12C
  • the TS may be obtained and used for controlling the air conditioner of the ship 10C. Controlling the air conditioner based on the total amount of solar radiation gained CQ TS means controlling the air conditioner based on the total amount of solar radiation J TS for each of the plurality of windowpanes 12C.
  • the solar radiation estimation system, the air conditioning control system, the air conditioner, the vehicle, the building, the solar radiation estimation method, the air conditioning control method, and the solar radiation estimation program according to the exemplary embodiment of the present invention have been described above.
  • the present invention is not limited to the disclosed embodiments, and various modifications and changes can be made without departing from the scope of claims.
  • the acquisition unit that acquires the first total amount of solar radiation received on the reference plane on the outer surface of the object, Based on the solar constant, atmospheric transmittance, sky solar radiation transmittance, and the first relational value related to the orientation of the reference plane, the second total on the reference plane at the time of acquisition of the first total solar radiation amount.
  • a calculation unit that calculates the amount of solar radiation
  • a correction unit that corrects the atmospheric transmittance and the sky solar radiation transmittance, respectively, to obtain the corrected atmospheric transmittance and the corrected sky solar radiation transmittance so that the second total solar radiation amount becomes equal to the first total solar radiation amount.
  • the acquisition unit (A) A pyranometer that measures the first total amount of solar radiation on the reference plane. (B) A pyranometer for measuring the total amount of solar radiation transmitted through the transparent wall body of the object, and the above-mentioned first solar radiation amount measured by the solar radiation meter based on the optical characteristics of the transparent wall body.
  • Item 1 to 8 which is any one of (1) an acquisition unit having a conversion unit for converting to total solar radiation, and
  • C a reception unit for receiving data representing the first total solar radiation.
  • the solar radiation estimation system according to any one of items 1 to 9 and Air conditioner and An air conditioning control system including an air conditioning control unit that controls air conditioning of the air conditioning device based on the estimated total solar radiation amount estimated by the solar radiation estimation system.
  • the object is a vehicle
  • the solar radiation estimation system acquires data from any one or more of an electronic control device, a position data acquisition unit, a navigation control unit, and a communication unit mounted on the vehicle, and uses the acquired data to obtain the data.
  • the air conditioning control system according to item 11, which calculates an estimated total amount of solar radiation.
  • An air conditioner including the calculation unit, the correction unit, and a control device having the estimation unit of the solar radiation estimation system according to any one of items 1 to 9.
  • (Item 14) A vehicle including the solar radiation estimation system according to any one of items 1 to 9.
  • (Item 15) A ship comprising the solar radiation estimation system according to any one of items 1 to 9.
  • (Item 16) An aircraft comprising the solar radiation estimation system according to any one of items 1 to 9.
  • (Item 17) A building that includes the solar radiation estimation system according to any one of items 1 to 9.
  • (Item 18) Obtain the first total amount of solar radiation received on the reference plane on the outer surface of the object, Based on the solar constant, atmospheric transmittance, sky solar radiation transmittance, and the first relational value related to the orientation of the reference plane, the second total on the reference plane at the time of acquisition of the first total solar radiation amount.
  • the calculation is to calculate the solar heat gained amount for each of the plurality of predetermined outer surfaces, and calculate the total solar heat gained amount which is the sum of the plurality of solar heat gained amounts, according to item 24.
  • Insolation estimation method (Item 26) Obtaining the first total insolation amount is (A) To acquire the total amount of solar radiation measured by a pyranometer provided on the reference plane as the first total amount of solar radiation. (B) The total amount of solar radiation transmitted through the transmission wall of the object is measured, and the measured total amount of solar radiation is measured based on the optical characteristics of the transmission wall. Item 18 which is one of obtaining the first total insolation amount by converting to a quantity and (C) receiving data representing the first total insolation amount.
  • the first total solar radiation is based on the solar constant, the corrected atmospheric transmittance, the corrected sky solar radiation transmittance, and a second relational value related to the orientation of a predetermined outer surface of the outer surface of the object.
  • a solar radiation estimation program that causes a computer to perform processing including calculating an estimated total solar radiation amount, which is an estimated value of the total solar radiation amount on the predetermined outer surface at the time of acquisition of the amount.

Landscapes

  • Engineering & Computer Science (AREA)
  • Mechanical Engineering (AREA)
  • General Engineering & Computer Science (AREA)
  • Chemical & Material Sciences (AREA)
  • Environmental & Geological Engineering (AREA)
  • Physics & Mathematics (AREA)
  • Combustion & Propulsion (AREA)
  • Atmospheric Sciences (AREA)
  • Ecology (AREA)
  • Environmental Sciences (AREA)
  • Biodiversity & Conservation Biology (AREA)
  • Thermal Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Air-Conditioning For Vehicles (AREA)
  • Photometry And Measurement Of Optical Pulse Characteristics (AREA)

Abstract

車両等の内部に入る日射の全天日射量を高精度に推定できる日射推定システム等を提供する。　日射推定システムは、対象物の外面にある基準面で受ける日射の第１全天日射量を取得する取得部と、太陽定数と、大気透過率と、天空日射透過率と、前記基準面の向きに関係する第１関係値とに基づいて、前記第１全天日射量の取得時の前記基準面における第２全天日射量を算出する算出部と、前記第２全天日射量が前記第１全天日射量に等しくなるように、前記大気透過率及び前記天空日射透過率をそれぞれ補正して補正大気透過率及び補正天空日射透過率を求める補正部と、前記太陽定数と、前記補正大気透過率及び前記補正天空日射透過率と、前記対象物の外面のうちの所定外面の向きに関係する第２関係値とに基づいて、前記第１全天日射量の取得時の前記所定外面における全天日射量の推定値である推定全天日射量を算出する推定部とを含む。

Description

日射推定システム、空調制御システム、空調装置、車両、建造物、日射推定方法、空調制御方法、及び日射推定プログラム
　本発明は、日射推定システム、空調制御システム、空調装置、車両、建造物、日射推定方法、空調制御方法、及び日射推定プログラムに関する。
　従来より、ナビゲーション装置を備えた車両に設けられ車両用空調装置を制御するエアコン制御装置がある。エアコン制御装置は、車両前方に配設された日射センサと、太陽位置検索マップから求めた太陽位置情報と道路周辺の建造物の情報を含む道路情報とに基づいて車両に対する日射の有無および侵入方向を演算する日射方向演算手段と、日射方向演算手段の演算結果と日射センサの検出値とに基づいて日射量を求める日射量演算手段とを備える。日射量演算手段は、車両独自の情報に基づいて、日射量を補正する（例えば、特許文献１参照）。
特開２００２－３６２１２９号公報
　ところで、従来のエアコン制御装置は、天空日射を考慮していないため、最終的に求められる日射量は、実際の全天日射量との差が大きく、精度が低い。
　そこで、車両、船舶、航空機、又は建造物の内部に入る日射の全天日射量を高精度に推定できる、日射推定システム、空調制御システム、空調装置、車両、建造物、日射推定方法、空調制御方法、及び日射推定プログラムを提供することを目的とする。
　本発明の実施の形態の日射推定システムは、対象物の外面にある基準面で受ける日射の第１全天日射量を取得する取得部と、太陽定数と、大気透過率と、天空日射透過率と、前記基準面の向きに関係する第１関係値とに基づいて、前記第１全天日射量の取得時の前記基準面における第２全天日射量を算出する算出部と、前記第２全天日射量が前記第１全天日射量に等しくなるように、前記大気透過率及び前記天空日射透過率をそれぞれ補正して補正大気透過率及び補正天空日射透過率を求める補正部と、前記太陽定数と、前記補正大気透過率及び前記補正天空日射透過率と、前記対象物の外面のうちの所定外面の向きに関係する第２関係値とに基づいて、前記第１全天日射量の取得時の前記所定外面における全天日射量の推定値である推定全天日射量を算出する推定部とを含む。
　車両、船舶、航空機、又は建造物の内部に入る日射の全天日射量を高精度に推定できる、日射推定システム、空調制御システム、空調装置、車両、建造物、日射推定方法、空調制御方法、及び日射推定プログラムを提供できる。
実施の形態で用いる直達日射、天空日射、地面反射日射を説明する図である。 車両のルーフと窓ガラスの配置を示す図である。 車両のルーフと窓ガラスの配置を示す図である。 太陽方位角と太陽高度角を示す図である。 直達日射の入射角θを示す図である。 形態係数ＶＦＳを説明する図である。 窓ガラスに入射する日射を示す図である。 実施の形態の日射推定システムを車両に搭載した構成の一例を示す図である。 日射推定プログラムによって実現される処理を示すフローチャートである。 日射推定システムでの実験結果を示す図である。 日射推定システムでの実験結果を示す図である。 日射推定システムでの実験結果を示す図である。 日射推定システムでの実験結果を示す図である。 日射推定システムでの実験結果を示す図である。 日射推定システムでの実験結果を示す図である。 日射推定システムでの実験結果を示す図である。 日射推定システムでの実験結果を示す図である。 日射推定システムでの実験結果を示す図である。 日射推定システムでの実験結果を示す図である。 日射推定システムでの実験結果を示す図である。 日射推定システムでの実験結果を示す図である。 日射推定システムでの実験結果を示す図である。 日射推定システムでの実験結果を示す図である。 実施の形態の変形例の日射推定システムを車両に搭載した構成の一例を示す図である。 建造物、航空機、船舶を示す図である。
　以下、本発明の日射推定システム、空調制御システム、空調装置、車両、建造物、日射推定方法、空調制御方法、及び日射推定プログラムを適用した実施の形態について説明する。
　＜実施の形態＞
　図１は、実施の形態で用いる直達日射、天空日射、地面反射日射を説明する図である。図１では、地面１を平面、天空２を半球状の空間として考える。直達日射は、太陽Ｓから放射される太陽放射（日射）のうち、直接地上に到達する成分である。天空日射は、地球の大気圏外に到達した太陽放射のうち、大気中の分子や浮遊粒子によって散乱されて地上に到達する成分である。地面反射日射とは、直達日射及び天空日射が地面で反射する成分である。
　以下では、直達日射、天空日射、地面反射日射による日射量をそれぞれ直達日射量ＪＤ、天空日射量ＪＳ、地面反射日射量ＪＧとする。図１に示すように、一例として車両１０の窓を介して車内に入る日射量を計算する際には、直達日射量ＪＤ、天空日射量ＪＳ、地面反射日射量ＪＧを考慮する。なお、地上に到達するすべての日射量を合わせた全天日射量のうち、直達日射量ＪＤは約７割、天空日射量ＪＳは約３割、地面反射日射量ＪＧは数パーセント程度である。
　図２及び図３は、車両１０のルーフ１１、窓ガラス１２ＷＳ、１２ＦＲ、１２ＲＲ、１２ＦＬ、１２ＲＬ、１２ＲＧ（以下、窓ガラス１２ＷＳ～１２ＲＧと記す）の配置を示す図である。窓ガラス１２ＷＳ～１２ＲＧの各々は、日射を透過する透過壁体の一例である。図２には、車両１０の左斜め前方から、ルーフ１１、窓ガラス１２ＷＳ、１２ＦＬ、１２ＲＬ、１２ＲＧを立体的に示し、図３には、ルーフ１１及び窓ガラス１２ＷＳ～１２ＲＧを平面的に示す。
　ここで、車両１０の座標系としてＸＹＺ座標系を定義する。Ｘ方向は、車両１０の前後方向（長さ方向）であり、後ろ向きが＋Ｘ方向である。Ｙ方向は、車両１０の横方向（幅方向）であり、右方向が＋Ｙ方向である。Ｚ方向は、車両１０の上下方向であり、上方向が＋Ｚ方向である。
　ルーフ１１は、車両１０の屋根であり、窓ガラス１２ＷＳ(Windshield)は車両のフロントガラスであり、窓ガラス１２ＦＲは車両１０の右前のサイドガラス、窓ガラス１２ＲＲは車両１０の右後ろのサイドガラスである。窓ガラス１２ＦＬは車両１０の左前のサイドガラス、窓ガラス１２ＲＬは車両１０の左後ろのサイドガラス、窓ガラス１２ＲＧは車両１０のリアガラスである。図２では、窓ガラス１２ＦＲ、１２ＲＲがルーフ１１で見えないが、各矢印は法線ベクトルを示す。ルーフ１１及び窓ガラス１２ＷＳ～１２ＲＧの法線ベクトルは、一例として、それぞれ、ルーフ１１及び窓ガラス１２ＷＳ～１２ＲＧの外表面の中心点における法線ベクトルである。車両１０は、ルーフ１１及び窓ガラス１２ＷＳ～１２ＲＧの法線ベクトルの座標と向きを車両１０のＸＹＺ座標系で表すデータを有する。
　図４は、太陽方位角Ａと太陽高度角ｈを示す図である。図４に示すように東西南北を定義する場合に、太陽方位角Ａ（度）は、ある日付のある時刻（年月日時分秒）における太陽Ｓの方位角を表し、南方向に対して東側を正の値、西側を負の値で表す。太陽高度角ｈ（度）は、太陽Ｓが存在する方位における高度角である。
　図５は、面２０における直達日射の入射角θを示す図である。図５には、図４と同様に太陽方位角Ａと太陽高度角ｈを示す。面２０の傾斜角をβ（度）、方位角をＡ’とする。面２０における直達日射の入射角θ（度）は、面２０の法線ベクトルｎの方向と、太陽方位角Ａ及び太陽高度角ｈによって特定される方向とがなす角度である。
　図６は、形態係数ＶＦＳを説明する図である。形態係数ＶＦＳは、天空２から面２０に到達する日射の割合を表し、０から１．０の値を取る。図６（Ａ）に示すように、面２０が水平であり、傾斜角βが０度のときには、形態係数ＶＦＳは１．０である。図６（Ｂ）に示すように傾斜角βが９０度で面２０が地面１に対して垂直である場合には、形態係数ＶＦＳは０．５である。
　図７は、窓ガラス３０に入射する日射（直達日射、天空日射、地面反射日射の合計）を示す図である。窓ガラス３０に入射する日射は、窓ガラス３０を透過する分（透過量）、窓ガラス３０の表面で反射される分（反射量）、窓ガラス３０に吸収される分（吸収量）に分けられる。吸収量は、窓ガラス３０の温度を上昇させて、窓ガラス３０から外側及び内側に再放熱される。
　ここで、実施の形態の日射推定方法によるルーフ１１における全天日射量の算出と、窓ガラス１２ＷＳ～１２ＲＧにおける日射熱取得量の計算について説明する。ここでは、一例として、ルーフ１１を基準面とする。車両１０は、対象物の一例であり、ルーフ１１は、対象物の外面にある基準面の一例である。また、窓ガラス１２ＷＳ～１２ＲＧの各々は、対象物の外面のうちの所定外面の一例である。
　実施の形態の日射推定方法では、窓ガラス１２ＷＳ～１２ＲＧにおける日射熱取得量を計算する前に、基準面であるルーフ１１における全天日射量を求める。また、全天日射量を求める際に、直達日射量及び地面反射日射量に加えて、天空日射量を考慮する。その際に、直達日射量及び地面反射日射量に天空日射量をそのまま加算するのではなく、天空日射量に天空日射透過率ＳＲＴを乗じた値を加算する。
　天空日射透過率ＳＲＴとは、天空日射が大気圏を透過して地面１に届く割合（透過率）である。雲が無い青空のときと、雲がある曇り空のときとでは、天空日射が地面１に届く割合は異なる。天空日射透過率ＳＲＴは、このような天空日射の透過率を表し、０から１．０の値を取る。
　そして、ルーフ１１について求めた全天日射量が、日射計で得られる全天日射量と略一致するように収束計算を行う。収束計算では、天空日射透過率ＳＲＴと、大気透過率Ｐとを補正し、日射計で得られる全天日射量と略一致する正確な全天日射量を与える補正天空日射透過率ＣＳＲＴと補正大気透過率ＣＰとを求める。なお、大気透過率Ｐは、大気圏に入射する前の日射のうち、大気圏を通過する割合を表す。
　そして、基準面であるルーフ１１について求めた補正天空日射透過率ＣＳＲＴと補正大気透過率ＣＰを窓ガラス１２ＷＳ～１２ＲＧの各々における全天日射量の計算に用いる。ルーフ１１と、測定対象面である窓ガラス１２ＷＳ～１２ＲＧとは、位置が非常に近いため、共通の補正天空日射透過率ＣＳＲＴ及び補正大気透過率ＣＰを用いてよいと考えられるからである。補正天空日射透過率ＣＳＲＴ及び補正大気透過率ＣＰを用いて、ルーフ１１と窓ガラス１２ＷＳ～１２ＲＧの法線ベクトルの方向の違い等を考慮して、窓ガラス１２ＷＳ～１２ＲＧの各々における正確な全天日射量を求める。
　そして、窓ガラス１２ＷＳ～１２ＲＧの各々について求めた全天日射量を用いて、窓ガラス１２ＷＳ～１２ＲＧの各々における日射熱取得量の合計値である総日射熱取得量を求める。以下、具体的に説明する。
　基準面であるルーフ１１における全天日射量ＪＴＳは、次式（１）で表される。
Figure JPOXMLDOC01-appb-M000001
　ルーフ１１における直達日射量ＪＤは、次式（２）で求めることができる。
Figure JPOXMLDOC01-appb-M000002
　式（２）において、Ｉ０は太陽定数、Ｒは地球の動径、Ｐは直達日射の大気透過率、ｈは太陽高度角である。また、θは、ルーフ１１における直達日射の入射角であり、図５に示す面２０をルーフ１１に置き換えて考えればよい。太陽定数Ｉ０は、一例として理科年表等に記載された値を用いればよい。太陽定数Ｉ０は、値が更新される場合があるため、更新された場合は更新値を用いればよい。
　地球の動径Ｒ、太陽高度角ｈの計算には、例えば、文献（１）「（株）気象データシステム技術解説,2016. 太陽位置の計算、松本 真一著」に記載された計算方法を用いればよい。また、大気透過率Ｐは、一例として、文献（２）「直散分離と斜面日射量の計算,（株）気象データシステム技術解説,2017. 赤坂 裕 著」に記載された手法で求めればよいが、ここでは、一例として、適当な初期値に設定する。大気透過率Ｐの適当な初期値は、一例として０．７である。また、ｃｏｓθは、後述する式（５）で求めればよい。
　ルーフ１１における天空日射量ＪＳは、天空日射透過率ＳＲＴを用いて次式（３）で求めることができる。基準面としてのルーフ１１について天空日射量ＪＳを求める際には、天空日射透過率ＳＲＴは、適当な初期値に設定すればよい。天空日射透過率ＳＲＴの適当な初期値は、一例として、１．０である。
Figure JPOXMLDOC01-appb-M000003
　式（３）において、形態係数ＶＦＳは、ルーフ１１の全天空（天空２の全体）に対する形態係数ＶＦＳであり、一例として、ルーフ１１の傾斜角βを用いて後述する式（６）で求めればよい。また、太陽定数Ｉ０、地球の動径Ｒ、太陽高度角ｈの計算には、例えば、式（２）と同様の計算方法を用いればよい。
　ルーフ１１における地面反射日射量ＪＧは、次式（４）で求めることができる。
Figure JPOXMLDOC01-appb-M000004
　式（４）において、形態係数ＶＦＧは、ルーフ１１の地面１に対する形態係数ＶＦＧであり、ルーフ１１の傾斜角βを用いて後述する式（７）で求めることができる。また、式（４）において、ρは、地面反射率である。地面反射率ρは、一例として、文献（３）「最新・建築環境工学[改訂3版], 井上書院, 2006. 田中 俊六,武田 仁,足立 哲夫,土屋 喬雄,寺尾 道仁 著」に記載された手法で求めればよい。
　また、式（４）において、天空日射量ＪＳは式（３）で求まる値を用いればよい。形態係数ＶＦＳは、ルーフ１１の全天空（天空２の全体）に対する形態係数ＶＦＳであり、一例として、ルーフ１１の傾斜角βを用いて後述する式（６）で求めればよい。また、太陽定数Ｉ０、地球の動径Ｒ、太陽高度角ｈは、一例として、式（２）における値と同一値を用いればよい。
　ルーフ１１への直達日射の入射角θの余弦（ｃｏｓθ）は、太陽高度角ｈ、ルーフ１１の傾斜角β、太陽の方位角Ａ、ルーフ１１の方位角Ａ’を用いると、次式（５）で求めることができる。
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
　窓ガラス１２ＷＳ～１２ＲＧの各々について日射熱取得率ＴＴＳを求める際に、式（１６）における日射透過率ＴＤＳと、式（１７）における日射吸収率α、室外側熱伝達係数ｈｅ、及び室内側熱伝達係数ｈｉとの導出は、一例として、文献（４）「Road vehicles-Safety glazing materials-Method for the determination of solar transmittance、ISO13837,2008」に記載された手法に従って行えばよい。
　車両１０の車室内に入り込む熱量を表す総日射熱取得量ＣＱＴＳ（Ｗ）は、窓ガラス１２ＷＳ～１２ＲＧの各々について求めた日射熱取得量ＱＴＳに、窓ガラス１２ＷＳ～１２ＲＧの各々の面積Ｓをそれぞれ乗算した積の総和として、次式（１８）で表されるように算出される。
Figure JPOXMLDOC01-appb-M000018
　なお、窓ガラス１２ＷＳ～１２ＲＧの各々の日射透過率ＴＤＳ及び日射吸収率αが、窓ガラス１２ＷＳ～１２ＲＧの各々への日射の入射角度によって変化する入射角度特性を有する場合は、文献（５）「窓及びドアの熱性能-日射熱取得率の計算, JIS A 2103, 2014」に記載された計算方法に従って、入射角度特性を考慮して計算を行えばよい。
　図８は、実施の形態の日射推定システム１００を車両１０に搭載した構成の一例を示す図である。図８には、車両１０に搭載される構成要素として、ＨＶＡＣ(Heating Ventilation and Air Conditioning)－ＥＣＵ(Electronic Control Unit：電子制御装置)２００、パワートレインＥＣＵ３００Ａ、ナビゲーションＥＣＵ３００Ｂ、及びＥＣＵ３００Ｃ、通信部３００Ｄ、日射計１０１、外気温センサ４１０、室温センサ４２０、湿度センサ４３０、車速センサ４４０、ＧＰＳ(Global Positioning System)ユニット４５０、加速度センサ４６０、及びＨＶＡＣユニット４７０を示す。パワートレインＥＣＵ３００Ａは、電子制御装置の一例であり、ナビゲーションＥＣＵ３００Ｂは、ナビゲーション制御部の一例であり、ＥＣＵ３００Ｃは、電子制御装置の一例であり、通信部３００Ｄは、通信部の一例であり、受信部の一例であり、取得部の一例である。
　ここで、ＨＶＡＣ－ＥＣＵ２００、パワートレインＥＣＵ３００Ａ、ナビゲーションＥＣＵ３００Ｂ、及びＥＣＵ３００Ｃは、ＣＰＵ(Central Processing Unit)、ＲＡＭ(Random Access Memory)、ＲＯＭ(Read Only Memory)、入出力インターフェース、及び内部バス等を含むコンピュータによって実現される。また、通信部３００Ｄは、一例として、４Ｇ(Fourth Generation)、ＬＴＥ(Long Term Evolution)、５Ｇ(Fifth Generation)等の公衆通信回線を介して無線通信を行う通信機である。ＨＶＡＣ－ＥＣＵ２００、パワートレインＥＣＵ３００Ａ、ナビゲーションＥＣＵ３００Ｂ、ＥＣＵ３００Ｃ、及び通信部３００Ｄは、ＣＡＮ(Controller Area Network)バス５０によって互いに通信可能に接続されている。
　ＨＶＡＣ－ＥＣＵ２００は、車両１０のＨＶＡＣユニット４７０の制御を行うＥＣＵである。ＨＶＡＣユニット４７０は、空調装置の一例であり、ヒータ、ベンチレーション装置、及びエアコンディショニング装置等を含む。ＨＶＡＣ－ＥＣＵ２００には、日射計１０１、外気温センサ４１０、室温センサ４２０、及び湿度センサ４３０が接続されている。
　ＨＶＡＣ－ＥＣＵ２００は、主制御部１１０、算出部１２０、補正部１３０、推定部１４０、計算部１５０、メモリ１６０、制御部２１０、及びメモリ２２０を有する。主制御部１１０、算出部１２０、補正部１３０、推定部１４０、計算部１５０、及び制御部２１０は、ＨＶＡＣ－ＥＣＵ２００が実行するプログラムの機能（ファンクション）を機能ブロックとして示したものである。また、メモリ１６０及び２２０は、ＨＶＡＣ－ＥＣＵ２００のメモリを機能的に表したものである。
　ここで、主制御部１１０、算出部１２０、補正部１３０、推定部１４０、計算部１５０、メモリ１６０、及び日射計１０１は、日射推定システム１００を構築するため、これらの符号には括弧書きで符号１００を記す。実施の形態の日射推定方法は、日射計１０１、主制御部１１０、算出部１２０、補正部１３０、推定部１４０、及び計算部１５０によって行われる。より具体的には、実施の形態の日射推定方法は、日射計１０１による全天日射量の取得と、主制御部１１０、算出部１２０、補正部１３０、推定部１４０、及び計算部１５０によって実行される処理とによって実現される。
　また、日射推定システム１００の構成要素のうち、主制御部１１０、算出部１２０、補正部１３０、推定部１４０、計算部１５０、及びメモリ１６０は、実施の形態の日射推定プログラムを実現する機能ブロックであり、ここでは、一例として、ＨＶＡＣ－ＥＣＵ２００に含まれる。主制御部１１０、算出部１２０、補正部１３０、推定部１４０、計算部１５０、及びメモリ１６０がＨＶＡＣ－ＥＣＵ２００に含まれるのは、一例として、実施の形態の日射推定プログラムをＨＶＡＣ－ＥＣＵ２００にインストールしたからである。ただし、主制御部１１０、算出部１２０、補正部１３０、推定部１４０、計算部１５０、及びメモリ１６０は、ＨＶＡＣ－ＥＣＵ２００とは異なる独自のＥＣＵ（日射推定ＥＣＵ）にインストールされていてもよく、この場合は、ＣＡＮバス５０を介して日射推定ＥＣＵをＨＶＡＣ－ＥＣＵ２００、パワートレインＥＣＵ３００Ａ、ナビゲーションＥＣＵ３００Ｂ、ＥＣＵ３００Ｃ、及び通信部３００Ｄに接続すればよい。
　主制御部１１０、算出部１２０、補正部１３０、推定部１４０、及び計算部１５０は、実施の形態の日射推定プログラムを実現するための処理を行い、メモリ１６０には、実施の形態の日射推定プログラムを実現するために利用するデータが格納される。ここでは、主制御部１１０、算出部１２０、補正部１３０、推定部１４０、計算部１５０、及びメモリ１６０について説明する前に、車両１０に搭載される他の構成要素について先に説明する。
　制御部２１０は、空調制御部の一例であり、計算部１５０によって計算される総日射熱取得量ＣＱＴＳを用いて、ＨＶＡＣユニット４７０の制御を行う。メモリ２２０には、ＨＶＡＣユニット４７０による車両１０の室内の温度や湿度を調整するために利用するプログラムやデータ等が格納されている。ここで、日射推定システム１００、制御部２１０、及びＨＶＡＣユニット４７０は、実施の形態の空調制御システムを構築する。また、算出部１２０、補正部１３０、及び推定部１４０を有するＨＶＡＣ－ＥＣＵ２００と、ＨＶＡＣユニット４７０とは、実施の形態の空調装置を構築する。この空調装置のＨＶＡＣ－ＥＣＵ２００には、さらに計算部１５０が含まれてもよい。また、実施の形態の空調制御システム、及び、実施の形態の空調装置におけるＨＶＡＣユニット４７０の制御方法は、実施の形態の空調制御方法である。
　パワートレインＥＣＵ３００Ａは、車両１０のパワートレインの制御を行う。パワートレインは、車両１０がＥＶ(Electric Vehicle)である場合には、動力源であるモータ、モータの駆動に用いるインバータ、及びバッテリの充放電等の制御を行う。パワートレインＥＣＵ３００Ａは、車両１０がＨＶ(Hybrid Vehicle)又はプラグインＨＶである場合は、モータ、インバータ、及びバッテリの充放電等の制御に加えて、エンジンの点火タイミングや燃料噴射料等の制御を行う。また、パワートレインＥＣＵ３００Ａは、車両１０がエンジンを動力源とする場合には、エンジンの点火タイミングや燃料噴射料等の制御を行う。
　パワートレインＥＣＵ３００Ａには、車速センサ４４０が接続されている。パワートレインＥＣＵ３００Ａは、車速センサ４４０から取得する車速データをＣＡＮバス５０を介してＨＶＡＣ－ＥＣＵ２００及びナビゲーションＥＣＵ３００Ｂに伝送する。車速データは、ＥＣＵ３００Ｃ、及び通信部３００Ｄにも伝送されてもよい。
　ナビゲーションＥＣＵ３００Ｂには、ＧＰＳユニット４５０が接続されている。ナビゲーションＥＣＵ３００Ｂは、ＧＰＳユニット４５０から取得する位置データを用いて、経路探索、経路案内、及び到着時刻の予想等の処理を行う。位置データは、ナビゲーションＥＣＵ３００ＢからＣＡＮバス５０を介してＨＶＡＣ－ＥＣＵ２００に伝送される。なお、ナビゲーションＥＣＵ３００Ｂは、ナビゲーション制御部の一例であり、ＧＰＳユニット４５０は、位置データ取得部の一例である。
　ＥＣＵ３００Ｃは、車両１０に搭載される、パワートレインＥＣＵ３００Ａ及びナビゲーションＥＣＵ３００Ｂ以外のＥＣＵを示したものである。一例として、ＥＣＵ３００には、加速度センサ４６０が接続されており、ＥＣＵ３００が加速度センサから取得する加速度データは、ＣＡＮバス５０を介してＨＶＡＣ－ＥＣＵ２００に伝送される。加速度データは、車両１０の進行方向と、前後方向及び横方向の傾きとを表す。
　通信部３００Ｄは、インターネットを介して車両１０の外部のサーバ等から各種データをダウンロードするとともに、車両１０のデータをサーバにアップロードする。通信部３００Ｄは、受信したデータをＣＡＮバス５０に出力する。なお、通信部３００Ｄは、ビーコン信号等を受信する受信部であってもよい。
　外気温センサ４１０は、室温センサ４２０、及び湿度センサ４３０は、それぞれ、車両１０の外部の気温（外気温）、車両１０の室内の温度、及び車両１０の室内の湿度を表すデータを取得し、ＨＶＡＣ－ＥＣＵ２００に伝送する。車速センサ４４０及びＧＰＳユニット４５０は、それぞれ、車速及び位置（緯度、経度）を検出し、車速データ及び位置データをパワートレインＥＣＵ３００Ａ及びナビゲーションＥＣＵ３００Ｂに伝送する。なお、位置データには、位置データの取得時の日付及び時刻（年月日時分秒）を表す時刻データが含まれる。
　日射計１０１は、取得部の一例であり、一例として、車両１０の外面にある基準面であるルーフ１１（図２、３参照）に配置され、ルーフ１１における日射量を取得する。日射計１０１が取得する日射量は、ルーフ１１における直達日射、天空日射、及び地面反射日射をすべて含むため、ルーフ１１における全天日射量の実測値である。日射計１０１が取得する全天日射量の実測値は、第１全天日射量の一例である。なお、車両１０の外面とは、車両１０の外側にあって、日射を受ける面である。
　主制御部１１０は、日射推定プログラムの処理を統括する処理部であり、算出部１２０、補正部１３０、推定部１４０、及び計算部１５０が実行する処理以外の処理を実行する。例えば、主制御部１１０は、日射計１０１から全天日射量を取得し、メモリ１６０に格納する。
　算出部１２０は、日射計１０１による全天日射量（実測値）の取得時に、式（１）～（７）を用いてルーフ１１における全天日射量ＪＴＳを算出する。算出部１２０がルーフ１１について算出する全天日射量ＪＴＳは、第２全天日射量の一例である。算出部１２０が実測値の取得時におけるルーフ１１での全天日射量ＪＴＳを算出するのは、実測値の取得時と同一時刻における条件下で算出した算出値と実測値とを比較するためである。
　算出部１２０は、式（１）に従って、ルーフ１１における全天日射量ＪＴＳを以下のように算出する。
　算出部１２０は、式（２）を用いて直達日射量ＪＤを算出する際には、太陽定数Ｉ０、地球の動径Ｒ、大気透過率Ｐ、太陽高度角ｈ、及びｃｏｓθとして以下のデータを用いる。
　算出部１２０は、太陽定数Ｉ０として、一例として理科年表等に記載された値を用いるが、更新された場合は更新値を用いればよい。また、算出部１２０は、大気透過率については、大気透過率Ｐの適当な初期値を用いる。適当な初期値は、一例として、１．０である。
　また、算出部１２０は、地球の動径Ｒ及び太陽高度角ｈについては、日射計１０１による全天日射量の取得時（日付及び時刻（年月日時分秒））に応じた値を用いればよい。地球の動径Ｒと太陽高度角ｈを求める際に利用する時刻データは、ナビゲーションＥＣＵ３００Ｂから取得する位置データに含まれる時刻データ、又は、ＨＶＡＣ－ＥＣＵ２００が利用する時刻データを用いればよい。地球の動径Ｒ及び太陽高度角ｈと日付及び時刻（年月日時分秒）とを関連付けたテーブルデータをメモリ１６０に格納しておき、時刻データを用いて、地球の動径Ｒ及び太陽高度角ｈを求めればよい。
　また、算出部１２０は、ｃｏｓθについては、日射計１０１による全天日射量の取得時におけるルーフ１１についての値を用いる。ｃｏｓθは、日付及び時刻によって変化するからである。ｃｏｓθを求める際に利用する位置データ（時刻データを含む）については、ナビゲーションＥＣＵ３００Ｂから取得する位置データに含まれる時刻データ、又は、ＨＶＡＣ－ＥＣＵ２００が利用する時刻データを用いればよい。
　式（５）に基づいてｃｏｓθを求める際に利用する、太陽高度角ｈ、傾斜角β、太陽方位角Ａ、及び、方位角Ａ’のうち、太陽高度角ｈについては、上述のように求めた値を用いればよい。傾斜角β及び方位角Ａ’は、ルーフ１１についての値である。
　ルーフ１１の傾斜角βについては、加速度センサ４６０によって取得される加速度データと、車両１０のＸＹＺ座標系におけるルーフ１１の法線ベクトルを表すデータとに基づいて求めればよい。ルーフ１１の傾斜角βは、基準面の向きに関係する第１関係値の一例である。
　加速度データは、車両１０の進行方向と、前後方向及び横方向の傾きとを表す。傾斜角βを求める際に用いる車両１０のＸＹＺ座標系におけるルーフ１１の法線ベクトルを表すデータは、メモリ１６０に格納しておけばよい。
　太陽方位角Ａについては、日射計１０１による全天日射量の取得時（日付及び時刻（年月日時分秒））に応じた値を用いればよい。太陽方位角Ａを求める際に利用する時刻データは、地球の動径Ｒと太陽高度角ｈを求める際に用いる時刻データを用いればよい。太陽方位角Ａと日付及び時刻（年月日時分秒）とを関連付けたテーブルデータをメモリ１６０に格納しておき、時刻データを用いて、太陽方位角Ａを求めればよい。
　ルーフ１１の方位角Ａ’については、加速度センサ４６０によって取得される加速度データと、車両１０のＸＹＺ座標系におけるルーフ１１の法線ベクトルを表すデータとに基づいて求めればよい。
　算出部１２０は、以上のようにして求めた値を用いて、式（２）及び式（５）を用いてルーフ１１における直達日射量ＪＤを算出する。算出部１２０によって算出されたルーフ１１における直達日射量ＪＤは、第１直達日射量の一例である。
　また、算出部１２０は、式（３）に従ってルーフ１１について天空日射量ＪＳを算出する。このときに、太陽定数Ｉ０、地球の動径Ｒ、太陽高度角ｈ、大気透過率Ｐについては、直達日射量ＪＤを求める際の値を用いればよい。形態係数ＶＦＳについては、直達日射量ＪＤを求める際のルーフ１１の傾斜角βを用いて式（６）に従って求めればよい。また、天空日射透過率ＳＲＴは、適当な初期値に設定すればよい。天空日射透過率ＳＲＴの適当な初期値は、一例として、１．０である。算出部１２０が算出する天空日射量ＪＳは、第１天空日射量の一例である。
　また、算出部１２０は、式（４）に従ってルーフ１１についての地面反射日射量ＪＧを算出する。このときに、太陽定数Ｉ０、地球の動径Ｒ、太陽高度角ｈ、大気透過率Ｐ、形態係数ＶＦＳについては、天空日射量ＪＳを求める際の値を用いればよい。形態係数ＶＦＧについては、直達日射量ＪＤを求める際のルーフ１１の傾斜角βを用いて式（７）に従って求めればよい。また、地面反射率ρは、上述した文献（３）に記載された手法で求めればよい。算出部１２０が算出する地面反射日射量ＪＧは、第１地面反射日射量の一例である。
　算出部１２０は、上述のようにして算出した、直達日射量ＪＤ、天空日射量ＪＳ、及び地面反射日射量ＪＧを式（１）に従って加算することで、ルーフ１１における全天日射量ＪＴＳを算出する。算出部１２０がルーフ１１について算出する全天日射量ＪＴＳは、第２全天日射量の一例である。
　補正部１３０は、算出部１２０が算出した全天日射量ＪＴＳ（算出値）が、日射計１０１で取得される全天日射量（実測値）と略一致するように収束計算を行う。収束計算では、天空日射透過率ＳＲＴと大気透過率Ｐとを補正し、日射計１０１で得られる全天日射量と略一致する正確な全天日射量ＪＴＳを与える補正天空日射透過率ＣＳＲＴと補正大気透過率ＣＰとを求める。収束計算としては、一例として、二分法を用いる。また、略一致するとは、一例として、実測値に対する算出値の誤差が±１％以下であることをいう。なお、二分法以外の収束計算方法を用いてもよい。
　推定部１４０は、補正部１３０によって求められる補正天空日射透過率ＣＳＲＴ及び補正大気透過率ＣＰを用いて、式（８）～（１４）に従って窓ガラス１２ＷＳ～１２ＲＧの各々における全天日射量ＪＴＳを算出する。推定部１４０が窓ガラス１２ＷＳ～１２ＲＧの各々について算出する全天日射量ＪＴＳは、推定全天日射量の一例である。
　推定部１４０は、式（８）に従って、窓ガラス１２ＷＳ～１２ＲＧの各々における全天日射量ＪＴＳを以下のように算出する。
　推定部１４０は、式（９）を用いて直達日射量ＪＤを算出する際には、太陽定数Ｉ０、地球の動径Ｒ、補正大気透過率ＣＰ、太陽高度角ｈ、及びｃｏｓθとして以下のデータを用いる。推定部１４０が算出する直達日射量ＪＤは、第２直達日射量の一例であり、補正大気透過率ＣＰを用いて求めた推定値である。
　推定部１４０は、太陽定数Ｉ０としては、算出部１２０がルーフ１１における全天日射量ＪＴＳを算出する際に用いる値と同一の値を用いればよい。また、推定部１４０は、大気透過率については、補正部１３０によって求められる補正大気透過率ＣＰを用いる。
　また、推定部１４０は、地球の動径Ｒ及び太陽高度角ｈについては、算出部１２０がルーフ１１における全天日射量ＪＴＳを算出する際に用いる値と同一値を用いればよい。
　また、推定部１４０は、ｃｏｓθについては、日射計１０１による全天日射量の取得時における窓ガラス１２ＷＳ～１２ＲＧの各々についての固有の値を用いる。ｃｏｓθは、日時及び時刻によって変化し、窓ガラス１２ＷＳ～１２ＲＧは互いに向きが異なる（法線ベクトルの方向が異なる）からである。
　ｃｏｓθを求める際に利用する位置データ（時刻データを含む）については、ナビゲーションＥＣＵ３００Ｂから取得する位置データに含まれる時刻データ、又は、ＨＶＡＣ－ＥＣＵ２００が利用する時刻データを用いればよい。
　式（１２）に基づいてｃｏｓθを求める際に利用する、太陽高度角ｈ、傾斜角β、太陽方位角Ａ、及び、方位角Ａ’のうち、太陽高度角ｈについては、算出部１２０がルーフ１１における全天日射量ＪＴＳを算出する際に用いる値と同一値を用いればよい。傾斜角β及び方位角Ａ’は、窓ガラス１２ＷＳ～１２ＲＧの各々について異なる。
　窓ガラス１２ＷＳ～１２ＲＧの各々の傾斜角βについては、加速度センサ４６０によって取得される加速度データと、車両１０のＸＹＺ座標系における窓ガラス１２ＷＳ～１２ＲＧの各々の法線ベクトルを表すデータとに基づいて求めればよい。窓ガラス１２ＷＳ～１２ＲＧの傾斜角βは、対象物の外面のうちの所定外面の向きに関係する第２関係値の一例である。
　車両１０のＸＹＺ座標系における窓ガラス１２ＷＳ～１２ＲＧの各々の法線ベクトルを表すデータは、メモリ１６０に格納しておけばよい。太陽方位角Ａについては、算出部１２０がルーフ１１における全天日射量ＪＴＳを算出する際に用いる値と同一の値を用いればよい。
　窓ガラス１２ＷＳ～１２ＲＧの各々の方位角Ａ’については、加速度センサ４６０によって取得される加速度データと、車両１０のＸＹＺ座標系における窓ガラス１２ＷＳ～１２ＲＧの各々の法線ベクトルを表すデータとに基づいて求めればよい。
　推定部１４０は、以上のようにして求めた値を用いて、式（９）及び式（１２）を用いて窓ガラス１２ＷＳ～１２ＲＧの各々についての直達日射量ＪＤを算出する。
　また、推定部１４０は、式（１０）に従って窓ガラス１２ＷＳ～１２ＲＧの各々について天空日射量ＪＳを算出する。推定部１４０が算出する天空日射量ＪＳは、第２天空日射量の一例であり、補正大気透過率ＣＰと補正天空日射透過率ＣＳＲＴを用いて求めた推定値である。
　推定部１４０は、太陽定数Ｉ０、地球の動径Ｒ、太陽高度角ｈについては、窓ガラス１２ＷＳ～１２ＲＧの各々についての直達日射量ＪＤを求める際の値を用いればよい。太陽定数Ｉ０、地球の動径Ｒ、太陽高度角ｈは、窓ガラス１２ＷＳ～１２ＲＧについて共通の値である。また、大気透過率については、補正大気透過率ＣＰを用いればよい。
　形態係数ＶＦＳについては、窓ガラス１２ＷＳ～１２ＲＧの各々について直達日射量ＪＤを求める際の窓ガラス１２ＷＳ～１２ＲＧの各々の傾斜角βをそれぞれ用いて式（１３）に従って求めればよい。また、補正部１３０によって補正された補正天空日射透過率ＣＳＲＴを用いればよい。
　また、推定部１４０は、式（１１）に従って窓ガラス１２ＷＳ～１２ＲＧの各々についての地面反射日射量ＪＧを算出する。このときに、太陽定数Ｉ０、地球の動径Ｒ、太陽高度角ｈ、補正大気透過率ＣＰ、形態係数ＶＦＳについては、窓ガラス１２ＷＳ～１２ＲＧの各々についての天空日射量ＪＳを求める際の値をそれぞれ用いればよい。太陽定数Ｉ０、地球の動径Ｒ、太陽高度角ｈ、補正大気透過率ＣＰについては、窓ガラス１２ＷＳ～１２ＲＧについて共通の値であり、形態係数ＶＦＳについては、窓ガラス１２ＷＳ～１２ＲＧの各々について別々の値である。窓ガラス１２ＷＳ～１２ＲＧの各々で傾斜角βが異なるからである。
　また、形態係数ＶＦＧについては、窓ガラス１２ＷＳ～１２ＲＧの各々についての傾斜角βを用いて式（１４）に従って求めればよい。また、地面反射率ρは、算出部１２０がルーフ１１についての地面反射日射量ＪＧを算出する際の値を用いればよい。推定部１４０が算出する地面反射日射量ＪＧは、第２地面反射日射量の一例である。
　推定部１４０は、窓ガラス１２ＷＳ～１２ＲＧの各々について、上述のようにして算出した、直達日射量ＪＤ、天空日射量ＪＳ、及び地面反射日射量ＪＧを式（８）に従って加算することで、全天日射量ＪＴＳを算出する。推定部１４０が窓ガラス１２ＷＳ～１２ＲＧの各々について算出する全天日射量ＪＴＳは、推定全天日射量の一例であり、補正大気透過率ＣＰと補正天空日射透過率ＣＳＲＴを用いて求めた推定値である。推定部１４０は、窓ガラス１２ＷＳ～１２ＲＧについての６つの全天日射量ＪＴＳを算出する。
　計算部１５０は、式（１５）に従って窓ガラス１２ＷＳ～１２ＲＧの各々の日射熱取得量ＱＴＳを計算する際に、次のような計算処理を行う。まず、計算部１５０は、式（１６）に従って窓ガラス１２ＷＳ～１２ＲＧの各々についての日射熱取得率ＴＴＳを求める。
　日射熱取得率ＴＴＳを求める際には、日射透過率ＴＤＳと、式（１７）で表される係数ｑｉとを用いる。係数ｑｉを求めるには、日射吸収率α、室外側熱伝達係数ｈｅ、及び室内側熱伝達係数ｈｉの値を用いる。日射透過率ＴＤＳと日射吸収率αは、窓ガラス１２ＷＳ～１２ＲＧの各々についての値が決まっているため、メモリ１６０に格納しておけばよい。
　室内側熱伝達係数ｈｉと室外側熱伝達係数ｈｅは、窓ガラス１２ＷＳ～１２ＲＧについて、文献（４）に記載された計算方法で求めればよい。なお、窓ガラス１２ＷＳ～１２ＲＧの各々のガラス温度は、外気温センサ４１０によって取得される外気温と、室温センサ４２０によって取得される室温とに基づいて求めることができる。
　計算部１５０は、メモリ１６０から読み出した日射吸収率αと、上述のようにして求まる室内側熱伝達係数ｈｉ及び室外側熱伝達係数ｈｅとを用いて、式（１７）に従って係数ｑｉを求める。そして、計算部１５０は、メモリ１６０から読み出した日射透過率ＴＤＳと、係数ｑｉとを用いて、式（１６）に従って日射熱取得率ＴＴＳを求める。
　計算部１５０は、式（１５）に従って、窓ガラス１２ＷＳ～１２ＲＧの各々について求めた全天日射量ＪＴＳに、窓ガラス１２ＷＳ～１２ＲＧの各々の日射熱取得率ＴＴＳを乗じて、窓ガラス１２ＷＳ～１２ＲＧの各々の日射熱取得量ＱＴＳを算出する。以上により、窓ガラス１２ＷＳ～１２ＲＧの各々についての日射熱取得量ＱＴＳが求まる。
　計算部１５０は、式（１８）に従って、窓ガラス１２ＷＳ～１２ＲＧの各々について求めた日射熱取得量ＱＴＳに、窓ガラス１２ＷＳ～１２ＲＧの各々の面積Ｓをそれぞれ乗算した積の総和を求めることによって、車両１０の室内の総日射熱取得量ＣＱＴＳを計算する。窓ガラス１２ＷＳ～１２ＲＧの各々の面積Ｓを表すデータは、メモリ１６０に格納しておけばよい。以上により、日射によって車両１０の室内に入り込む総日射熱取得量ＣＱＴＳが求まる。
　図９は、日射推定プログラムによって実現される処理を示すフローチャートである。日射推定プログラムは、図９に示すスタートからエンドまでの処理を一例として１秒おきに繰り返し実行する。
　主制御部１１０は、処理がスタートすると、現在時刻（日付及び時刻（年月日時分秒））、車両１０の位置（緯度、経度）、車両１０の進行方向の方位角、車両１０の傾き角度を求める（ステップＳ１）。現在時刻と車両１０の位置（緯度、経度）は、ＧＰＳユニット４５０で取得される位置データ（時刻データを含む）から取得すればよい。なお、現在時刻は、ＨＶＡＣ－ＥＣＵ２００が利用する時刻データを用いてもよい。車両１０の進行方向の方位角は、位置データと、加速度センサ４６０によって取得される加速度データとに基づいて求めればよい。車両１０の傾き角度については、加速度データが表す前後方向の傾き角度及び横方向の傾き角度を取得すればよい。
　主制御部１１０は、メモリ１６０からルーフ１１及び窓ガラス１２ＷＳ～１２ＲＧの各々についての法線ベクトルのデータを読み出し、車両１０の位置、車両１０の方位角、車両１０の傾き角度を用いて、すべての法線ベクトルの座標を絶対座標系における法線ベクトルの座標に変換する（ステップＳ２）。ここで、絶対座標系とは、地球上における位置と方向を表す座標系であり、一例として、ＧＰＳで用いる測位系の座標を用いることができる。
　主制御部１１０は、メモリ１６０に格納されている太陽方位角Ａと太陽高度角ｈの計算用のテーブルデータを用いて、現在時刻における太陽方位角Ａと太陽高度角ｈを求める（ステップＳ３）。
　算出部１２０は、式（１）に従って、ルーフ１１における全天日射量ＪＴＳを算出する（ステップＳ４）。より具体的には、式（２）～（７）に従って、ルーフ１１における直達日射量ＪＤ、天空日射量ＪＳ、地面反射日射量ＪＧを求め、式（１）に従って加算することで全天日射量ＪＴＳを算出する。
　補正部１３０は、収束計算を行って天空日射透過率ＳＲＴと大気透過率Ｐとを補正し、日射計１０１で得られる全天日射量と略一致する正確な全天日射量ＪＴＳを与える補正天空日射透過率ＣＳＲＴと補正大気透過率ＣＰとを求める（ステップＳ５）。
　推定部１４０は、補正天空日射透過率ＣＳＲＴ及び補正大気透過率ＣＰを用いて、式（８）～（１４）に従って窓ガラス１２ＷＳ～１２ＲＧの各々における全天日射量ＪＴＳを算出する（ステップＳ６）。より具体的には、式（９）～（１４）に従って、窓ガラス１２ＷＳ～１２ＲＧの各々における直達日射量ＪＤ、天空日射量ＪＳ、地面反射日射量ＪＧを求め、式（８）に従って窓ガラス１２ＷＳ～１２ＲＧの各々について別々に加算することで、窓ガラス１２ＷＳ～１２ＲＧの各々における全天日射量ＪＴＳを算出する。
　計算部１５０は、式（１５）に従って窓ガラス１２ＷＳ～１２ＲＧの各々の日射熱取得量ＱＴＳを計算する（ステップＳ７）。日射熱取得率ＴＴＳは、式（１６）に従って求めればよい。
　計算部１５０は、式（１８）に従って、車両１０の室内における総日射熱取得量ＣＱＴＳを計算する（ステップＳ８）。計算部１５０は、計算した総日射熱取得量ＣＱＴＳをメモリ１６０に格納する。
　以上で日射推定プログラムによる一連の処理が終了する（エンド）。主制御部１１０は、メモリ１６０から総日射熱取得量ＣＱＴＳを読み出し、制御部２１０に伝送する。制御部２１０は、総日射熱取得量ＣＱＴＳをメモリ２２０に格納し、ＨＶＡＣユニット４７０の制御に用いる。
　図１０乃至図２３は、日射推定システム１００での実験結果を示す図である。ここでは、日本国内のある場所に車両１０を駐車し、全天日射量等の計算値と測定値との比較結果について説明する。車両１０は、真南を向くように水平面に配置した。すなわち、図２及び図３に示す車両１０のＸＹＺ座標系の－Ｘ方向が真南を向き、ＸＹ平面が水平面と平行になるように配置した。
　図１０乃至図２０及び図２３において、横軸は、２０１９年のある日における午前から午後にかけての時刻を２４時間表示で示し、縦軸は、全天日射量（図１０、１２～１４、１８、２０、２３）、日射量（図１５～１７）、又は、補正大気透過率ＣＰ及び補正天空日射透過率ＣＳＲＴ（図１１、１９）を示す。日射推定システム１００による全天日射量ＪＴＳの計算は、一例として、２秒毎に行った。
　図１０には、２０１９年１０月３１日において、日射推定システム１００でルーフ１１について求めた全天日射量ＪＴＳ（計算値）と、日射計１０１で求めた全天日射量の実測値とを示す。計算値と実測値は完全に一致しており、実測値と等しい全天日射量ＪＴＳの計算値が得られることを確認できた。当日の天気は晴れであり、雲が発生しなければ全天日射量の計算値（ＪＴＳ）と測定値は連続的な値を取ると考えられるが、実際には時々雲が発生するため、時々刻々と変化する特性を示した。
　図１１には、図１０に示す全天日射量の計算値（ＪＴＳ）と測定値を取得した際の補正大気透過率ＣＰ及び補正天空日射透過率ＣＳＲＴを示す。例えば、１０時５１分０６秒のあたりで補正大気透過率ＣＰ及び補正天空日射透過率ＣＳＲＴがともに低下している。特に、補正大気透過率ＣＰはゼロになっている。これは、太陽が雲に隠れて直達日射がゼロになったことを示している。図１１では、日射推定システム１００が求める補正大気透過率ＣＰ及び補正天空日射透過率ＣＳＲＴが時々刻々と変化していることが分かる。
　図１２には、２０１９年１０月３１日において、日射推定システム１００で窓ガラス１２ＷＳについて求めた全天日射量ＪＴＳ（計算値）と、日射計１０１とは別に窓ガラス１２ＷＳに配置した日射計で求めた全天日射量の実測値とを示す。図１２に示す計算値は、日射推定システム１００が日射計１０１の測定値を用いて窓ガラス１２ＷＳについて求めた全天日射量ＪＴＳ（計算値）である。
　日射計の取付位置に誤差があったため、全天日射量ＪＴＳ（計算値）と、日射計で求めた全天日射量の実測値とに少し差があるが、窓ガラス１２ＷＳについて全天日射量ＪＴＳ（計算値）を求めることができていることが分かる。
　図１３は、２０１９年１０月３１日において、日射推定システム１００で窓ガラス１２ＦＲについて求めた全天日射量ＪＴＳ（計算値）と、日射計１０１とは別に窓ガラス１２ＦＲに配置した日射計で求めた全天日射量の実測値とを示す。図１３に示す計算値は、日射推定システム１００が日射計１０１の測定値を用いて窓ガラス１２ＦＲについて求めた全天日射量ＪＴＳ（計算値）である。
　窓ガラス１２ＦＲは右前のサイドガラスであり、実験時には午前中には日射があたらず、午後になって太陽が西側に移動してから日射があたるため、午前中は全天日射量の計算値ＪＴＳと実測値が低く、午後になってから増大している。午前には窓ガラス１２ＦＲに日射が当たらないため、午前に得られている全天日射量は、天空日射による日射量が主体的であると考えられる。このように、窓ガラス１２ＦＲについて全天日射量ＪＴＳ（計算値）を求めることができていることが分かる。
　図１４は、２０１９年１０月３１日において、日射推定システム１００で窓ガラス１２ＦＬについて求めた全天日射量ＪＴＳ（計算値）と、日射計１０１とは別に窓ガラス１２ＦＬに配置した日射計で求めた全天日射量の実測値とを示す。図１３に示す計算値は、日射推定システム１００が日射計１０１の測定値を用いて窓ガラス１２ＦＬについて求めた全天日射量ＪＴＳ（計算値）である。
　窓ガラス１２ＦＬは左前のサイドガラスであり、実験時には午前中には日射（朝日）があたり、午後になると陰になって日射が殆どあたらなくなるため、午前中は全天日射量の計算値ＪＴＳと実測値が高く、昼に向かった徐々に低下し、午後は低い値になっている。午後には窓ガラス１２ＦＬに日射が当たらないため、午後に得られている全天日射量は、天空日射による日射量が主体的であると考えられる。このように、窓ガラス１２ＦＬについて全天日射量ＪＴＳ（計算値）を求めることができていることが分かる。
　図１５には、日射量として、２０１９年１０月３１日において、日射推定システム１００で窓ガラス１２ＷＳについて求めた直達日射量ＪＤ（計算値）と、日射計１０１とは別に窓ガラス１２ＷＳに配置した日射計で求めた全天日射量の実測値とを示す。図１５に示す実測値は、図１２に示す実測値と同一である。図１５では、直達日射量ＪＤ（計算値）と全天日射量の実測値とを比較している。
　天空日射量ＪＳと地面反射日射量ＪＧとの分がなく、直達日射量ＪＤ（計算値）のみであるため、全天日射量の実測値とは差がある。１１時５５分０６秒から１２時２０分４２秒までの区間のように、直達日射量ＪＤ（計算値）が殆どゼロになっているときに、全天日射量の実測値は約２００Ｗ以上ある。これは、天空日射量ＪＳが約２００Ｗ／ｍ２から約３００Ｗ／ｍ２あることを示している。
　図１６には、日射量として、２０１９年１０月３１日において、日射推定システム１００で窓ガラス１２ＦＲについて求めた直達日射量ＪＤ（計算値）と、日射計１０１とは別に窓ガラス１２ＦＲに配置した日射計で求めた全天日射量の実測値とを示す。図１６に示す実測値は、図１３に示す実測値と同一である。図１６では、直達日射量ＪＤ（計算値）と全天日射量の実測値とを比較している。
　天空日射量ＪＳと地面反射日射量ＪＧとの分がなく、直達日射量ＪＤ（計算値）のみであるため、全天日射量の実測値とは差がある。窓ガラス１２ＦＲに殆ど日射が当たらない午前中は、直達日射量ＪＤ（計算値）が略ゼロであり、全天日射量の実測値は約１００Ｗ程度である。また、１１時４８分４２秒から１３時４３分５４秒までの区間のように、直達日射量ＪＤ（計算値）が殆どゼロになっているときに、全天日射量の実測値は約２００Ｗ以上ある。このように、天空日射量ＪＳが約１００Ｗ／ｍ２から約２００Ｗ／ｍ２あることが分かり、全天日射量ＪＴＳを求める際に、天空日射量ＪＳを考慮する重要性が分かる。
　図１７には、日射量として、２０１９年１０月３１日において、日射推定システム１００で窓ガラス１２ＦＬについて求めた直達日射量ＪＤ（計算値）と、日射計１０１とは別に窓ガラス１２ＦＬに配置した日射計で求めた全天日射量の実測値とを示す。図１７に示す実測値は、図１４に示す実測値と同一である。図１７では、直達日射量ＪＤ（計算値）と全天日射量の実測値とを比較している。
　天空日射量ＪＳと地面反射日射量ＪＧとの分がなく、直達日射量ＪＤ（計算値）のみであるため、全天日射量の実測値とは差がある。午前中には日射（朝日）があたり、午後になると陰になって日射が殆どあたらなくなるため、午前中は直達日射量ＪＤ（計算値）と全天日射量の実測値が高く、昼に向かった徐々に低下し、午後は直達日射量ＪＤ（計算値）は殆どゼロになり、全天日射量の計算値は低い値（約１００Ｗ／ｍ２～２００Ｗ／ｍ２）になっている。午後には窓ガラス１２ＦＬに日射が当たらないため、午後に得られている全天日射量は、天空日射による日射量が主体的であると考えられる。このように、天空日射量ＪＳが約１００Ｗ／ｍ２から約２００Ｗ／ｍ２あることが分かり、全天日射量ＪＴＳを求める際に、天空日射量ＪＳを考慮する重要性が分かる。
　図１８には、２０１９年８月２９日において、日射推定システム１００でルーフ１１について求めた全天日射量ＪＴＳ（計算値）と、日射計１０１で求めた全天日射量の実測値とを示す。図１８に示す計算値は、日射推定システム１００が日射計１０１の測定値を用いてルーフ１１について求めた全天日射量ＪＴＳ（計算値）である。なお、８月２９日の天気は曇りで時々雨が降る天気であった。
　全天日射量ＪＴＳ（計算値）と、日射計で求めた全天日射量の実測値とが殆ど一致しており、ルーフ１１について全天日射量ＪＴＳ（計算値）を求めることができていることが分かる。
　図１９は、図１８に示す全天日射量の計算値（ＪＴＳ）と測定値を取得した際の補正大気透過率ＣＰ及び補正天空日射透過率ＣＳＲＴを示す。計測開始から、例えば、１２時２７分０６秒のあたりまでは、補正大気透過率ＣＰはゼロであり、補正天空日射透過率ＣＳＲＴは、約０．５から約１．０の間で推移している。このときは、直達日射が無い状態であった。また、１２時２７分０６秒よりも後では、直達日射がある時間帯があるため、補正大気透過率ＣＰの値が部分的に増大するときがあり、補正天空日射透過率ＣＳＲＴの値は、１２時２７分０６秒よりも前のときの値に比べて高くなった。図１９では、日射推定システム１００が求める補正大気透過率ＣＰ及び補正天空日射透過率ＣＳＲＴが時々刻々と変化していることが分かる。
　図２０には、２０１９年８月２９日において、日射推定システム１００で車両１０のダッシュボードの上面について求めた全天日射量ＪＴＳ（計算値）と、ダッシュボードに配置した日射計１０１で求めた全天日射量の実測値とを示す。ダッシュボードの上面について求めた全天日射量ＪＴＳ（計算値）は、ルーフ１１を基準面として算出部１２０がルーフ１１における全天日射量ＪＴＳ（計算値）を求めてから、補正部１３０が補正大気透過率ＣＰと補正天空日射透過率ＣＳＲＴを求め、推定部１４０が窓ガラス１２ＷＳについて求めた全天日射量ＪＴＳに、窓ガラス１２ＷＳの透過率を乗算することによって求めた値である。
　ダッシュボードの上面について求めた全天日射量ＪＴＳ（計算値）は、ダッシュボードに配置した日射計１０１で求めた全天日射量の実測値と殆ど一致しており、車両１０の室内に日射計１０１を配置しても、ルーフ１１に配置したときと同様に計算できることが確認できた。
　図２１は、２０１９年１０月３１日において窓ガラス１２ＷＳ～１２ＲＧの各々について求めた全天日射量ＪＴＳ、直達日射量ＪＤ、天空日射量ＪＳ、地面反射日射量ＪＧ、及び、日射熱取得量ＱＴＳ×Ｓ（Ｗ）を示す。日射熱取得量ＱＴＳ×Ｓは、窓ガラス１２ＷＳ～１２ＲＧの日射熱取得量ＱＴＳに窓ガラス１２ＷＳ～１２ＲＧの各々の面積Ｓを乗算した値である。
　また、図２１には、車両１０の室内における全天日射量ＪＴＳ、直達日射量ＪＤ、天空日射量ＪＳ、地面反射日射量ＪＧ、及び総日射熱取得量ＣＱＴＳも示す。車両１０の室内について求めた値は、窓ガラス１２ＷＳ～１２ＲＧの各々について求めた値の合計値である。図２１に示す各値は、２０１９年１０月３１日において全天日射量ＪＴＳ（計算値）が最も大きい時点（１２時２７分０４秒）の値である。
　図２１に示すように、全天日射量ＪＴＳ、直達日射量ＪＤ、天空日射量ＪＳは、窓ガラス１２ＷＳが最も大きく、車両１０の右側に位置する窓ガラス１２ＦＲ、１２ＲＲが２番目、３番目に大きく、窓ガラス１２ＦＬ、１２ＲＬ、１２ＲＧが４番目から６番目のいずれかになる傾向を示した。また、地面反射日射量ＪＧは、車両１０のボンネットに近い窓ガラス１２ＷＳが最も小さく、その他の窓ガラス１２ＦＲ、１２ＦＬ、１２ＲＲ、１２ＲＬ、１２ＲＧの方が大きい値であった。
　日射熱取得量ＱＴＳ×Ｓは、窓ガラス１２ＷＳが最も大きく、車両１０の右側に位置する窓ガラス１２ＦＲ、１２ＲＲが２番目、３番目に大きく、窓ガラス１２ＦＬ、１２ＲＧ、１２ＲＬが４番目から６番目であった。総日射熱取得量ＣＱＴＳは、９３０．４Ｗであった。以上のように、車両１０の室内に入り込む熱量である総日射熱取得量ＣＱＴＳを計算できることが確認できた。
　図２２は、２０１９年１０月３１日において窓ガラス１２ＷＳ～１２ＲＧの各々について求めた全天日射量ＪＴＳ、直達日射量ＪＤ、天空日射量ＪＳ、地面反射日射量ＪＧ、及び、日射熱取得量ＱＴＳ×Ｓを示す。
　また、図２２には、車両１０の室内における全天日射量ＪＴＳ、直達日射量ＪＤ、天空日射量ＪＳ、地面反射日射量ＪＧ、及び総日射熱取得量ＣＱＴＳも示す。車両１０の室内について求めた値は、窓ガラス１２ＷＳ～１２ＲＧの各々について求めた値の合計値である。図２２に示す各値は、２０１９年１０月３１日において直達日射量ＪＤが０．０Ｗの時点（１２時０８分０４秒）の値である。
　直達日射量ＪＤが０．０Ｗでも、天空日射量ＪＳと地面反射日射量ＪＧは、それぞれ、５６．７Ｗ～３３２．２Ｗ、１．９Ｗ～４．４Ｗの値が得られており、全天日射量ＪＴＳにおける天空日射量ＪＳと地面反射日射量ＪＧを考慮する意味があることが確認できた。また、直達日射が無いため、左右の前の窓ガラス１２ＦＲ、１２ＦＬの全天日射量ＪＴＳ、天空日射量ＪＳ、地面反射日射量ＪＧの値は略同一であり、左右の後の窓ガラス１２ＲＲ、１２ＲＬの全天日射量ＪＴＳ、天空日射量ＪＳ、地面反射日射量ＪＧの値は同一であった。総日射熱取得量ＣＱＴＳは、４０９．８Ｗであった。以上のように、車両１０の室内に入り込む熱量である総日射熱取得量ＣＱＴＳを計算できることが確認できた。
　図２３は、２０１９年１０月３１日において、日射推定システム１００で窓ガラス１２ＷＳについて求めた全天日射量ＪＴＳ（計算値）と、補正大気透過率ＣＰの代わりに、固定値（０．７４）の大気透過率を用いて窓ガラス１２ＷＳについて計算した全天日射量ＪＴＳ（０．７４）とを示す。図２３に示す全天日射量ＪＴＳ（計算値）は、図１２に示す値と同一である。大気透過率の０．７４という値は、２０１９年１０月において、理科年表に掲載されている値であり、過去３０年の平均値である。このような大気透過率は、主に晴れの日に測定した値の平均値であると考えられる。
　図２３に示すように、全天日射量ＪＴＳ（０．７４）は、測定開始の０９時２１分３０秒から緩やかに増大し、１２時前後で最大値を取り、測定終了の１３時５６分４２秒に向けて緩やかに低下する特性を示した。これに対して、日射推定システム１００は、天空日射量ＪＳを考慮するとともに、補正大気透過率ＣＰ及び補正天空日射透過率ＣＳＲＴを２秒毎の計算の度に求めることにより、日射計１０１の測定値と略同一の全天日射量を求めることができ、車両１０の室内に入り込む熱量である総日射熱取得量ＣＱＴＳを正確に求めることができることが分かった。
　以上のように、日射推定システム１００は、車両１０の内部に入る全天日射量ＪＴＳを高精度に求めることができる。全天日射量ＪＴＳは、補正大気透過率ＣＰ及び補正天空日射透過率ＣＳＲＴを用いて求めた推定値である。
　したがって、車両１０の内部に入る日射の全天日射量を高精度に推定できる日射推定システム１００、日射推定方法、及び日射推定プログラムを提供できる。
　また、１つの日射計１０１を用いて、車両１０の室内に入ってくる総日射熱取得量を高精度に算出できるため、ＨＶＡＣユニット４７０の制御を非常に精緻に行うことができる。
　また、車両１０の外気温と室温を用いてＨＶＡＣユニット４７０を制御する場合には、外気温と室温の上昇又は低下を検出してからＨＶＡＣユニット４７０を作動させる。このときに、窓ガラス１２ＷＳ～１２ＲＧや室内の空気の熱容量の影響で、車両１０の室内に日射が当たってから室温が上昇するまでには遅れが生じる。このため、車両１０の外気温と室温を用いてＨＶＡＣユニット４７０を制御する場合には、日射が当たってから５分程度の遅れが生じる。
　これに対して、日射推定システム１００を用いれば、日射よる車両１０の室内の総日射熱取得量を瞬時に正確に推定できるため、ＨＶＡＣユニット４７０を早期に稼働させることができ、遅れを抑制でき、室温を素早く適温に制御できる。また、ＨＶＡＣユニット４７０を早期に稼働させて室温の上昇度合が少ないうちに室温を制御できるので、ＨＶＡＣユニット４７０の稼働量を低減でき、車両１０の燃費や航続距離の向上に貢献できる。また、車両１０がＥＶである場合には、地面反射日射がバッテリの寿命に与える影響をより正確に把握できる。
　なお、以上では、直達日射量ＪＤと、天空日射量ＪＳと、地面反射日射量ＪＧとを加算して全天日射量ＪＴＳを求める形態について説明したが、３つの日射量のうち、地面反射日射量ＪＧの割合は、約数％から約１０％程度であるため、地面反射日射量ＪＧを求めずに、直達日射量ＪＤと、天空日射量ＪＳとを加算して全天日射量ＪＴＳを求めてもよい。
　また、以上では、日射計１０１で基準面の全天日射量を取得する形態について説明したが、日射計１０１から全天日射量を取得する代わりに、通信部３００Ｄを介して全天日射量を表すデータを取得してもよい。例えば、気象台が全天日射量を定期的に（例えば１時間おきに）算出する場合には、通信部３００Ｄを介して気象台にアクセスし、気象台から取得し、基準面の全天日射量として用いてもよい。また、この場合に、気象台が算出する全天日射量に所定の係数等を乗算することによって、基準面の全天日射量に変換してもよい。この場合には、日射計１０１が不要になるので、システム構成を簡略化できる。また、気象台が算出する全天日射量を用いることにより、計算の精度を向上させることができる。なお、この場合に、通信部３００Ｄは、受信部の一例であり、取得部の一例である。
　また、以上では、透過壁体の一例が窓ガラス１２ＷＳ～１２ＲＧである形態について説明したが、透過壁体はガラス製に限らず、例えばポリカーボネート等の樹脂製であってもよい。
　また、以上では、基準面がルーフ１１であり、ルーフ１１に日射計１０１を配置し、算出部１２０がルーフ１１における全天日射量ＪＳＴを算出する形態について説明した。しかしながら、基準面は、車両１０の外面であれば、ルーフ１１以外のどこであってもよい。算出部１２０は、車両１０の基準面における全天日射量ＪＳＴを算出すればよい。
　また、基準面は車両１０の外面であるが、日射計１０１が車両１０の室内に設けられ、窓ガラス１２ＷＳ～１２ＲＧのいずれかを透過した日射の全天日射量を測定する形態であってもよい。例えば、日射計１０１が車両１０のダッシュボードの上に設けられる場合には、日射計１０１は、窓ガラス１２ＷＳを透過した日射の全天日射量を測定する。このような場合には、基準面を窓ガラス１２ＷＳの外表面に設定し、車両１０のダッシュボードの上に設けられた日射計１０１によって測定された全天日射量を窓ガラス１２ＷＳの光学特性に基づいて、窓ガラス１２ＷＳの外表面において測定された全天日射量に変換する変換部を用いればよい。
　図２４は、実施の形態の変形例の日射推定システム１００Ａを車両１０に搭載した構成の一例を示す図である。日射推定システム１００Ａは、図８に示す日射推定システム１００に対して、変換部１２５を追加した構成を有する。日射計１０１及び変換部１２５は、取得部の一例である。また、メモリ１６０には、日射計１０１によって測定された全天日射量を窓ガラス１２ＷＳの外表面において測定された全天日射量に変換する変換率を表す変換率データを格納しておけばよい。変換率は、光学特性の一例である。
　変換率データは、一例として、窓ガラス１２ＷＳの透過率の逆数で与えられる変換率を表すデータであればよい。窓ガラス１２ＷＳの透過率は、窓ガラス１２ＷＳの外表面における日射量に対する透過量の割合を表し、図７に示す窓ガラス３０の外表面における日射量に対する透過量の割合と同様である。
　変換部１２５は、日射計１０１によって測定される全天日射量に、変換率データが表す変換率を乗算した全天日射量を求め、補正部１３０は、変換部１２５によって変換率が乗算された全天日射量と、算出部１２０によって算出される全天日射量ＪＳＴとを比較すればよい。また、この場合には、算出部１２０は、式（１）に従って、基準面としての窓ガラス１２ＷＳの外表面における全天日射量ＪＴＳを算出すればよい。推定部１４０は、窓ガラス１２ＦＲ、１２ＲＲ、１２ＦＬ、１２ＲＬ、１２ＲＧの各々についての全天日射量ＪＴＳを算出し、計算部１５０は、窓ガラス１２ＦＲ、１２ＲＲ、１２ＦＬ、１２ＲＬ、１２ＲＧの各々についての全天日射量ＪＴＳに、それぞれの面積Ｓを乗算した積の合計値を５つの窓ガラス１２ＦＲ、１２ＲＲ、１２ＦＬ、１２ＲＬ、１２ＲＧについての総日射熱取得量ＣＱＴＳとして求めればよい。車両１０の室内に入るすべての総日射熱取得量を計算する場合は、５つの窓ガラス１２ＦＲ、１２ＲＲ、１２ＦＬ、１２ＲＬ、１２ＲＧについての総日射熱取得量ＣＱＴＳと、窓ガラス１２ＷＳの全天日射量ＪＴＳに面積Ｓを乗算した積とを加算すればよい。
　また、以上では、日射推定システム１００を車両１０に搭載する形態について説明したが、車両１０の代わりに、建造物、航空機、又は船舶に実装してもよい。図２５は、建造物１０Ａ、航空機１０Ｂ、船舶１０Ｃを示す図である。
　図２５（Ａ）に示す建造物１０Ａの屋上には日射計１０１が配置されている。主制御部１１０、算出部１２０、補正部１３０、推定部１４０、計算部１５０、メモリ１６０を有する電子制御装置は、一例として、建造物１０Ａの内部に配置すればよい。一例として、基準面を建造物１０Ａの屋上の表面に設定し、建造物１０Ａの窓ガラス１２Ａの日射熱取得量ＱＴＳを求め、複数の窓ガラス１２Ａから内部に入り込む熱量を表す総日射熱取得量ＣＱＴＳを求めて、建造物１０Ａの空調装置の制御に用いればよい。また、空調装置の制御に加えて、電動ブラインド等の開閉制御を行ってもよい。なお、建造物１０Ａは、図２５（Ａ）に示すようなビルディングに限らず、マンションや戸建て住宅、体育館や倉庫等であってもよい。なお、総日射熱取得量ＣＱＴＳに基づいて空調装置の制御を行うことは、複数の窓ガラス１２Ａの各々についての全天日射量ＪＴＳに基づいて空調装置の制御を行うことである。
　ここでは、建造物１０Ａの屋上に日射計１０１が配置される形態について説明したが、日射計１０１は、建造物１０Ａの内部に配置してもよい。１つの日射計１０１を用いて、建造物１０Ａの室内に入ってくる総日射熱取得量を高精度に算出できるため、空調装置の制御を非常に精緻に行うことができる。
　図２５（Ｂ）に示す航空機１０Ｂの屋根の上には日射計１０１が配置されている。主制御部１１０、算出部１２０、補正部１３０、推定部１４０、計算部１５０、メモリ１６０を有する電子制御装置は、一例として、航空機１０Ｂの内部に配置すればよい。一例として、基準面を航空機１０Ｂの屋根の表面に設定し、航空機１０Ｂの窓ガラス１２Ｂの日射熱取得量ＱＴＳを求め、複数の窓ガラス１２Ｂから内部に入り込む熱量を表す総日射熱取得量ＣＱＴＳを求めて、航空機１０Ｂの空調装置の制御に用いればよい。なお、総日射熱取得量ＣＱＴＳに基づいて空調装置の制御を行うことは、複数の窓ガラス１２Ｂの各々についての全天日射量ＪＴＳに基づいて空調装置の制御を行うことである。
　ここでは、航空機１０Ｂの屋根の上に日射計１０１が配置される形態について説明したが、日射計１０１は、航空機１０Ｂの内部に配置してもよい。１つの日射計１０１を用いて、航空機１０Ｂの室内に入ってくる総日射熱取得量を高精度に算出できるため、空調装置の制御を非常に精緻に行うことができる。
　図２５（Ｃ）に示す船舶１０Ｃの屋根の上には日射計１０１が配置されている。主制御部１１０、算出部１２０、補正部１３０、推定部１４０、計算部１５０、メモリ１６０を有する電子制御装置は、一例として、船舶１０Ｃの内部に配置すればよい。一例として、基準面を船舶１０Ｃの屋根の表面に設定し、船舶１０Ｃの窓ガラス１２Ｃの日射熱取得量ＱＴＳを求め、複数の窓ガラス１２Ｃから内部に入り込む熱量を表す総日射熱取得量ＣＱＴＳを求めて、船舶１０Ｃの空調装置の制御に用いればよい。なお、総日射熱取得量ＣＱＴＳに基づいて空調装置の制御を行うことは、複数の窓ガラス１２Ｃの各々についての全天日射量ＪＴＳに基づいて空調装置の制御を行うことである。
　ここでは、船舶１０Ｃの屋根の上に日射計１０１が配置される形態について説明したが、日射計１０１は、船舶１０Ｃの内部に配置してもよい。１つの日射計１０１を用いて、船舶１０Ｃの室内に入ってくる総日射熱取得量を高精度に算出できるため、空調装置の制御を非常に精緻に行うことができる。
　以上、本発明の例示的な実施の形態の日射推定システム、空調制御システム、空調装置、車両、建造物、日射推定方法、空調制御方法、及び日射推定プログラムについて説明したが、本発明は、具体的に開示された実施の形態に限定されるものではなく、特許請求の範囲から逸脱することなく、種々の変形や変更が可能である。
　なお、本国際出願は、２０２０年３月２７日に出願した日本国特許出願２０２０－０５７９５３号に基づく優先権を主張するものであり、その全内容は本国際出願にここでの参照により援用されるものとする。
　なお、以下の項目を開示する。
（項目１）
　対象物の外面にある基準面で受ける日射の第１全天日射量を取得する取得部と、
　太陽定数と、大気透過率と、天空日射透過率と、前記基準面の向きに関係する第１関係値とに基づいて、前記第１全天日射量の取得時の前記基準面における第２全天日射量を算出する算出部と、
　前記第２全天日射量が前記第１全天日射量に等しくなるように、前記大気透過率及び前記天空日射透過率をそれぞれ補正して補正大気透過率及び補正天空日射透過率を求める補正部と、
　前記太陽定数と、前記補正大気透過率と、前記補正天空日射透過率と、前記対象物の外面のうちの所定外面の向きに関係する第２関係値とに基づいて、前記第１全天日射量の取得時の前記所定外面における全天日射量の推定値である推定全天日射量を算出する推定部と
　を含む、日射推定システム。
（項目２）
　前記算出部は、
　　　前記太陽定数と、前記大気透過率と、前記第１関係値とに基づいて第１直達日射量
　　を算出し、
　　　前記太陽定数と、前記大気透過率と、前記天空日射透過率と、前記第１関係値とに
　　基づいて第１天空日射量を算出し、
　　　前記第１直達日射量と前記第１天空日射量とを加算することで、前記第２全天日射
　　量を算出し、
　前記推定部は、
　　　前記太陽定数と、前記補正大気透過率と、前記第２関係値とに基づいて第２直達日
　　射量を算出し、
　　　前記太陽定数と、前記補正大気透過率と、前記補正天空日射透過率と、前記第２関
　　係値とに基づいて第２天空日射量を算出し、
　　　前記第２直達日射量と前記第２天空日射量とを加算することで、前記推定全天日射
　　量を算出する、
　項目１記載の日射推定システム。
（項目３）
　前記第１直達日射量の算出に用いられる前記第１関係値は、前記基準面における直達日射の入射角の余弦であり、
　前記第２直達日射量の算出に用いられる前記第２関係値は、前記所定外面における直達日射の入射角の余弦である、項目２記載の日射推定システム。
（項目４）
　前記第１天空日射量の算出に用いられる前記第１関係値は、前記基準面の傾斜角であり、
　前記第２天空日射量の算出に用いられる前記第２関係値は、前記所定外面の傾斜角である、項目２又は３記載の日射推定システム。
（項目５）
　前記算出部は、さらに、前記第１直達日射量と、前記第１天空日射量と、前記第１関係値と、地面反射率とに基づいて、前記基準面における第１地面反射日射量を算出し、前記第１直達日射量と、前記第１天空日射量と、前記第１地面反射日射量とを加算することで、前記第２全天日射量を算出し、
　前記推定部は、さらに、前記第１直達日射量と、前記第２天空日射量と、前記第２関係値と、前記地面反射率とに基づいて、前記所定外面における第２地面反射日射量を算出し、前記第２直達日射量と、前記第２天空日射量と、前記第２地面反射日射量とを加算することで、前記推定全天日射量を算出する、項目２乃至４のいずれか一項記載の日射推定システム。
（項目６）
　前記第１地面反射日射量の算出に用いられる前記第１関係値は、前記基準面の傾斜角であり、
　前記第２地面反射日射量の算出に用いられる前記第２関係値は、前記所定外面の傾斜角である、項目５記載の日射推定システム。
（項目７）
　前記推定全天日射量に前記所定外面の日射熱取得率を乗算して前記所定外面における日射熱取得量を計算する計算部をさらに含む、項目１乃至６のいずれか一項記載の日射推定システム。
（項目８）
　前記計算部は、複数の前記所定外面の各々について前記日射熱取得量を計算し、複数の前記日射熱取得量の総和である総日射熱取得量を計算する、項目７記載の日射推定システム。
（項目９）
　前記取得部は、
　（Ａ）前記第１全天日射量を前記基準面で測定する日射計、
　（Ｂ）前記対象物の透過壁体を透過した日射の全天日射量を測定する日射計と、前記日射計によって測定された全天日射量を前記透過壁体の光学特性に基づいて前記第１全天日射量に変換する変換部とを有する取得部、及び
　（Ｃ）前記第１全天日射量を表すデータを受信する受信部
　のうちのいずれか１つである、項目１乃至８のいずれか一項記載の日射推定システム。
（項目１０）
　前記対象物は、車両、船舶、航空機、及び建造物のうちのいずれか１つである、項目１乃至９のいずれか一項記載の日射推定システム。
（項目１１）
　項目１乃至９のいずれか一項記載の日射推定システムと、
　空調装置と、
　前記日射推定システムによって推定される前記推定全天日射量に基づいて、前記空調装置の空調制御を行う空調制御部と
　を含む、空調制御システム。
（項目１２）
　前記対象物は車両であり、
　前記日射推定システムは、前記車両に搭載される、電子制御装置、位置データ取得部、ナビゲーション制御部、及び通信部のうちのいずれか１つ以上からデータを取得し、取得したデータを用いて前記推定全天日射量を算出する、項目１１記載の空調制御システム。
（項目１３）
　項目１乃至９のいずれか一項記載の日射推定システムの前記算出部、前記補正部、及び前記推定部を有する制御装置を含む、空調装置。
（項目１４）
　項目１乃至９のいずれか一項記載の日射推定システムを含む、車両。
（項目１５）
　項目１乃至９のいずれか一項記載の日射推定システムを含む、船舶。
（項目１６）
　項目１乃至９のいずれか一項記載の日射推定システムを含む、航空機。
（項目１７）
　項目１乃至９のいずれか一項記載の日射推定システムを含む、建造物。
（項目１８）
　対象物の外面にある基準面で受ける日射の第１全天日射量を取得し、
　太陽定数と、大気透過率と、天空日射透過率と、前記基準面の向きに関係する第１関係値とに基づいて、前記第１全天日射量の取得時の前記基準面における第２全天日射量を算出し、
　前記第２全天日射量が前記第１全天日射量に等しくなるように、前記大気透過率及び前記天空日射透過率をそれぞれ補正して補正大気透過率及び補正天空日射透過率を求め、
　前記太陽定数と、前記補正大気透過率と、前記補正天空日射透過率と、前記対象物の外面のうちの所定外面の向きに関係する第２関係値とに基づいて、前記第１全天日射量の取得時の前記所定外面における全天日射量の推定値である推定全天日射量を算出する
　日射推定方法。
（項目１９）
　前記第２全天日射量を算出することは、
　　　前記太陽定数と、前記大気透過率と、前記第１関係値とに基づいて第１直達日射量
　　を算出し、
　　　前記太陽定数と、前記大気透過率と、前記天空日射透過率と、前記第１関係値とに
　　基づいて第１天空日射量を算出し、
　　　前記第１直達日射量と前記第１天空日射量とを加算することで、前記第２全天日射
　　量を算出することであり、
　前記推定全天日射量を算出することは、
　　　前記太陽定数と、前記補正大気透過率と、前記第２関係値とに基づいて第２直達日
　　射量を算出し、
　　　前記太陽定数と、前記補正大気透過率と、前記補正天空日射透過率と、前記第２関
　　係値とに基づいて第２天空日射量を算出し、
　　　前記第２直達日射量と前記第２天空日射量とを加算することで、前記推定全天日射
　　量を算出することである、
　項目１８記載の日射推定方法。
（項目２０）
　前記第１直達日射量の算出に用いられる前記第１関係値は、前記基準面における直達日射の入射角の余弦であり、
　前記第２直達日射量の算出に用いられる前記第２関係値は、前記所定外面における直達日射の入射角の余弦である、項目１９記載の日射推定方法。
（項目２１）
　前記第１天空日射量の算出に用いられる前記第１関係値は、前記基準面の傾斜角であり、
　前記第２天空日射量の算出に用いられる前記第２関係値は、前記所定外面の傾斜角である、項目１９又は２０記載の日射推定方法。
（項目２２）
　前記第２全天日射量を算出することは、さらに、前記第１直達日射量と、前記第１天空日射量と、前記第１関係値と、地面反射率とに基づいて、前記基準面における第１地面反射日射量を算出し、前記第１直達日射量と、前記第１天空日射量と、前記第１地面反射日射量とを加算することで、前記第２全天日射量を算出することであり、
　前記推定全天日射量を算出することは、さらに、前記第１直達日射量と、前記第２天空日射量と、前記第２関係値と、地面反射率とに基づいて、前記所定外面における第２地面反射日射量を算出し、前記第２直達日射量と、前記第２天空日射量と、前記第２地面反射日射量とを加算することで、前記推定全天日射量を算出することである、項目１３乃至２１のいずれか一項記載の日射推定方法。
（項目２３）
　前記第１地面反射日射量の算出に用いられる前記第１関係値は、前記基準面の傾斜角であり、
　前記第２地面反射日射量の算出に用いられる前記第２関係値は、前記所定外面の傾斜角である、項目２２記載の日射推定方法。
（項目２４）
　前記推定全天日射量に前記所定外面の日射熱取得率を乗算して前記所定外面における日射熱取得量を計算することをさらに含む、項目１８乃至２３のいずれか一項記載の日射推定方法。
（項目２５）
　前記計算することは、複数の前記所定外面の各々について前記日射熱取得量を計算し、複数の前記日射熱取得量の総和である総日射熱取得量を計算することである、項目２４記載の日射推定方法。
（項目２６）
　前記第１全天日射量を取得することは、
　（Ａ）前記基準面に設けられた日射計で測定される全天日射量を前記第１全天日射量として取得すること、
　（Ｂ）前記対象物の透過壁体を透過した日射の全天日射量を測定することと、前記測定された全天日射量を前記透過壁体の光学特性に基づいて前記第１全天日射量に変換することとによって、前記第１全天日射量を取得すること、及び
　（Ｃ）前記第１全天日射量を表すデータを受信すること
　のうちのいずれか１つである、項目１８乃至２５のいずれか一項記載の日射推定方法。
（項目２７）
　前記対象物は、車両、船舶、航空機、及び建造物のうちのいずれか１つである、項目１８乃至２６のいずれか一項記載の日射推定方法。
（項目２８）
　項目１８乃至２６のいずれか一項記載の日射推定方法によって推定される前記推定全天日射量に基づいて、空調制御を行う、空調制御方法。
（項目２９）
　前記対象物は車両であり、
　前記車両に搭載される、電子制御装置、位置データ取得部、ナビゲーション制御部、及び通信部のうちのいずれか１つ以上からデータを取得し、取得したデータを用いて前記推定全天日射量を算出する、項目２８記載の空調制御方法。
（項目３０）
　太陽定数と、大気透過率と、天空日射透過率と、基準面の向きに関係する第１関係値とに基づいて、対象物の外面にある基準面で受ける日射の第１全天日射量の取得時の前記基準面における第２全天日射量を算出することと、
　前記第２全天日射量が前記第１全天日射量に等しくなるように、前記大気透過率及び前記天空日射透過率をそれぞれ補正して補正大気透過率及び補正天空日射透過率を求めることと、
　前記太陽定数と、前記補正大気透過率と、前記補正天空日射透過率と、前記対象物の外面のうちの所定外面の向きに関係する第２関係値とに基づいて、前記第１全天日射量の取得時の前記所定外面における全天日射量の推定値である推定全天日射量を算出することと
　を含む処理をコンピュータに実行させる、日射推定プログラム。
　１０　車両（対象物の一例）
　１０Ａ　建造物（対象物の一例）
　１０Ｂ　航空機（対象物の一例）
　１０Ｃ　船舶（対象物の一例）
　１１　ルーフ（基準面の一例）
　１２ＷＳ、１２ＦＲ、１２ＲＲ、１２ＦＬ、１２ＲＬ、１２ＲＧ、１２Ａ、１２Ｂ、１２Ｃ　窓ガラス（透過壁体の一例）
　１００　日射推定システム
　１０１　日射計
　１２０　算出部
　１３０　補正部
　１４０　推定部
　１５０　計算部
　２１０　制御部（空調制御部の一例）
　３００Ａ　パワートレインＥＣＵ（電子制御装置の一例）
　３００Ｂ　ナビゲーションＥＣＵ（ナビゲーション制御部の一例）
　３００Ｃ　ＥＣＵ（電子制御装置の一例）
　３００Ｄ　通信部（通信部の一例、受信部の一例）
　４７０　ＨＶＡＣユニット (空調装置の一例)

Claims (22)

  1. 　対象物の外面にある基準面で受ける日射の第１全天日射量を取得する取得部と、
    　太陽定数と、大気透過率と、天空日射透過率と、前記基準面の向きに関係する第１関係値とに基づいて、前記第１全天日射量の取得時の前記基準面における第２全天日射量を算出する算出部と、
    　前記第２全天日射量が前記第１全天日射量に等しくなるように、前記大気透過率及び前記天空日射透過率をそれぞれ補正して補正大気透過率及び補正天空日射透過率を求める補正部と、
    　前記太陽定数と、前記補正大気透過率と、前記補正天空日射透過率と、前記対象物の外面のうちの所定外面の向きに関係する第２関係値とに基づいて、前記第１全天日射量の取得時の前記所定外面における全天日射量の推定値である推定全天日射量を算出する推定部と
    　を含む、日射推定システム。
  2. 　前記算出部は、
    　　　前記太陽定数と、前記大気透過率と、前記第１関係値とに基づいて第１直達日射量
    　　を算出し、
    　　　前記太陽定数と、前記大気透過率と、前記天空日射透過率と、前記第１関係値とに
    　　基づいて第１天空日射量を算出し、
    　　　前記第１直達日射量と前記第１天空日射量とを加算することで、前記第２全天日射
    　　量を算出し、
    　前記推定部は、
    　　　前記太陽定数と、前記補正大気透過率と、前記第２関係値とに基づいて第２直達日
    　　射量を算出し、
    　　　前記太陽定数と、前記補正大気透過率と、前記補正天空日射透過率と、前記第２関
    　　係値とに基づいて第２天空日射量を算出し、
    　　　前記第２直達日射量と前記第２天空日射量とを加算することで、前記推定全天日射
    　　量を算出する、
    　請求項１記載の日射推定システム。
  3. 　前記算出部は、さらに、前記第１直達日射量と、前記第１天空日射量と、前記第１関係値と、地面反射率とに基づいて、前記基準面における第１地面反射日射量を算出し、前記第１直達日射量と、前記第１天空日射量と、前記第１地面反射日射量とを加算することで、前記第２全天日射量を算出し、
    　前記推定部は、さらに、前記第１直達日射量と、前記第２天空日射量と、前記第２関係値と、前記地面反射率とに基づいて、前記所定外面における第２地面反射日射量を算出し、前記第２直達日射量と、前記第２天空日射量と、前記第２地面反射日射量とを加算することで、前記推定全天日射量を算出する、請求項２記載の日射推定システム。
  4. 　前記推定全天日射量に前記所定外面の日射熱取得率を乗算して前記所定外面における日射熱取得量を計算する計算部をさらに含む、請求項１乃至３のいずれか一項記載の日射推定システム。
  5. 　前記計算部は、複数の前記所定外面の各々について前記日射熱取得量を計算し、複数の前記日射熱取得量の総和である総日射熱取得量を計算する、請求項４記載の日射推定システム。
  6. 　前記取得部は、
    　（Ａ）前記第１全天日射量を前記基準面で測定する日射計、
    　（Ｂ）前記対象物の透過壁体を透過した日射の全天日射量を測定する日射計と、前記日射計によって測定された全天日射量を前記透過壁体の光学特性に基づいて前記第１全天日射量に変換する変換部とを有する取得部、及び
    　（Ｃ）前記第１全天日射量を表すデータを受信する受信部
    　のうちのいずれか１つである、請求項１乃至５のいずれか一項記載の日射推定システム。
  7. 　前記対象物は、車両、船舶、航空機、及び建造物のうちのいずれか１つである、請求項１乃至６のいずれか一項記載の日射推定システム。
  8. 　請求項１乃至６のいずれか一項記載の日射推定システムと、
    　空調装置と、
    　前記日射推定システムによって推定される前記推定全天日射量に基づいて、前記空調装置の空調制御を行う空調制御部と
    　を含む、空調制御システム。
  9. 　前記対象物は車両であり、
    　前記日射推定システムは、前記車両に搭載される、電子制御装置、位置データ取得部、ナビゲーション制御部、及び通信部のうちのいずれか１つ以上からデータを取得し、取得したデータを用いて前記推定全天日射量を算出する、請求項８記載の空調制御システム。
  10. 　請求項１乃至６のいずれか一項記載の日射推定システムの前記算出部、前記補正部、及び前記推定部を有する制御装置を含む、空調装置。
  11. 　請求項１乃至６のいずれか一項記載の日射推定システムを含む、車両。
  12. 　請求項１乃至６のいずれか一項記載の日射推定システムを含む、建造物。
  13. 　対象物の外面にある基準面で受ける日射の第１全天日射量を取得し、
    　太陽定数と、大気透過率と、天空日射透過率と、前記基準面の向きに関係する第１関係値とに基づいて、前記第１全天日射量の取得時の前記基準面における第２全天日射量を算出し、
    　前記第２全天日射量が前記第１全天日射量に等しくなるように、前記大気透過率及び前記天空日射透過率をそれぞれ補正して補正大気透過率及び補正天空日射透過率を求め、
    　前記太陽定数と、前記補正大気透過率と、前記補正天空日射透過率と、前記対象物の外面のうちの所定外面の向きに関係する第２関係値とに基づいて、前記第１全天日射量の取得時の前記所定外面における全天日射量の推定値である推定全天日射量を算出する
    　日射推定方法。
  14. 　前記第２全天日射量を算出することは、
    　　　前記太陽定数と、前記大気透過率と、前記第１関係値とに基づいて第１直達日射量
    　　を算出し、
    　　　前記太陽定数と、前記大気透過率と、前記天空日射透過率と、前記第１関係値とに
    　　基づいて第１天空日射量を算出し、
    　　　前記第１直達日射量と前記第１天空日射量とを加算することで、前記第２全天日射
    　　量を算出することであり、
    　前記推定全天日射量を算出することは、
    　　　前記太陽定数と、前記補正大気透過率と、前記第２関係値とに基づいて第２直達日
    　　射量を算出し、
    　　　前記太陽定数と、前記補正大気透過率と、前記補正天空日射透過率と、前記第２関
    　　係値とに基づいて第２天空日射量を算出し、
    　　　前記第２直達日射量と前記第２天空日射量とを加算することで、前記推定全天日射
    　　量を算出することである、
    　請求項１３記載の日射推定方法。
  15. 　前記第２全天日射量を算出することは、さらに、前記第１直達日射量と、前記第１天空日射量と、前記第１関係値と、地面反射率とに基づいて、前記基準面における第１地面反射日射量を算出し、前記第１直達日射量と、前記第１天空日射量と、前記第１地面反射日射量とを加算することで、前記第２全天日射量を算出することであり、
    　前記推定全天日射量を算出することは、さらに、前記第１直達日射量と、前記第２天空日射量と、前記第２関係値と、地面反射率とに基づいて、前記所定外面における第２地面反射日射量を算出し、前記第２直達日射量と、前記第２天空日射量と、前記第２地面反射日射量とを加算することで、前記推定全天日射量を算出することである、請求項１４記載の日射推定方法。
  16. 　前記推定全天日射量に前記所定外面の日射熱取得率を乗算して前記所定外面における日射熱取得量を計算することをさらに含む、請求項１４又は１５記載の日射推定方法。
  17. 　前記計算することは、複数の前記所定外面の各々について前記日射熱取得量を計算し、複数の前記日射熱取得量の総和である総日射熱取得量を計算することである、請求項１６記載の日射推定方法。
  18. 　前記第１全天日射量を取得することは、
    　（Ａ）前記基準面に設けられた日射計で測定される全天日射量を前記第１全天日射量として取得すること、
    　（Ｂ）前記対象物の透過壁体を透過した日射の全天日射量を測定することと、前記測定された全天日射量を前記透過壁体の光学特性に基づいて前記第１全天日射量に変換することとによって、前記第１全天日射量を取得すること、及び
    　（Ｃ）前記第１全天日射量を表すデータを受信すること
    　のうちのいずれか１つである、請求項１４乃至１７のいずれか一項記載の日射推定方法。
  19. 　前記対象物は、車両、船舶、航空機、及び建造物のうちのいずれか１つである、請求項１４乃至１８のいずれか一項記載の日射推定方法。
  20. 　請求項１４乃至１８のいずれか一項記載の日射推定方法によって推定される前記推定全天日射量に基づいて、空調制御を行う、空調制御方法。
  21. 　前記対象物は車両であり、
    　前記車両に搭載される、電子制御装置、位置データ取得部、ナビゲーション制御部、及び通信部のうちのいずれか１つ以上からデータを取得し、取得したデータを用いて前記推定全天日射量を算出する、請求項２０記載の空調制御方法。
  22. 　太陽定数と、大気透過率と、天空日射透過率と、基準面の向きに関係する第１関係値とに基づいて、対象物の外面にある基準面で受ける日射の第１全天日射量の取得時の前記基準面における第２全天日射量を算出することと、
    　前記第２全天日射量が前記第１全天日射量に等しくなるように、前記大気透過率及び前記天空日射透過率をそれぞれ補正して補正大気透過率及び補正天空日射透過率を求めることと、
    　前記太陽定数と、前記補正大気透過率と、前記補正天空日射透過率と、前記対象物の外面のうちの所定外面の向きに関係する第２関係値とに基づいて、前記第１全天日射量の取得時の前記所定外面における全天日射量の推定値である推定全天日射量を算出することと
    　を含む処理をコンピュータに実行させる、日射推定プログラム。
PCT/JP2021/010883 2020-03-27 2021-03-17 日射推定システム、空調制御システム、空調装置、車両、建造物、日射推定方法、空調制御方法、及び日射推定プログラム WO2021193296A1 (ja)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2022510009A JPWO2021193296A1 (ja) 2020-03-27 2021-03-17

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP2020057953 2020-03-27
JP2020-057953 2020-03-27

Publications (1)

Publication Number Publication Date
WO2021193296A1 true WO2021193296A1 (ja) 2021-09-30

Family

ID=77890254

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2021/010883 WO2021193296A1 (ja) 2020-03-27 2021-03-17 日射推定システム、空調制御システム、空調装置、車両、建造物、日射推定方法、空調制御方法、及び日射推定プログラム

Country Status (2)

Country Link
JP (1) JPWO2021193296A1 (ja)
WO (1) WO2021193296A1 (ja)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20220396116A1 (en) * 2021-06-11 2022-12-15 Hyundai Motor Company Air conditioning control apparatus and method

Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH10108486A (ja) * 1996-09-26 1998-04-24 Kansai Electric Power Co Inc:The 太陽光発電システムの発電量予測方法
JP2002362129A (ja) * 2001-02-02 2002-12-18 Denso Corp 車両の日射検出装置及びそれを用いた車両用空調装置
JP2003021688A (ja) * 2001-07-06 2003-01-24 Honda Motor Co Ltd 日射センサ
US20070131782A1 (en) * 2005-12-08 2007-06-14 Ziehr Lawrence P Vehicle position data enhanced solar sensing for vehicle HVAC applications
JP2011163973A (ja) * 2010-02-10 2011-08-25 Osaka Gas Co Ltd 日射量推定装置、太陽光発電量推定装置及びシステム
JP2011247785A (ja) * 2010-05-28 2011-12-08 Denso Corp 車両用日射検出装置
WO2012111215A1 (ja) * 2011-02-18 2012-08-23 Ｊｘ日鉱日石エネルギー株式会社 日射評価装置、日射評価方法及び日射評価プログラム
JP2015059923A (ja) * 2013-09-20 2015-03-30 株式会社東芝 日射量見積装置、日射量見積方法および日射量見積プログラム
US20150096733A1 (en) * 2013-10-07 2015-04-09 GM Global Technology Operations LLC System and method for compensating for solar load
KR20180032943A (ko) * 2016-09-23 2018-04-02 현대자동차주식회사 시간 정보를 활용한 차량 입사 일사량 검출 및 보상제어 방법 및 차량

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH10108486A (ja) * 1996-09-26 1998-04-24 Kansai Electric Power Co Inc:The 太陽光発電システムの発電量予測方法
JP2002362129A (ja) * 2001-02-02 2002-12-18 Denso Corp 車両の日射検出装置及びそれを用いた車両用空調装置
JP2003021688A (ja) * 2001-07-06 2003-01-24 Honda Motor Co Ltd 日射センサ
US20070131782A1 (en) * 2005-12-08 2007-06-14 Ziehr Lawrence P Vehicle position data enhanced solar sensing for vehicle HVAC applications
JP2011163973A (ja) * 2010-02-10 2011-08-25 Osaka Gas Co Ltd 日射量推定装置、太陽光発電量推定装置及びシステム
JP2011247785A (ja) * 2010-05-28 2011-12-08 Denso Corp 車両用日射検出装置
WO2012111215A1 (ja) * 2011-02-18 2012-08-23 Ｊｘ日鉱日石エネルギー株式会社 日射評価装置、日射評価方法及び日射評価プログラム
JP2015059923A (ja) * 2013-09-20 2015-03-30 株式会社東芝 日射量見積装置、日射量見積方法および日射量見積プログラム
US20150096733A1 (en) * 2013-10-07 2015-04-09 GM Global Technology Operations LLC System and method for compensating for solar load
KR20180032943A (ko) * 2016-09-23 2018-04-02 현대자동차주식회사 시간 정보를 활용한 차량 입사 일사량 검출 및 보상제어 방법 및 차량

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20220396116A1 (en) * 2021-06-11 2022-12-15 Hyundai Motor Company Air conditioning control apparatus and method

Also Published As

Publication number Publication date
JPWO2021193296A1 (ja) 2021-09-30

Similar Documents

Publication Publication Date Title
CN106864205B (zh) 一种车辆用太阳能空调的温度补偿方法及系统
CN102338869B (zh) 下行短波辐射和光合有效辐射数据的反演方法及系统
CN101450716B (zh) 一种地球同步转移轨道在轨卫星故障光电探测方法
WO2021193296A1 (ja) 日射推定システム、空調制御システム、空調装置、車両、建造物、日射推定方法、空調制御方法、及び日射推定プログラム
US20130066585A1 (en) Predicted mean vote estimating device and computer program product
CN106932101A (zh) Hj‑1b卫星中红外和热红外通道联合的地表温度反演方法
CN103076616A (zh) 大气对流层电离层电波折射误差一体化实时修正装置
Vasarevičius et al. Solar irradiance model for solar electric panels and solar thermal collectors in Lithuania
US20140207431A1 (en) Urban mapping technique for photovoltaic potential of rooftops
CN116045970B (zh) 一种基于外部条件约束的多平台信息协同导航增强方法
JP2003021688A (ja) 日射センサ
CN106646430A (zh) 一种基于地面探测器的激光足印中心确定方法
CN114266984B (zh) 一种利用高分辨遥感影像计算建筑物屋顶光伏可改造区减碳量的方法
CN110103829B (zh) 车载显示屏的显示方法、装置、车载显示屏和车辆
Suttles et al. Angular radiation models for earth-atmosphere system. Volume 2: Longwave radiation
Korablev et al. AOST: Fourier spectrometer for studying Mars and Phobos
Weiss et al. Albedo of the ice covered Weddell and Bellingshausen Seas
CN105082936B (zh) 用于控制气候控制系统的系统和方法
CN103529861B (zh) 智能窗窗页基于太阳光线的转角方法
CN113787878B (zh) 基于气象数据演算整车热负荷的氢能汽车空调控制方法
Panão et al. Determining the shading correction factor using a smartphone camera with a fisheye lens
CN106324275B (zh) 一种基于双视角光学遥感图像的海面风速探测方法
Escobar et al. Solar Energy Resource Assessment in Chile: Satellite estimation and ground station measurement
CN109375287B (zh) 与空调负荷相关的综合气象因素实感温度修正方法
Inoue et al. Aircraft observations of air-mass modification over the sea of Okhotsk during sea-ice growth

Legal Events

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

Ref document number: 21774307

Country of ref document: EP

Kind code of ref document: A1

ENP Entry into the national phase

Ref document number: 2022510009

Country of ref document: JP

Kind code of ref document: A

NENP Non-entry into the national phase

Ref country code: DE

122 Ep: pct application non-entry in european phase

Ref document number: 21774307

Country of ref document: EP

Kind code of ref document: A1