(57)【要約】【課題】熱流束センサを用いる必要なく、座席と着座者の間の熱の流れを推定する技術を提供する。【解決手段】熱流束推定装置であるエアコンＥＣＵ１４は、赤外線カメラ４の出力に基づいて、人が車両１の座席２、３に着座する際の当該座席の温度である初期温度Ｔｓｉを推定する。そして、エアコンＥＣＵ１４は、人が当該座席に着座した着座時点ｔ＝１の後の期間における、当該座席と人との間の熱流束の推移を算出する。そしてエアコンＥＣＵ１４は、熱流束の算出において、着座時点を起点とする経過時間ｔに依存して変化する時間依存要素ｌｎ（ｔ）と初期温度Ｔｓｉに依存する初期温度依存要素ωとを含む第１計算式Σ１により、着座時点より後の期間における熱流束の推移を算出する。【選択図】図１
2020063029.tif 000002
【特許請求の範囲】
【請求項１】
センサ（４、１６、１７、１８）の出力に基づいて、人が車両（１）の座席（２、３）に着座する際の前記座席の温度である初期温度（Ｔｓｉ）を推定する温度推定部（Ｓ１３０、Ｓ１３５）と、
前記人が前記座席に着座した着座時点の後の期間における、前記座席と前記人との間の熱流束の推移を算出する熱流算出部（Ｓ１４０）と、を備え、
前記熱流算出部は、前記着座時点を起点とする経過時間に依存して変化する時間依存要素（ｌｎ（ｔ））と前記初期温度に依存する初期温度依存要素（ω）とを含む第１計算式（Σ１）により、前記着座時点より後の期間における前記熱流束の推移を算出する、熱流束推定装置。
【請求項２】
前記初期温度依存要素の前記第１計算式に対する寄与度は、前記着座時点以外の時点における前記座席の温度の推定値を含む要素の、前記第１計算式に対する寄与度よりも、高い、請求項１に記載の熱流束推定装置。
【請求項３】
前記熱流算出部は、前記着座時点における前記座席と前記人との間の熱流束である着座時熱流束（Ｒｉｎｉｔ）を前記初期温度に基づいて推定し、
前記時間依存要素は、前記経過時間が長くなるほど前記熱流束が前記着座時熱流束からゼロ側に大きく離れるよう、前記第１計算式に寄与する請求項１に記載の熱流束推定装置。
【請求項４】
前記時間依存要素は、前記着座時点を起点とする経過時間が長くなるほど前記熱流束の変化速度の絶対値が小さくなるよう、前記第１計算式に寄与する請求項３に記載の熱流束推定装置。
【請求項５】
前記初期温度依存要素は、前記初期温度と基準温度（ａ１）との差の絶対値が大きいほど、前記第１計算式に対する前記時間依存要素の寄与を増幅させる請求項１ないし４のいずれか１つに記載の熱流束推定装置。
【請求項６】
前記車両の前記座席には、温度調整装置（５、６、７、８）が配置されており、
前記温度調整装置は、前記座席に着座する人を暖める、前記座席に着座する人を冷やす、前記座席に着座する人に送風する、のうち少なくとも１つを行い、
前記第１計算式は、さらに前記温度調整装置の作動状態に依存する温調依存要素（ε）を含む請求項１ないし５のいずれか１つに記載の熱流束推定装置。
【請求項７】
前記熱流算出部は、前記第１計算式により前記熱流束の推移を算出すると共に、前記第１計算式とは異なる第２計算式（Σ２）により前記熱流束の推移を算出し、
前記第２計算式は、前記座席の温度を含まず、時間経過と共に前記熱流束の経時変化量の絶対値を減少させ、
前記熱流算出部は、前記着座時点の後の切替タイミング（ｔｃ）よりも前においては、前記第１計算式で算出された前記熱流束の推移を採用し、前記切替タイミングの後においては、前記第２計算式で算出された前記熱流束の推移を採用し、
前記切替タイミングは、前記第１計算式で算出された前記熱流束と前記第２計算式で算出された前記熱流束との差の絶対値が所定の許容値より小さくなるタイミングである、請求項１ないし６のいずれか１つに記載の熱流束推定装置。
【請求項８】
前記センサは、前記座席を前記座席の前側から撮影する赤外線カメラである、請求項１ないし７のいずれか１つに記載の熱流束推定装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、熱流束推定装置に関するものである。
【背景技術】
【０００２】
特許文献１には、車両の座席の比較的表面側に埋め込まれた熱流束センサによって、座席と着座者の間の熱の流れを測定する技術が知られている。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】特許第３２２０４０２号公報
【発明の概要】
【発明が解決しようとする課題】
【０００４】
しかしながら、発明者の検討によれば、座席と着座者の間の熱の流れを推定するために、その用途に特化した機能を有する熱流束センサを利用すると、そのセンサを他の用途に
転用することが困難になる。したがって、他の用途への転用の余地が熱流束センサよりも広いセンサを利用して、座席と着座者の間の熱の流れを推定できることが望ましい。本発
明は上記点に鑑み、熱流束センサを用いる必要なく、座席と着座者の間の熱の流れを推定する技術を提供することを目的とする。
【課題を解決するための手段】
【０００５】
上記目的を達成するための請求項１に記載の発明は、センサ（４、１６、１７、１８）の出力に基づいて、人が車両（１）の座席（２、３）に着座する際の前記座席の温度である初期温度（Ｔｓｉ）を推定する温度推定部（Ｓ１３０、Ｓ１３５）と、前記人が前記座席に着座した着座時点の後の期間における、前記座席と前記人との間の熱流束の推移を算出する熱流算出部（Ｓ１４０）と、を備え、前記熱流算出部は、前記着座時点を起点とする経過時間に依存して変化する時間依存要素（ｌｎ（ｔ））と前記初期温度に依存する初期温度依存要素（ω）とを含む第１計算式（Σ１）により、前記着座時点より後の期間における前記熱流束の推移を算出する、熱流束推定装置である。
【０００６】
このように、熱流束推定装置は、人が車両の座席に着座する際の座席の温度である初期温度と着座時点を起点とする経過時間に基づいて、その後の熱流束の推移を算出する。より具体的には、経過時間に依存して変化する時間依存要素と初期温度に依存する初期温度依存要素とを含む第１計算式によって熱流束の推移を算出する。このような第１計算式を用いることで、実際の熱流束に近い熱流束を推定することができた。したがって、結果的に初期温度が推定できるセンサを用いれば、熱流束センサを用いる必要なく、座席と着座者の間の熱の流れを推定することができる。
【０００７】
なお、各構成要素等に付された括弧付きの参照符号は、その構成要素等と後述する実施形態に記載の具体的な構成要素等との対応関係の一例を示すものである。
【図面の簡単な説明】
【０００８】
【図１】車室内に搭載される装置の構成を示す図である。
【図２】車室内に搭載される装置の電気的接続関係を示す図である。
【図３】エアコンＥＣＵが実行する熱流速推定処理のフローチャートである。
【図４】熱流束推定処理のフローチャートである。
【図５】算出された熱流束の推移の一例を表すグラフである。
【図６】ヒータクーラの作動状態とεとの関係を示す表である。
【図７】算出された熱流束と、熱流束センサで実際に測定された熱流束との比較結果を示すグラフである。
【図８】算出された熱流束と、熱流束センサで実際に測定された熱流束との比較結果を示すグラフである。
【図９】第２実施形態においてエアコンＥＣＵが実行する熱流速推定処理のフローチャートである。
【図１０】変形例における熱流束の推移等を示すグラフである。
【図１１】変形例においてヒータクーラの作動状態とεとの関係を示す表である。
【発明を実施するための形態】
【０００９】
（第１実施形態）
以下、本発明の一実施形態について説明する。図１に示すように、本実施形態に係る車両１の車室内には、運転席２、後席３、および不図示の助手席、不図示の他の後席が搭載されている。
【００１０】
運転席２の背もたれ部には、運転席背部ヒータクーラ５が埋め込まれて配置されている。運転席２の座部（すなわち、シートクッション）には、運転席座部ヒータクーラ６が埋め込まれて配置されている。後席３の背もたれ部（すなわち、シートバック）には、後席背部ヒータクーラ７が埋め込まれて配置されている。また、後席３の座部には、後席座部ヒータクーラ８が埋め込まれて配置されている。運転席２、後席３の座部とは、運転席２、後席３に着座した人の臀部および大腿部を支持する部分である。
【００１１】
運転席背部ヒータクーラ５、運転席座部ヒータクーラ６、後席背部ヒータクーラ７、後席座部ヒータクーラ８の各々は、取り付け先の座席２、３に着座する人を、あるときには暖め、またあるときには冷やす温度調整装置である。
【００１２】
例えば、ヒータクーラ５、６、７、８の各々は、ペルチェ素子から構成されていてもよい。取り付け先の座席２、３に着座する人を暖める場合は、ペルチェ素子の両面のうち、当該人側の面が放熱面となるよう、ペルチェ素子に電流が印加される。取り付け先の座席２、３に着座する人を冷やす場合は、ペルチェ素子の両面のうち、当該人側の面が吸熱面となるよう、ペルチェ素子に電流が印加される。
【００１３】
運転席背部ヒータクーラ５の主な温度調整対象（すなわち加熱対象および冷却対象）は、運転席２の背もたれ部および運転席２に着座する人の背中である。運転席座部ヒータクーラ６の主な温度調整対象は、運転席２の座部および運転席２に着座する人の臀部と大腿部である。後席背部ヒータクーラ７の主な温度調整対象は、後席３の背もたれ部および後席３に着座する人の背中である。後席座部ヒータクーラ８の主な温度調整対象は、後席３の座部および後席３に着座する人の臀部と大腿部である。
【００１４】
また、車両のダッシュボード内には、温度調整された風（すなわち空調風）を車室内に吹き出すための空調ユニット９が配置されている。空調ユニット９は、ケーシング１０、エバポレータ１１、ヒータコア１２、エアミックスドア１３等を含んでいる。
【００１５】
ケーシング１０は、空調ユニット９の外殻となる部材であり、このケーシング１０内に導入された空気が、ケーシング１０内で加熱または冷却され、その後車室内に吹き出される。
【００１６】
エバポレータ１１は、ケーシング１０内に配置され、その内部を冷媒が流通する。空調ユニット９内に導入された空気がエバポレータ１１を通過する際、当該空気と冷媒とが熱交換することで、冷媒が蒸発すると共に空気が冷却される。
【００１７】
ヒータコア１２は、ケーシング１０内において、エバポレータ１１よりも空気流れ下流側に配置され、その内部をエンジン冷却水等の熱媒体が流通する。エバポレータ１１を通過した後の空気がヒータコア１２を通過する際、当該空気と熱媒体とが熱交換することで、熱媒体が冷却されると共に空気が加熱される。
【００１８】
エアミックスドア１３は、ケーシング１０内においてエバポレータ１１よりも空気流れ下流側かつヒータコア１２よりも空気流れ上流側に配置される可動部材である。エアミックスドア１３は、位置が変化することにより、エアミックス比率が変化する。エアミックス比率は、エバポレータ１１を通った後の空気がヒータコア１２を通過する空気の流量と、あるいはヒータコア１２を迂回してケーシング１０内のバイパス通路を通過する空気の流量との比率である。
【００１９】
また、ダッシュボード内には、運転席背部ヒータクーラ５、運転席座部ヒータクーラ６、後席背部ヒータクーラ７、後席座部ヒータクーラ８、空調ユニット９の作動を制御するエアコンＥＣＵ１４が配置されている。
【００２０】
また、車室内には、赤外線カメラ４が配置されている。赤外線カメラ４は、あらかじめ決められた固定の撮影範囲内から放射される赤外線を取得し、取得した赤外線に基づいて、当該撮影範囲内の位置毎の表面温度を画素値として表す画像を生成して出力するセンサである。当該画像はサーモグラフィである。
【００２１】
この撮影範囲内には、運転席２、後席３等、その車室内のその他すべての座席の各々の、背もたれ部の表面および座部の表面が含まれる。赤外線カメラ４は、これらすべての座席を前側から撮影する。ここで、背もたれ部の表面および座部の表面とは、背もたれ部および座部のうち、人が着座した場合における人側の面をいう。つまり、背もたれ部の表面および座部の表面は、背もたれ部の前面および座部の上面である。
【００２２】
人が着座していない座席については、背もたれ部の表面および座部の表面の全体が赤外線カメラ４によって撮影可能であるが、人が着座している座席については、背もたれ部の表面および座部の表面の一部のみが赤外線カメラ４によって撮影可能である。
【００２３】
このような撮影範囲を実現する車室内の赤外線カメラ４の位置としては、車両の左右方向中央部分における、フロントガラスの上端部近傍等がある。この部分にルームミラーが設置されている場合は、赤外線カメラ４はルームミラーに取り付けられてもよい。なお、本実施形態においては、運転席２、後席３のいずれにも、熱流束センサは設置されていない。
【００２４】
また、図２に示すように、車両には、ＡＭドアアクチュエータ１５、外気温センサ１６、内気温センサ１７、日射センサ１８、運転席着座センサ１９、後席着座センサ２０等も搭載されている。
【００２５】
ＡＭドアアクチュエータ１５は、エアミックスドア１３を駆動してエアミックスドア１３の位置を変化させるモータ等のアクチュエータである。外気温センサ１６は、車室外の空気（すなわち外気）の温度に応じた検出信号を出力するセンサである。内気温センサ１７は、車室内の空気（すなわち内気）の温度に応じた検出信号を出力するセンサである。日射センサ１８は、車両への日射量に応じた検出信号を出力するセンサである。
【００２６】
運転席着座センサ１９は、運転席２への乗員の着座の有無に応じた検出信号を出力するセンサである。後席着座センサ２０は、後席３への乗員の着座の有無に応じた検出信号を出力するセンサである。運転席着座センサ１９、後席着座センサ２０は、それぞれ、例えば座部に設置される周知のメンブレン式のスイッチであってもよいし、座部に設置される周知の静電容量型のセンサであってもよいし、座部に設置される周知の圧力センサであってもよい。
【００２７】
エアコンＥＣＵ１４は、ＣＰＵ、ＲＡＭ、ＲＯＭ等を有するＣＰＵが、ＲＯＭに記憶されたプログラムを実行することで、後述する各種処理を実現し、それら処理において、ＲＡＭをワークメモリとして使用する。
【００２８】
エアコンＥＣＵ１４には、赤外線カメラ４から繰り返し定期的に（例えば１秒に１回）サーモグラフィが入力される。また、エアコンＥＣＵ１４には、外気温センサ１６、内気温センサ１７、日射センサ１８、運転席着座センサ１９、後席着座センサ２０からの検出信号が入力される。
【００２９】
また、エアコンＥＣＵ１４は、必要に応じて、運転席背部ヒータクーラ５、運転席座部ヒータクーラ６、後席背部ヒータクーラ７、後席座部ヒータクーラ８の作動、非作動および作動時の作動モードを切り替える制御を行う。
【００３０】
これらヒータクーラ５、６、７、８の作動モードとしては、低出力加熱モード、高出力加熱モード、低出力冷却モード、高出力冷却モードの４つがある。低出力加熱モードと高出力加熱モードは、いずれも着座した人を加熱するモードであり、高出力加熱モードの方が低出力加熱モードよりも放熱量が多く、より高温の熱が放出される。つまり、高出力加熱モードの方が低出力加熱モードよりも加熱能力が高い。低出力冷却モードと高出力冷却モードは、いずれも着座した人を冷却するモードであり、高出力冷却モードの方が低出力冷却モードよりも吸熱量が多く、着座した人をより強く冷やす。つまり、高出力冷却モードの方が低出力冷却モードよりも冷却能力が高い。
【００３１】
また、エアコンＥＣＵ１４は、必要に応じて、ＡＭドアアクチュエータ１５の作動を制御することで、エアミックスドア１３の位置によって決まるエアミックス比率を調節する。
【００３２】
より具体的には、エアコンＥＣＵ１４は、ＲＯＭに記録されているプログラムを実行することで、熱流速推定処理１４Ａと車室内温度調整処理１４Ｂをマルチタスクで同時並行的に実行する。
【００３３】
エアコンＥＣＵ１４は、熱流速推定処理１４Ａを実行することで、運転席２、後席３のうち人が着座している席（以下、着座席という）と当該着座席に着座する人（以下、着座者という）との間の熱流束の推移（すなわち経時変化）を推定する。着座席は運転席２、後席３のうち１つだけの場合もあれば、両方の場合もある。エアコンＥＣＵ１４は、熱流速推定処理１４Ａを実行することで、熱流束推定装置として機能する。
【００３４】
推定対象の熱流束は、各着座席の背もたれ部と人の背中との間の熱流束、および、各着座席の座部と人の臀部との間の熱流束である。つまり、人体の部位毎に、着座席との間の熱流束が算出される。ここで、熱流束は、座席とその着座者の境界面の単位面積を単位時間に横切る熱量であり、単位はＷ／ｍ２であり、座席から人へ熱が流れる場合に正となり、人から座席へ熱が流れる場合に負となる。
【００３５】
エアコンＥＣＵ１４は、熱流速推定処理１４Ａにおいて、推定対象の熱流束の各々を繰り返し定期的（例えば１秒に１回）に算出することで、各熱流束の推移すなわち複数時点の各々における熱流束を推定する。熱流速推定処理１４Ａにおける熱流束の算出方法の詳細については後述する。
【００３６】
また、エアコンＥＣＵ１４は、車室内温度調整処理１４Ｂを実行することで、空調ユニット９から車室内に吹き出される空気の温度を調整すると共に、ヒータクーラ５、６、７、８の作動を制御する。
【００３７】
具体的には、エアコンＥＣＵ１４は、車室内温度調整処理１４Ｂにおいて、車両にいる人が不図示の温度設定スイッチを操作して設定した設定温度Ｔｓｅｔと、外気温センサ１６、内気温センサ１７、日射センサ１８の検出信号に基づいた外気温Ｔａｍ、内気温Ｔｒ、および日射量Ｔｓとに基づいて、以下の式に従って目標吹出温度ＴＡＯを算出する。ＴＡＯ＝Ｋｓｅｔ×Ｔｓｅｔ−Ｋｒ×Ｔｒ−Ｋａｍ×Ｔａｍ−Ｋｓ×Ｔｓ＋Ｃ但し、Ｋｓｅｔ、Ｋｒ、Ｋａｍ、Ｋｓは、制御ゲインであり、Ｃは補正用の定数である。
【００３８】
そして、エアコンＥＣＵ１４は、算出した目標吹出温度ＴＡＯに応じて、当該目標吹出温度ＴＡＯと空調ユニット９から車室内に吹き出される空気の温度が同じになるよう、エアミックス比率を決定する。そしてエアコンＥＣＵ１４は、決定したエアミックス比率を実現するためにＡＭドアアクチュエータ１５に指令信号を出力する。ＡＭドアアクチュエータ１５は、この指令信号に従って、当該エアミックス比率を実現するためにエアミックスドア１３を駆動する。また、エアコンＥＣＵ１４は、算出した目標吹出温度ＴＡＯに応じて、不図示の送風ファンを制御することで、空調ユニット９から車室内に吹き出される空気の風量を変化させてもよい。
【００３９】
またエアコンＥＣＵ１４は、車室内温度調整処理１４Ｂにおいて、赤外線カメラ４から取得した最新のサーモグラフィに基づいて、運転席２、助手席、後席３、および他の後席に着席している乗員の体表面（例えば顔の表面）の温度を特定する。サーモグラフィから乗員の体表面の温度を算出する方法としては、例えば、以下のような方法がある。すなわち、エアコンＥＣＵ１４は、周知の画像認識処理によってサーモグラフィから乗員の顔の位置範囲を特定し、特定した位置範囲内の画素に基づいて、乗員の体表面の温度を特定してもよい。
【００４０】
そしてエアコンＥＣＵ１４は、車室内温度調整処理１４Ｂにおいて、特定した体表面の温度に基づいて、上述の目標吹出温度ＴＡＯを補正することにより、車室内に吹き出される空気を調整してもよい。例えば、体表面の温度が基準値よりも高い場合には、上述の式で算出された目標吹出温度ＴＡＯをより低い値に補正してもよい。また例えば、体表面の温度が基準値よりも低い場合には、上述の式で算出された目標吹出温度ＴＡＯをより高い値に補正してもよい。このように、赤外線カメラ４は、乗員の体表面の温度の検出にも用いられる。
【００４１】
また、エアコンＥＣＵ１４は、車室内温度調整処理１４Ｂにおいて、ヒータクーラ５、６、７、８の各々に設けられた不図示の操作スイッチの操作内容に応じて、ヒータクーラ５、６、７、８の各々の、作動オフ、作動オンを切り替える。それと共にエアコンＥＣＵ１４は、車室内温度調整処理１４Ｂにおいて、作動中のヒータクーラ５、６、７、８の各々について、低出力加熱モード、高出力加熱モード、低出力冷却モード、高出力冷却モードを切り替える。
【００４２】
具体的には、エアコンＥＣＵ１４は、不図示の運転席用オンオフスイッチに対する乗員の手動操作に従って、運転席背部ヒータクーラ５および運転席座部ヒータクーラ６の作動オン、作動オフを切り替える。また、エアコンＥＣＵ１４は、不図示の後席用オンオフスイッチに対する乗員の手動操作に従って、後席背部ヒータクーラ７および後席座部ヒータクーラ８の作動オン、作動オフを切り替える。
【００４３】
そしてエアコンＥＣＵ１４は、ヒータクーラ５−８の各々について、その作動中の作動モードを、低出力加熱モード、高出力加熱モード、低出力冷却モード、高出力冷却モードの間で切り替える。切り替えのために用いられるのが、運転席２とその着座者との間の熱流束および後席３とその着座者との間の熱流束である。
【００４４】
まず、運転席背部ヒータクーラ５の制御について説明する。エアコンＥＣＵ１４は、熱流速推定処理１４Ａによって算出した運転席２の背もたれ部と運転席２の着座者の背中との間の熱流束に基づいて、当該着座者の背中における温熱感を特定する。特定される温熱感の値としては、例えば、「寒い」、「快適」、「暑い」、の３つがある。
【００４５】
例えば、夏に乗員が座席に着座した直後から数分は、座席の温度が非常に高い。その結果、座席から乗員に多量の熱が流れ、乗員と座席の接触面が蒸れ、乗員が暑さを感じる。その後、車内の冷房が進行すれば、座席の温度が低下し、座席から乗員への熱の流れは大きく低下し、乗員が暑さを感じることもなくなる。
【００４６】
また例えば、冬に乗員が座席に着座した直後から数分は、座席の温度が非常に低い。その結果、乗員から座席に多量の熱が流れ、乗員が寒さを感じる。その後、車内の暖房が進行すれば、座席の温度が上昇し、乗員から座席への熱の流れは大きく低下し、乗員が寒さを感じることもなくなる。
【００４７】
このような理由から、エアコンＥＣＵ１４は、熱流束が所定の正の基準値Ｐ１よりも大きい場合に温熱感として「暑い」という値を特定する。またエアコンＥＣＵ１４は、熱流束が所定の負の基準値Ｐ２よりも小さい場合に温熱感として「寒い」という値を特定する。またエアコンＥＣＵ１４は、熱流束が基準値Ｐ１以下かつ基準値Ｐ２以上である場合に「快適」という値を特定する。
【００４８】
さらにエアコンＥＣＵ１４は、熱流速推定処理１４Ａによって算出した運転席２の座部と運転席２の着座者の臀部との間の熱流束に基づいて、当該着座者の臀部における温熱感を特定する。さらにエアコンＥＣＵ１４は、熱流速推定処理１４Ａによって算出した後席３の背もたれ部と後席３の着座者の背中との間の最新の（すなわち、現在の）熱流束に基づいて、当該着座者の背中における温熱感を特定する。さらにエアコンＥＣＵ１４は、熱流速推定処理１４Ａによって算出した後席３の座部と後席３の着座者の臀部との間の最新の熱流束に基づいて、当該着座者の臀部における温熱感を特定する。これら３つの温熱感の特定方法は、使用する熱流束が異なる以外は、その他は運転席２の着座者の背中における温熱感の特定方法と同じである。
【００４９】
続いてエアコンＥＣＵ１４は、特定した４つの温熱感に基づいて、ヒータクーラ５−８を制御する。具体的には、運転席２の着座者の背中における温熱感に基づいて運転席背部ヒータクーラ５を制御する。また、運転席２の着座者の臀部における温熱感に基づいて運転席座部ヒータクーラ６を制御する。また、後席３の着座者の背中における温熱感に基づいて後席背部ヒータクーラ７を制御する。また、後席３の着座者の臀部における温熱感に基づいて後席座部ヒータクーラ８を制御する。
【００５０】
制御の方法としては、例えば温熱感が「暑い」となっている部位については、その部位に対応するヒータクーラを高出力冷却モードに設定する。そして、温熱感が「寒い」となっている部位については、その部位に対応するヒータクーラを高出力加熱モードに設定する。そして、温熱感が「快適」となっている部位については、その部位に対応するヒータクーラの作動モードを低出力加熱モードまたは低出力冷却モードに設定する。低出力加熱モードと低出力冷却モードのどちらを選ぶかは、例えば、内気温センサ１７によって検出された車室内温度が閾値以上の場合に低出力冷却モードを選び、車室内温度が閾値未満の場合に低出力加熱モードを選んでもよい。このように、ヒータクーラ５−８によって、運転席２、後席３の乗員の背中および臀部の着座後の温熱感が早期に改善される。
【００５１】
ここで、熱流速推定処理１４Ａの詳細について説明する。上述の通り、エアコンＥＣＵ１４は、熱流速推定処理１４Ａにおいて、運転席２に人が着座していれば、その背もたれ部とその人の背中との間の熱流束の推移を算出し、それと共に、その座部とその人の臀部との間の熱流束の推移を算出する。またエアコンＥＣＵ１４は、後席３に人が着座していれば、その背もたれ部とその人の背中との間の熱流束の推移を算出し、それと共に、その座部とその人の臀部との間の熱流束の推移を算出する。
【００５２】
図３に、エアコンＥＣＵ１４が熱流速推定処理１４Ａを実現するために実行する処理のフローチャートを示す。このフローチャートに示す処理は、エアコンＥＣＵ１４が起動時に不図示のＲＯＭ中に記録された所定のプログラムを実行することで実現される。なお、エアコンＥＣＵ１４は、算出対象の熱流束毎に１つ、この図３の処理を実行する。
【００５３】
したがって、エアコンＥＣＵ１４は、運転席２の背もたれ部用の図３の処理、運転席２の座部用の図３の処理、後席３の背もたれ部用の図３の処理、および、後席３の座部用の図３の処理を、それぞれ独立に同時並行で実行する。つまり、４つの図３の処理が同時並行的に実行される。
【００５４】
運転席２の背もたれ部用の処理では、運転席２の背もたれ部と運転席２の着座者（すなわち運転者）の背中との間の熱流束が算出される。運転席２の座部用の処理では、運転席２の座部と運転席２の着座者の臀部との間の熱流束が算出される。後席３の背もたれ部用の処理では、後席３の背もたれ部と後席３の着座者の背中との間の熱流束が算出される。後席３の座部用の処理では、後席３の座部と後席３の着座者の臀部との間の熱流束が算出される。
【００５５】
以下、１つの図３の処理について説明する。以下の説明は、上記４つの図３の処理のいずれにもあてはまる。１つの図３の処理において、まずエアコンＥＣＵ１４は、ステップＳ１１０で、赤外線カメラ４から取得した最新のサーモグラフィに基づいて、対象の座席の対象の部分の表面温度を特定する。
【００５６】
なお、対象の座席は、着座者との間で熱流束を測定する対象の座席である。また、対象の部分は、対象の座席の背もたれ部と背中との間の熱流束を算出する場合は当該背もたれ部であり、対象の座席の座部と臀部との間の熱流束を算出する場合は当該座部である。また、対象人体部位は、対象の座席の背もたれ部と背中との間の熱流束を算出する場合は当該背中であり、対象の座席の座部と臀部との間の熱流束を算出する場合は当該臀部である。
【００５７】
車両１のメインスイッチがオフの場合もオンの場合も車両１のバッテリから赤外線カメラ４およびエアコンＥＣＵ１４に電力が供給されてそれらが作動するような、赤外線カメラ４およびエアコンＥＣＵ１４の設置形態がある。そのような設置形態による作動を、常時作動という。また、車両１のメインスイッチがオンの場合にのみ赤外線カメラ４およびエアコンＥＣＵ１４に電力が供給されてそれらが作動するような、赤外線カメラ４およびエアコンＥＣＵ１４の設置形態がある。そのような設置形態による作動を、以下、オン時限定作動という。
【００５８】
なお、メインスイッチは、車両の走行のための動力を発生する動力源（例えば、内燃機関、電気モータ等）の始動を可能な状態にするスイッチである。メインスイッチとしては、例えば、内燃機関の動力で走行する車両のイグニッションスイッチ、電気自動車に走行のための動力を供給する電気モータへの通電を許可する主電源スイッチがある。
【００５９】
エアコンＥＣＵ１４がステップＳ１１０を実行するタイミングでは、対象の座席に人は着座していない。常時作動の場合は、車両のメインスイッチがオフの場合にこのような状況が発生し得る。
【００６０】
オン時限定作動の場合は、運転席２については赤外線カメラ４およびエアコンＥＣＵ１４の起動直後のステップＳ１１０の段階で既に人が着座している可能性が高い。しかし、この場合であっても、後席３については、メインスイッチがオンであっても人が着座していない場合が多い。車両１がカーシェアに用いられる場合は特にその傾向が強い。
【００６１】
ステップＳ１１０において、エアコンＥＣＵ１４は、取得した最新のサーモグラフィ中で、対象の座席の対象の部分であるとしてあらかじめ定められた画素範囲内における複数の画素値（すなわち温度値）を抽出する。そして、抽出した複数の画素値の統計代表値（例えば平均値）を、対象の座席の対象の部分の表面温度として特定する。
【００６２】
続いてエアコンＥＣＵ１４は、ステップＳ１２０で、運転席着座センサ１９および後席着座センサ２０のうち対象の座席の着座センサからの信号に基づいて、対象の座席に人が着座しているか否かを判定する。着座していないと判定した場合は、ステップＳ１１０に戻る。したがって、対象の座席に人が着座するまでは、ステップＳ１１０、Ｓ１２０の処理が繰り返される。
【００６３】
ステップＳ１２０で着座したと判定した場合、エアコンＥＣＵ１４はステップＳ１２５に進み、着座時点の時刻をエアコンＥＣＵ１４のＲＡＭに記録する。着座時点は、対象の座席に着座者が着座していないと判定される状態から着座していると判定される状態へ切り替わる時点である。本実施形態では、直前のステップＳ１２０で着座したと判定された時刻が、着座時点の時刻として記録される。
【００６４】
続いてエアコンＥＣＵ１４は、ステップＳ１３０に進み、対象の座席に着座者が着座する直前の、対象の座席の対象の部分の表面温度を初期温度Ｔｓｉとして特定する。具体的には、最後にステップＳ１１０で特定した対象の座席の対象の部分の表面温度を、初期温度Ｔｓｉとして特定する。このような初期温度Ｔｓｉは、対象の座席に人が着座する際の当該座席の温度の推定値に相当する。
【００６５】
続いてエアコンＥＣＵ１４は、ステップＳ１４０に進み、ステップＳ１３０で特定した表面温度に基づいて、ステップＳ１４０の実行時点における対象の座席の対象の部分と対象人体部位との間の熱流束を推定する。熱流束推定の方法については、後で詳述する。ここで推定された熱流束が、車室内温度調整処理１４Ｂにおいて既に説明した通りに用いられる。
【００６６】
ステップＳ１４０に続いては、ステップＳ１５０で、運転席着座センサ１９および後席着座センサ２０のうち対象の座席の着座センサからの信号に基づいて、対象の座席から着座者が離席したか否かを判定する。着座していないと判定した場合、ステップＳ１４０に戻る。したがって、対象の座席から着座者が離れるまでは、ステップＳ１４０、Ｓ１５０の処理が繰り返され、対象の座席の対象の部分と対象人体部位との間の熱流束が逐次算出される。ステップＳ１５０で離席したと判定した場合、ステップＳ１１０に戻る。
【００６７】
ここで、ステップＳ１４０の熱流束推定処理の内容について詳述する。ステップＳ１４０では、エアコンＥＣＵ１４は、図４に示す処理を実行する。
【００６８】
図４の処理においてエアコンＥＣＵ１４は、まずステップＳ１４１で、現在の時刻に基づいて、着座時点を起点（すなわちｔ＝１）とする経過時間ｔを特定する。より具体的には、直前のステップＳ１２５を実行して以降初めてのステップＳ１４１の実行機会である場合は、ステップＳ１２５で記録した着座時点から所定の待機時間すなわち１秒が経過するまで待機する。そして、当該待機時間が経過した時点で、着座時点に上記待機時間の長さ分だけ加算した時間（すなわちｔ＝２）を、新たな経過時間とする。
【００６９】
また、直前のステップＳ１２５を実行して以降２回目以上のステップＳ１４１の実行機会である場合は、前回のステップＳ１４１の実行機会から上記待機時間が経過するまで待機する。そして、当該待機時間が経過した時点で、前回のステップＳ１４１で特定した経過時間に上記待機時間の長さ分だけ加算した時間を、新たな経過時間とする。例えば、前回のステップＳ１４１で特定した経過時間ｔが１００であった場合は、ｔ＝１０１を新たな経過時間とする。
【００７０】
続いてエアコンＥＣＵ１４は、ステップＳ１４２で、直前で算出した経過時間ｔにおける対象座席の対象部分と着座者の対象人体部位との間の熱流束を第１計算式Σ１で算出する。そして続くステップＳ１４３で、直前で算出した経過時間ｔにおける対象座席の対象部分と着座者の対象人体部位との間の熱流束を第２計算式Σ２で算出する。
【００７１】
第１計算式Σ１および第２計算式Σ２は、以下の通りである。
Σ１：Ｒｔ＝Ｒｉｎｉｔ−ω×ｌｎ（ｔ）＋ε×（ｌｎ（ｔ）−ｌｎ（ｔ−１））
Σ２：Ｒｔ＝Ｒｔ−１−（ｌｎ（ｔ）−ｌｎ（ｔ−１））×Ｒｔ−１＋ε×（ｌｎ（ｔ）
−ｌｎ（ｔ−１））
ここで、Ｒｔは経過時間ｔに対応する時点における算出対象の熱流束であり、その単位は［Ｗ／ｍ２］である。また、Ｒｔ−１は前回の経過時間ｔ−１に対応する時点における算出対象の熱流束である。またｌｎ（）は、括弧内の値の自然対数である。
【００７２】
また、Ｒｉｎｉｔは着座時点ｔ＝１における算出対象の熱流束であり、以下の式によって算出される。
Ｒｉｎｉｔ＝ｋ１×（Ｔｓｉ−ａ１）
ここで、ｋ１、ａ１は、算出される熱流束が実際の熱流束に近くなるよう、実験等によってあらかじめ決められた正の定数である。ａ１は、人の体温に近い温度である。以下、ａ１を第１基準温度という。
【００７３】
また、ωは以下の式によって算出される。
ω＝ｋ２×（Ｔｓｉ−ａ２）
ここで、ｋ２、ａ２は、算出される熱流束が実際の熱流束に近くなるよう、実験等によってあらかじめ決められた正の定数である。ａ２は、人の体温に近い温度である。以下、ａ２を第１基準温度という。ｋ２はｋ１よりも小さい。例えば、ｋ２はｋ１の約１／６である。
【００７４】
εは、ヒータクーラ５−８のうち対象のヒータクーラの現時点の作動状態に応じて決定される。対象のヒータクーラは、ヒータクーラ５−８のうち、主として対象の座席の対象の部分の温度調整を行うヒータクーラである。εの詳細については後述する。
【００７５】
ステップＳ１４３に続いては、エアコンＥＣＵ１４は、ステップＳ１４５に進み、第１計算式Σ１による熱流束Ｒｔの算出値と第２計算式Σ２による熱流束Ｒｔの算出値とが同等であるか否かを判定する。両方の算出値の差の絶対値が所定の許容値よりも小さい場合に同等であると判定され、当該許容値以上の場合は同等でないと判定される。
【００７６】
同等であると判定した場合、エアコンＥＣＵ１４はステップＳ１４６に進んで切替フラグをオンにして、ステップＳ１４７に進む。この切替フラグは、図３の処理毎にエアコンＥＣＵ１４のＲＡＭに設けられるフラグであり、図３の処理の開始時にオフに初期化される。ステップＳ１４５で同等でないと判定した場合、エアコンＥＣＵ１４はステップＳ１４６をバイパスしてステップＳ１４７に進む。
【００７７】
ステップＳ１４７では、切替フラグがオンであるか否かを判定し、オフであればステップＳ１４８に進み、オンであればステップＳ１４９に進む。ステップＳ１４８では、第１計算式Σ１による熱流束Ｒｔの算出値と第２計算式Σ２による熱流束Ｒｔの算出値のうち、前者を現在の熱流束として採用し、ステップＳ１５０に進む。ステップＳ１４９では、第１計算式Σ１による熱流束Ｒｔの算出値と第２計算式Σ２による熱流束Ｒｔの算出値のうち、後者を現在の熱流束として採用し、ステップＳ１５０に進む。
【００７８】
このような図４の処理により、経過時間ｔ＝２およびその後において、第１計算式Σ１による熱流束Ｒｔの算出値と第１計算式Σ１による熱流束Ｒｔの算出値とが同等であると判定されるまでは、以下のような作動になる。すなわち、経過時間ｔが待機時間（すなわち１秒）だけ増える度に、切替フラグがオフのままステップＳ１４１、Ｓ１４２、Ｓ１４３、Ｓ１４５、Ｓ１４７、Ｓ１４８が実行される。つまり、第１計算式Σ１による熱流束Ｒｔの算出値が現時点の熱流束として採用される。
【００７９】
その後、第１計算式Σ１による熱流束Ｒｔの算出値と第１計算式Σ１による熱流束Ｒｔの算出値とが同等になったと判定されると、ステップＳ１４１、Ｓ１４２、Ｓ１４３、Ｓ１４５に続いてステップＳ１４６が実行されて切替フラグがオンになる。すると、続いてステップＳ１４７、Ｓ１４９が実行される。すなわち、第２計算式Σ２による熱流束Ｒｔの算出値が現時点の熱流束として採用される。この時点が、切替タイミングに相当する。
【００８０】
さらにその後は、ステップＳ１４５の判定結果にかかわらず切替フラグがオンのままなので、経過時間ｔが待機時間だけ増える度に、ステップＳ１４１、Ｓ１４２、Ｓ１４３、Ｓ１４５、Ｓ１４７、Ｓ１４９が実行される。すなわち、第２計算式Σ２による熱流束Ｒｔの算出値が現時点の熱流束として採用される。
【００８１】
ここで、第１計算式Σ１、第２計算式Σ２によって算出される熱流束Ｒｔの特徴について、図５を参照して説明する。図５は、着座時点ｔ＝１における熱流束が負の場合、すなわち、対象の人体部位から対象の座席の対象の部分へと熱が流れる場合の、その後に算出される熱流束の推移の一例を表す。縦軸が熱流束であり、横軸が経過時間ｔである。経過時間ｔ＝ｔｃが、切替タイミングに相当する。図５において、第１計算式Σ１で算出された熱流束Ｒｔが実線５１で示され、第２計算式Σ２で算出された熱流束Ｒｔが一点鎖線５２で示される。
【００８２】
まず第１計算式Σ１について説明する。第１計算式Σ１の右辺は、３つの項から構成されている。第１項は、Ｒｉｎｉｔである。Ｒｉｎｉｔは、対象の座席に着座者が着座した時点すなわち着座時点ｔ＝１における、対象の座席の対象の部分と対象の人体部位との間の熱流束である。すなわち、Ｒｉｎｉｔは、着座時熱流束である。
【００８３】
上述の通り、Ｒｉｎｉｔ＝ｋ１×（Ｔｓｉ−ａ１）という式が成り立つので、着座時熱流束Ｒｉｎｉｔは初期温度Ｔｓｉによって決まる。これは、ある時点における２つの物体の間の熱流束はその時点における２つの物体の温度差に比例するという、熱伝達に関する周知のフーリエの法則に適っている。このように、着座地熱流束は、第１計算式Σ１および第２計算式Σ２による計算の過程で、初期温度Ｔｓｉに基づいて推定される。
【００８４】
第１計算式Σ１の右辺を構成する第２項は、−ω×ｌｎ（ｔ）である。この第２項は、因数としてωとｌｎ（ｔ）を含んでいる。ｌｎ（ｔ）は、着座時点を起点とする経過時間ｔに依存して変化すると共に初
期温度Ｔｓｉに依存しない時間依存要素である。ωは、初期温度Ｔｓｉに依存すると共に経過時間ｔに依存しない初期温度依存要素である。
【００８５】
ｌｎ（ｔ）はｔ＝２以降は常に正なので、−ω×ｌｎ（ｔ）の正負は、ω＝ｋ２×（Ｔｓｉ−ａ２）の正負によって決まる。また、ωの正負は、Ｒｉｎｉｔ＝ｋ１×（Ｔｓｉ−ａ１）の正負と殆どの場合で同じである。これは、人の体温に相当する第１基準温度ａ１と第２基準温度ａ２がほぼ同じ値だからである。
【００８６】
したがって、殆どの場合で、Ｒｉｎｉｔが正の場合は−ω×ｌｎ（ｔ）は負となり、Ｒｉｎｉｔが負の場合は−ω×ｌｎ（ｔ）は正となる。また、ｋ１はｋ２よりも十分大きい。また、ｌｎ（ｔ）は経過時間ｔの増加関数である。したがって、第１計算式Σ１の右辺を構成する第１項と第２項すなわちＲｉｎｉｔ−ω×ｌｎ（ｔ）は、経過時間ｔが増大するほど、ｔ＝１における熱流束Ｒｔ＝Ｒｉｎｉｔに比べて小さくなる。
【００８７】
つまり、時間依存要素であるｌｎ（ｔ）は、ある期間においては、経過時間ｔが長くなるほど熱流束Ｒｔの絶対値が小さくなるよう、第１計算式Σ１に寄与する。以下、この期間を減少期間という。図５の例においては、ｔ＝２から熱流束Ｒｔがゼロになるまでの期間Ｘ１が減少期間に相当する。
【００８８】
別の観点から見れば、時間依存要素であるｌｎ（ｔ）は、ｔ＝２以降、経過時間ｔが長くなるほど熱流束Ｒｔが着座時熱流束Ｒｉｎｉｔからゼロ側に大きく離れるよう、第１計算式Σ１に寄与する。実際、図５においても、ｔ＝２以降、経過時間ｔが長くなるほど熱流束ＲｔがＲｉｎｉｔからゼロ側に（すなわち、熱流束が上昇する側に）大きく離れるよう、第１計算式Σ１に寄与する。時間依存要素であるｌｎ（ｔ）がこのように機能することで、第１計算式Σ１によって推定された熱流束Ｒｔの、経過時間ｔの変化に応じた推移が、実際の熱流束に近くなる。
【００８９】
また、時間依存要素であるｌｎ（ｔ）の時間微分１／ｔは、経過時間ｔの減少関数である。したがって、ｌｎ（ｔ）は、経過時間ｔが長くなるほど熱流束の変化速度の絶対値が小さくなるよう、第１計算式Σ１に寄与する。これは、減少期間Ｘ１においても、それ以外の期間においても、同様である。時間依存要素であるｌｎ（ｔ）がこのように機能することで、第１計算式Σ１によって推定された熱流束Ｒｔの、経過時間ｔの変化に応じた推移が、さらに実際の熱流束に近くなる。
【００９０】
また、初期温度依存要素であるωの絶対値は、初期温度Ｔｓｉと第１基準温度ａ１との差の絶対値が大きいほど、大きくなる。またωは、時間依存要素であるｌｎ（ｔ）と共に、第１計算式Σ１の右辺を構成する第２項の因数となる。したがって、ωは、初期温度Ｔｓｉと第１基準温度ａ１との差の絶対値が大きいほど、第１計算式Σ１に対する時間依存要素ｌｎ（ｔ）の寄与を増幅させる。このようにすることで、経過時間ｔ＝２以降の初期における、第１計算式Σ１によって推定された熱流束Ｒｔの、経過時間ｔの変化に応じた変化速度が、初期温度Ｔｓｉに応じて異なることで、算出された熱流束が実際の熱流束に近くなる。
【００９１】
また、第１計算式Σ１に現れる対象の座席の対象の部分の表面温度の推定値は、着座時点の表面温度である初期温度Ｔｓｉのみであり、当該部分の他の時点の表面温度の推定値は現れない。このように、当該対象の部分の表面温度を着座後も常時測定する必要がないので、着座後の表面温度測定による誤差が発生しない。したがって、熱流束の推定が安定化する。
【００９２】
別の観点から言えば、第１計算式Σ１に対する初期温度依存要素の寄与度は、第１計算式Σ１に対する着座時点以外の時点の当該対象の部分の表面温度から成る要素（例えば項）の寄与度よりも高い。なぜなら、本実施形態においては、着座時点以外の時点の当該対象の部分の表面温度から成る項の係数がゼロだからである。このように、当該対象の部分の着座後の表面温度の寄与が小さいので、着座後の測定によるが発生する可能性が低い。したがって、熱流束の推定が安定化する。
【００９３】
第１計算式Σ１の右辺を構成する第３項は、ε×（ｌｎ（ｔ）−ｌｎ（ｔ−１））である。この第３項は、因数として温調依存要素であるεと時間依存要素である（ｌｎ（ｔ）−ｌｎ（ｔ−１））を含む。
【００９４】
ここで温調依存要素εは、上述の通り、ヒータクーラ５−８のうち対象のヒータクーラの現時点の作動状態に応じて決定される。具体的には、図６に示すように、対象のヒータクーラが非作動の場合は、εはゼロになる。したがって、εを含む第３項の熱流束への寄与はゼロになる。また、対象のヒータクーラが低出力冷却モードで作動中の場合はεが−Ｃ１になり、高出力冷却モードで作動中の場合はεが−Ｃ２になる。ここで、Ｃ１、Ｃ２は正の値であり、Ｃ１＜Ｃ２の関係が成立する。つまり、低出力冷却モードよりも高出力冷却モードの方が、εを含む第３項の熱流束への寄与が大きくなる。
【００９５】
また、対象のヒータクーラが低出力加熱モードで作動中の場合はεがＨ１になり、高出力加熱モードで作動中の場合はεがＨ２になる。ここで、Ｈ１、Ｈ２は正の値であり、Ｈ１＜Ｈ２の関係が成立する。つまり、低出力加熱モードよりも高出力加熱モードの方が、εを含む第３項の熱流束への寄与が大きくなる。
【００９６】
このように、第３項が対象のヒータクーラの作動状態に依存する温調依存要素εを含むことで、対象のヒータクーラの作動状態（すなわち能力）に起因する熱流束の変化が第１計算式Σ１に反映される。上述の通り、算出された熱流束に基づいて対象のヒータクーラの能力が調整されるので、対象のヒータクーラのフィードバック制御が実現する。
【００９７】
また、時間依存要素である（ｌｎ（ｔ）−ｌｎ（ｔ−１））は、ｔに関する減少関数であり、ｔ＝２以上の場合は常に正であり、ｔが増大するほどゼロに漸近的に近付く。したがって、第３項の時間依存要素は、第１計算式Σ１で算出される熱流束Ｒｔに対するヒータクーラの作動状態の寄与度を、時間の経過と共に低下させていく。このようにすることで、対象のヒータクーラの加熱対象の人体部位が、着座時点後の初期の段階において、時間経過と共にヒータクーラの表面温度に急速に近付いていく現象を反映した熱流束計算を、エアコンＥＣＵ１４が行うことができる。
【００９８】
次に第２計算式Σ２について説明する。第２計算式Σ２の右辺は、３つの項から構成されている。第１項は、Ｒｔ−１である。Ｒｔ−１は、前回のステップＳ１４３において第２計算式Σ２によって算出された熱流束である。つまり、Ｒｔ−１は、現時点よりも１秒前の熱流束の算出値であると共に、現時点よりも１回前の熱流束の算出値である。したがって、第２計算式Σ２の右辺の第２項および第３項は、現時点よりも１回前の熱流束の算出値Ｒｔ−１から現時点で算出する熱流束Ｒｔへの変化量に相当する。
【００９９】
第２計算式Σ２の右辺の第２項は、−（ｌｎ（ｔ）−ｌｎ（ｔ−１））×Ｒｔ−１である。この第２項には、−１、（ｌｎ（ｔ）−ｌｎ（ｔ−１））、Ｒｔ−１という因数が含まれている。これらのうち、時間依存要素である（ｌｎ（ｔ）−ｌｎ（ｔ−１））は経過時間ｔ＝２以降では正である。したがって、第２項は、今回算出する熱流束Ｒｔを前回算出の熱流束Ｒｔ−１よりもゼロ側に変動させるよう第２計算式Σ２に寄与する。そして、時間依存要素である（ｌｎ（ｔ）−ｌｎ（ｔ−１））は経過時間ｔ＝２以降ではｔが増大するほど漸近的にゼロに近付く。したがって、前回算出の熱流束Ｒｔ−１に対して今回算出する熱流束Ｒｔを第２項が変化させる量の絶対値は、ｔが増大するほど小さくなっていき、漸近的にゼロに近付く。
【０１００】
そして、ｔ＝２以降において、第１項と第２項の和におけるＲｔ−１の係数は、１より小さい正数である。したがって、第１項と第２項は、Ｒｔの絶対値を時間経過と共に減少させていくよう、第２計算式Σ２に寄与する。
【０１０１】
第２計算式Σ２の右辺の第３項は、ε×（ｌｎ（ｔ）−ｌｎ（ｔ−１））であり、第１計算式Σ１の右辺を構成する第３項と同じである。したがって、第２計算式Σ２の右辺の第３項自体の振る舞いは第１計算式Σ１の右辺を構成する第３項と同じである。なお、第２計算式Σ２の右辺の第３項は、現時点よりも１回前の熱流束の算出値Ｒｔ−１から現時点で算出する熱流束Ｒｔへの変化量として寄与する点は、第１計算式Σ１の右辺を構成する第３項とは異なる。
【０１０２】
第２項も第３項も（ｌｎ（ｔ）−ｌｎ（ｔ−１））という因数を含んでいる。したがって、第２計算式Σ２は、時間経過と共に前記熱流束Ｒｔの経時変化量の絶対値を減少させる。
【０１０３】
なお、この第２計算式Σ２はｔ＝２以降で用いられるが、ｔ＝１における熱流束は、第１計算式Σ１と同様に、Ｒｉｎｉｔが採用される。第２計算式Σ２の特徴は、どの時点の表面温度も要素として含まれないことである。
【０１０４】
殆どの場合、図５に例示するように、ｔ＝２からｔ＝ｔｃ−１までの期間では、熱流束Ｒｔの着座時熱流束Ｒｉｎｉｔに対するゼロ側への乖離量（すなわち差分）は、第１計算式Σ１で算出された値５１よりも第２計算式Σ２で算出された値５２の方が、大きい。そして、発明者の実験によれば、第１計算式Σ１で算出された値５１の方が、実際に発生している熱流束に近い。
【０１０５】
そして、殆どの場合、図５に示すように、ｔ＝ｔｃの時点において、ｔ＝２以降で初めて値５１、５２の差の絶対値が所定の許容値よりも小さくなる。例えば、ｔ＝ｔｃの時点において値５１、５２が一致する。ｔ＝２以降で初めて値５１、５２の差の絶対値が所定の許容値よりも小さくなるタイミングは、初期温度Ｔｓｉに依存して変動する。
【０１０６】
そして、殆どの場合、図５に例示するように、ｔ＝ｔｃよりも後の期間では、熱流束Ｒｔの着座時熱流束Ｒｉｎｉｔに対するゼロ側への乖離量（すなわち差分）は、第１計算式Σ１で算出された値５１よりも第２計算式Σ２で算出された値５２の方が、小さい。そして、発明者の実験によれば、第２計算式Σ２で算出された値５２の方が、実際に発生している熱流束に近い。
【０１０７】
したがって、上述の通り、ｔ＝２以降かつｔ＝ｔｃの直前までは、第１計算式Σ１による熱流束Ｒｔの算出値５１が現時点の熱流束として採用され、ｔ＝ｔｃ以降は、第２計算式Σ２による熱流束Ｒｔの算出値５２が現時点の熱流束として採用される。このようになっていることで、第１計算式Σ１のみで熱流束を算出する場合や、第２計算式Σ２のみで熱流束を算出する場合に比べて、精度の高い熱流束算出を行うことができる。すなわち、前期は第１計算式Σ１によって初期温度に応じた熱流束を算出できる。そして、後期は座席の温度にかかわらず時間経過と共に経時変化量の絶対値が低下していく熱流束の特性を第２計算式Σ２で再現することができる。
【０１０８】
ここで、所定の条件において上述の通りにエアコンＥＣＵ１４で算出された熱流束と、熱流束センサで実際に測定された熱流束との比較結果を、図７、図８に示す。図７は、初期温度Ｔｓｉが人の体温よりも低い１７℃の場合の結果を、図８は、初期温度Ｔｓｉが人の体温よりも高い４５℃の場合の結果を、示す。両図において、実線がエアコンＥＣＵ１４によって上述の通りに算出された熱流束を示し、破線が熱流束センサで実際に測定された熱流束を示す。
【０１０９】
図７の例では、Ｒｉｎｉｔもωも負となるので、ｔ＝２から時間経過と共に熱流束は上昇していく。図８の例では、Ｒｉｎｉｔもωも正となるので、ｔ＝２から時間経過と共に熱流束は低下していく。いずれの場合も、熱流束センサで計測された熱流束に近い値が算出されている。
【０１１０】
以上説明した通り、エアコンＥＣＵ１４は、赤外線カメラ４の出力に基づいて、座席２、３の背もたれ部および座部の各々について、人が当該座席に着座する際の当該座席の当該部分の温度である初期温度Ｔｓｉを推定する。そしてエアコンＥＣＵ１４は、当該人が当該座席に着座した着座時点ｔ＝１の後の期間における、当該座席の当該部分と当該人との間の熱流束Ｒｔの推移（すなわち経時変化）を算出する。
【０１１１】
そしてエアコンＥＣＵ１４は、当該推移を算出する際、時間依存要素ｌｎ（ｔ）と初期温度依存要素ωとを含む第１計算式Σ１により、ｔ＝２からｔ＝ｔｃまでの期間における熱流束の推移を算出する。
【０１１２】
このように、エアコンＥＣＵ１４は、初期温度Ｔｓｉと着座時点ｔ＝１を起点とする経過時間に基づいて、その後の熱流束を算出する。したがって、結果的に初期温度Ｔｓｉが推定できるセンサを用いれば、熱流束センサを用いる必要なく、座席と着座者の間の熱の流れを推定することができる。
【０１１３】
従来、車両の座席の比較的表面側に埋め込まれた熱流束センサによって、座席と着座者の間の熱の流れを測定する技術が知られている。しかしながら、発明者の検討によれば、座席と着座者の間の熱の流れを推定するために、その用途に特化した機能を有する熱流束センサを利用すると、そのセンサを他の用途に転用することが困難になる。したがって、他の用途への転用の余地が熱流束センサよりも広いセンサを利用して、座席と着座者の間の熱の流れを推定できることが望ましい。したがって、本実施形態のような第１計算式Σ１を用いた熱流束の推定が有用である。
【０１１４】
また、赤外線カメラ４は、上述の通り、乗員の体表面の温度を特定する用途にも用いられるので、熱流束の推定以外の用途に転用することが容易なセンサである。また、通常、車両１に搭載される赤外線カメラ４は、座席に座った乗員の前側の表面の温度を特定するために用いられるのに対し、本実施形態では、乗員の背中または臀部と座席との間の、すなわち隠れた部分の、熱流束を推定するために用いられる。つまり、本実施形態における赤外線カメラ４の用途は、従来の用途とは観点が大きく異なっている。
【０１１５】
また、熱流束センサは非常に高価であるという問題がある。しかも、熱流束センサを使用する場合、運転席２の背もたれ部および座部、および、後席３の背もたれ部および座部の各々に熱流束センサを配置する必要があるので、部品点数の増大を招くと共に、高価であるということの問題がより顕著になる。これに対し、本実施形態で用いられる赤外線カメラ４は、１台で複数の座席の複数の部分の温度を検出可能なので、熱流束の測定対象が多い場合でも部品点数が低減されると共に設置費用を抑えることができる。また、本実施形態では、熱流束を推定するために、着座者の体温を検出することは行っていない。
【０１１６】
なお、本実施形態においては、エアコンＥＣＵ１４が、ステップＳ１３０を実行することで、温度推定部として機能し、ステップＳ１４０を実行することで熱流算出部として機能する。
【０１１７】
（第２実施形態）
次に第２実施形態について説明する。本実施形態は、第１実施形態に対して、エアコンＥＣＵ１４が実行する熱流速推定処理１４Ａの処理内容が異なり、その他は第１実施形態と同じである。
【０１１８】
具体的には、エアコンＥＣＵ１４は、車室内温度調整処理１４Ｂにおいて、図３の処理に代えて図９の処理を実行する。図９の処理では、図３の処理に対して、ステップＳ１１０が廃され、ステップＳ１３０に代えてステップＳ１３５が設けられている。
【０１１９】
したがって、本実施形態では、エアコンＥＣＵ１４は、１つの図９の処理の実行において、ステップＳ１２０で対象の座席に人が着座していないとの判定を繰り返す間は、対象の座席の対象の部分の表面温度を一度も特定していない。ステップＳ１２０における判定方法は、図３の処理と同じである。
【０１２０】
そして、エアコンＥＣＵ１４は、ステップＳ１２０で対象の座席に人が着座していると判定してステップＳ１２５を実行する。ステップＳ１２５では、図３の処理と同様、着座時点の時刻をＲＡＭに記録するが、記録する着座時点の特定方法が、図３の処理とは異なる。
【０１２１】
具体的には、赤外線カメラ４およびエアコンＥＣＵ１４が常時作動であるかオン時限定作動であるかによって、記録する着座時点の特定方法が異なる。常時作動の場合は、対象の座席に人が着座していない時点においてステップＳ１２０の処理が繰り返されているので、ステップＳ１２５では、直前のステップＳ１２０で着座したと判定した時刻を着座時点として記録する。
【０１２２】
オン時限定作動の場合は、エアコンＥＣＵ１４が図９の処理を開始した時点で既に人が対象の座席に着座している場合がある。特に、対象の座席が運転席２である場合は、その可能性が非常に高い。その場合、図９の処理が始まって最初のステップＳ１２０で着座していると判定されるが、実際の着座時点は、ステップＳ１２０で着座していると判定された時点ではなく、もっと前の時点（例えば１分前）である可能性がある。
【０１２３】
そのような場合のために、オン時限定作動におけるエアコンＥＣＵ１４は、ステップＳ１２５では、図９の処理が始まって最初のステップＳ１２０で着座していると判定している場合は、他の情報に基づいて着座時点を特定してもよい。
【０１２４】
例えば、車両１のすべてのドアの各々に、当該ドアの開閉を検出するドア開閉センサが設置されている場合は、これらのドア開閉センサを利用して着座時点を特定してもよい。ただしこの場合、これらドア開閉センサは常時作動であり、また、これらドア開閉センサの検出結果および時刻を逐次記録する不図示のドア開閉記録装置も常時作動である。
【０１２５】
この場合、エアコンＥＣＵ１４は、ステップＳ１２５では、ドア開閉記録装置の記録に基づいて、ドアが開いてから所定時間後（例えば３秒後）の時刻を着座時点として記録してもよい。
【０１２６】
なお、オン時限定作動であっても、対象の座席が後席３である場合は、エアコンＥＣＵ１４が図９の処理を開始した時点でまだ人が対象の座席に着座していない場合も多い。
【０１２７】
そのような場合のために、オン時限定作動におけるエアコンＥＣＵ１４は、ステップＳ１２５では、図９の処理が始まって最初のステップＳ１２０で着座していないと判定していた場合は、ステップＳ１２５では第１実施形態と同様の判定を行ってもよい。すなわち、直前のステップＳ１２０で着座したと判定した時刻を着座時点として記録してもよい。
【０１２８】
エアコンＥＣＵ１４は、ステップＳ１２５の後、ステップＳ１３５で、着座の際の対象の座席の対象の部分の表面温度を推定し、推定の結果得られた推定値を初期温度Ｔｓｉとする。ステップＳ１３５に続いては、ステップＳ１４０に進む。ステップＳ１４０以降の処理内容は、図３の処理と同じである。
【０１２９】
ここで、ステップＳ１３５における表面温度の推定について説明する。ステップＳ１３５のタイミングでは、対象の座席に着座者が着座しているのだから、対象の座席の対象の部分の大部分は、赤外線カメラ４に対して着座者の後側に隠れている。
【０１３０】
したがって、エアコンＥＣＵ１４は、車両１のすべての座席（すなわち、運転席２、助手席、後席３、およびその他の座席）のうち、まだ人が着座していない座席を特定する。そして、エアコンＥＣＵ１４は、特定した座席の対象の部分であるとしてあらかじめ定められた画素範囲内における複数の画素値を抽出する。そして、抽出した複数の画素値の統計代表値（例えば平均値）を、対象の座席への着座の際の、対象の座席の対象の部分の表
面温度の推定値とする。
【０１３１】
なお、本実施形態では、車両１の運転席２、後席３以外のすべての座席にも、着座センサ１９、２０と同様の着座センサが取り付けられており、エアコンＥＣＵ１４は、それら着座センサからの信号に基づいて、各座席における着座の有無を判定できる。
【０１３２】
ステップＳ１３５では、対象の座席に人が着座した直後に実行される場合もあれば、対象の座席に人が着座してからしばらく時間が経ってから実行される場合もある。しかし、どちらの場合であっても、人が座っていない座席の温度は、対象の座席への着座の際における対象の座席の表面温度に近い温度となっている。したがって、どちらの場合にステップＳ１３５が実行されたとしても、得られた推定値は、対象の座席への着座の際の、対象の座席の対象の部分の表面温度の推定値として有効である。
【０１３３】
また、エアコンＥＣＵ１４は、車両１のすべての座席に人が着座している場合、サーモグラフィのうち、人が映らない部分としてあらかじめ定められた画素範囲における複数の画素値の統計代表値を、対象の座席の対象の部分の表面温度の推定値とする。
【０１３４】
具体的には、エアコンＥＣＵ１４は、対象の座席の対象の部分のうち、人が着座していても赤外線カメラ４が撮影できる可能性が高いとしてあらかじめ定められた画素範囲における複数の画素値を抽出する。そして抽出した複数の画素値の統計代表値（例えば平均値）を、対象の座席の対象の部分の表面温度の推定値とする。そのような画素範囲としては、例えば、対象の座席の対象の部分のうち車両左右方向に面している側面部分、対象の座席のヘッドレストの側面部分、対象の座席の背もたれ部うち車両上方向に面している天面部分がある。また、そのような画素範囲としては、車両のセンターコンソール等の、座席でない部分を採用してもよい。
【０１３５】
これらの部分の温度も、人が座っていない座席の温度と同様の理由で、得られた推定値は、対象の座席への着座の際の、対象の座席の対象の部分の表面温度の推定値として有効である。
【０１３６】
なお、エアコンＥＣＵ１４は、サーモグラフィのうち、人が映らない部分としてあらかじめ定められた画素範囲における複数の画素値の統計代表値そのものを、対象の座席の対象の部分の表面温度の推定値としなくてもよい場合がある。すなわち、当該人が映らない部分の表面温度と対象の座席の対象の部分の表面温度との対応関係があらかじめ実験等により得られてエアコンＥＣＵ１４のＲＯＭに記録されている場合がある。そのような場合、エアコンＥＣＵ１４は、当該統計代表値を当該対応関係に適用することで、対象の座席の対象の部分の表面温度の推定値を得てもよい。
【０１３７】
このようにすることで、対象の座席の対象の部分の表面温度を直接測定できなくても、車両１内の他の部分の表面温度を用いることで、ある程度の正確性をもって初期温度Ｔｓｉを算出することができる。また、本実施形態では、熱流束を推定するために、着座者の体温を検出することは行っていない。
【０１３８】
なお、本実施形態においては、エアコンＥＣＵ１４が、ステップＳ１３５を実行することで、温度推定部として機能し、ステップＳ１４０を実行することで熱流算出部として機能する。
【０１３９】
（他の実施形態）
本発明は上記した実施形態に限定されるものではなく、適宜変更が可能である。また、上記各実施形態は、互いに無関係なものではなく、組み合わせが明らかに不可な場合を除き、適宜組み合わせが可能である。また、上記各実施形態において、実施形態を構成する要素は、特に必須であると明示した場合および原理的に明らかに必須であると考えられる場合等を除き、必ずしも必須のものではない。また、上記各実施形態において、実施形態の構成要素の個数、数値、量、範囲等の数値が言及されている場合、特に必須であると明示した場合および原理的に明らかに特定の数に限定される場合等を除き、その特定の数に限定されるものではない。また、上記実施形態において、センサから車両の外部環境情報（例えば車外の湿度）を取得することが記載されている場合、そのセンサを廃し、車両の外部のサーバまたはクラウドからその外部環境情報を受信することも可能である。あるいは、そのセンサを廃し、車両の外部のサーバまたはクラウドからその外部環境情報に関連する関連情報を取得し、取得した関連情報からその外部環境情報を推定することも可能である。特に、ある量について複数個の値が例示されている場合、特に別記した場合および原理的に明らかに不可能な場合を除き、それら複数個の値の間の値を採用することも可能である。また、上記各実施形態において、構成要素等の形状、位置関係等に言及するときは、特に明示した場合および原理的に特定の形状、位置関係等に限定される場合等を除き、その形状、位置関係等に限定されるものではない。また、本発明は、上記各実施形態に対する以下のような変形例および均等範囲の変形例も許容される。なお、以下の変形例は、それぞれ独立に、上記実施形態に適用および不適用を選択できる。すなわち、以下の変形例のうち任意の組み合わせを、上記実施形態に適用することができる。
【０１４０】
（変形例１）
上記実施形態では、運転席２、後席３について、背もたれ部および座部と着座者の間の熱流束の算出が行われている。しかし、車両１の運転席２、後席３以外の座席についても、同様に、背もたれ部および座部と着座者の間の熱流束の算出が行われてもよい。また、１つの座席について、右大腿部、左大腿部、臀部、背中のそれぞれについて、熱流束の算出が行われてもよい。
【０１４１】
（変形例２）
上記実施形態では、赤外線カメラ４の撮影範囲は固定となっていたが、必ずしもこのようになっておらずともよい。例えば、赤外線カメラ４は、あるタイミングでは運転席２を撮影し、別のタイミングでは後席３を撮影するという様に、撮影方向を周期的に変化させるようになっていてもよい。
【０１４２】
（変形例３）
上記実施形態では、ヒータクーラ５−８は非作動、高出力冷却モード、低出力冷却モード、低出力加熱モード、高出力加熱モードの５つの状態間を遷移するが、非作動、冷却モード、加熱モードの３つの状態のみを遷移してもよい。
【０１４３】
また、ヒータクーラ５−８は、加熱能力を有し冷却能力を有さないヒータに置き換えられてもよい。また、ヒータクーラ５−８は、冷却能力を有し加熱能力を有さないクーラに置き換えられてもよい。
【０１４４】
また、上記各実施形態において、ヒータクーラ５−８が廃されてもよい。その場合、エアコンＥＣＵ１４によって算出された熱流束は、ヒータクーラ５−８の制御以外の用途（例えば、車室内温度調整処理１４Ｂ）に利用されてもよいし、それ以外の用途（例えば単なる記録用途）に用いられてもよい。
【０１４５】
（変形例４）
上記実施形態において、赤外線カメラ４は、運転席２、後席３および車両１内の他の座席のすべてを撮影範囲としていたが、必ずしもこのようになっておらずともよい。例えば、赤外線カメラ４は、車両１内の座席のうち後席３のみを撮影範囲としてもよい。そのような撮影範囲を実現するための赤外線カメラ４の配置としては、例えば、車両１の車室内における前後方向中央かつ左右方向中央の天井部がある。この場合、エアコンＥＣＵ１４は、後席３の背もたれ部と着座者の背中の間の熱流束、および、後席３の座部と着座者の臀部の間の熱流束のみを算出してもよい。あるいは、エアコンＥＣＵ１４は、後席３の背もたれ部と着座者の背中の間の熱流束のみを算出してもよい。
【０１４６】
（変形例５）
上記実施形態では、ステップＳ１２０における対象の座席に人が着座したか否かの判定は、着座センサ１９、２０の検出結果に基づいて行われている。しかし、この判定は、必ずしも着座センサの検出結果に基づいて行われなくてもよい。すなわち、運転席着座センサ１９、後席着座センサ２０は必須ではない。
【０１４７】
例えば、車両１のすべてのドアの各々に、当該ドアの開閉を検出するドア開閉センサが設置されている場合は、これらのドア開閉センサを利用して上記判定を行うことが可能である。すなわち、エアコンＥＣＵ１４は、ステップＳ１２０においては、対象の座席に最も近いドアとしてあらかじめ定められたドアのドア開閉センサの検出結果に基づいて、当該ドアが開いてから所定時間が経過した時点で、人が当該座席に着座していると判定してもよい。所定時間は、例えば３秒でもよい。
【０１４８】
また例えば、車両１のすべての座席の各々に、当該座席におけるシートベルトの装着の有無を検出するシートベルト装着センサが設置されている場合は、これらのシートベルト装着センサを利用して上記判定を行うことが可能である。すなわち、エアコンＥＣＵ１４は、ステップＳ１２０においては、対象の座席のシートベルト装着センサの検出結果に基づいて、当該座席のシートベルトが装着された時点よりも所定時間前の時点に、人が当該座席に着座したと判定してもよい。所定時間は、例えば５秒でもよい。
【０１４９】
この場合は、図３のステップＳ１２５で記録する着座時点は、当該座席のシートベルトが装着された時点よりも当該所定時間前の時点の時刻となる。また、図３のステップＳ１３０で特定される初期温度Ｔｓｉは、当該着座時点に最も近い時点においてステップＳ１１０で特定された対象の座席の対象の部分の表面温度となる。
【０１５０】
また例えば、赤外線カメラ４が出力するサーモグラフィを利用して上記判定を行うことが可能である。すなわち、エアコンＥＣＵ１４は、対象の座席に該当するとあらかじめ定められたサーモグラフィ中の画素範囲の複数の画素の画素値の変化に基づいて、当該座席に人が着座したか否かを判定してもよい。例えば、当該複数の画素の統計代表値（例えば平均値）が基準時間内に基準量以上変化したときに、当該座席に人が着座したと判定してもよい。
【０１５１】
また例えば、可視光で車室内を撮影する可視光カメラが車両１に設置されている場合、当該可視光カメラが出力する撮影画像を利用して上記判定を行うことが可能である。すなわち、エアコンＥＣＵ１４は、対象の座席に該当するとあらかじめ定められた撮影画像中の画素範囲の複数の画素の画素値（すなわち輝度）の変化に基づいて、当該座席に人が着座したか否かを判定してもよい。例えば、当該複数の画素の統計代表値（例えば平均値）が基準時間内に基準量以上変化したときに、当該座席に人が着座したと判定してもよい。      
【０１５２】
（変形例６）
上記第１実施形態では、対象の座席の対象の部分の初期温度Ｔｓｉは、当該座席に人が着座する直前に当該部分について検出された表面温度とされている。また、上記第２実施形態では、対象の座席の対象の部分の初期温度Ｔｓｉは、当該座席以外の部分の表面温度または当該座席のうち赤外線カメラ４から着座者によって隠されない部分の表面温度とされている。
【０１５３】
（変形例７）
上記実施形態では、赤外線カメラ４から出力されたサーモグラフィを利用して初期温度Ｔｓｉを特定している。しかし、初期温度Ｔｓｉを特定するために赤外線カメラ４以外のものを利用してもよい。つまり、赤外線カメラ４は必須ではない。
【０１５４】
例えば、外気温センサ１６、内気温センサ１７、日射センサ１８のうち１つまたは複数の検出結果を利用して初期温度Ｔｓｉを特定してもよい。このようにすることで、車室内の空調制御に用いられることが多いセンサを熱流束の推定に流用することで、部品点数を低減させることができる。
【０１５５】
具体的には、エアコンＥＣＵ１４は、図９のステップＳ１３５において、現時点において内気温センサ１７が検出した温度を初期温度Ｔｓｉとして特定してもよい。あるいは、エアコンＥＣＵ１４は、図９のステップＳ１３５において、現時点において外気温センサ１６、内気温センサ１７、日射センサ１８が検出した温度を、所定の温度対応マップに入力することで得られた温度を、初期温度Ｔｓｉとして特定してもよい。ここで、温度対応マップは、対象の座席への着座の際における外気温センサ１６、内気温センサ１７、日射センサ１８の検出値に対する対象の座席の表面温度の対応関係が表されたデータ、エアコンＥＣＵ１４のＲＯＭに記録されている。温度対応マップに表される応関係は、あらかじめ実験等によって決められる。
【０１５６】
また例えば、車両１のすべての座席について、背もたれ部の表面近くおよび座部の表面
近くに温度センサが設置されている場合、これら温度センサの検出結果を利用して初期温
度Ｔｓｉを特定してもよい。
【０１５７】
具体的には、エアコンＥＣＵ１４は、図３のステップＳ１１０において、対象の座席の対象の部分の温度センサが検出した温度を、表面温度として特定してもよい。そしてエアコンＥＣＵ１４は、ステップＳ１３０において、最後にステップＳ１１０で特定した表面温度を、初期温度Ｔｓｉとして特定してもよい。
【０１５８】
また、エアコンＥＣＵ１４は、図９のステップＳ１５０において、現時点において（すなわち着座したと判定された直後に）対象の座席の対象の部分の温度センサが検出した温度を、初期温度Ｔｓｉとして特定してもよい。
【０１５９】
（変形例８）
上記実施形態では、Ｒｉｎｉｔおよびωを算出するための定数ａ１、ａ２、ｋ１、ｋ２の値のそれぞれは、対象の座席にも対象の部分にもよらず同じであった。しかし、定数ａ１、ａ２、ｋ１、ｋ２値のそれぞれは、対象の座席および対象の部分に応じて異なっていてよい。
【０１６０】
また、Ｒｉｎｉｔおよびωをより現実に近く算出するためには、エアコンＥＣＵ１４は、着座者の衣服の材質および厚みに応じて、定数ａ１、ａ２、ｋ１、ｋ２値を変化させてもよい。このためには、エアコンＥＣＵ１４は、着座者の衣服の材質および厚みの情報を、乗員の入力操作によって取得してもよい。そして、着座者の衣服の材質および厚みと上記定数との関係は、あらかじめ実験等によって得られて対応関係が既述されてＲＯＭに記録された対応マップに基づいて決定してもよい。
【０１６１】
また、Ｒｉｎｉｔおよびωをより現実に近く算出するためには、エアコンＥＣＵ１４は、着座者の背中と背もたれ部との間の圧力および着座者の臀部と座部との間の圧力（すなわち、座圧）に応じて、定数ａ１、ａ２、ｋ１、ｋ２の値を変化させてもよい。このためには、エアコンＥＣＵ１４は、上記圧力を、背もたれ部および背面部に設けられた圧力センサの検出結果に基づいて特定してもよい。そして、圧力と上記定数との関係は、あらかじめ実験等によって得られて対応関係が既述されてＲＯＭに記録された対応マップに基づいて決定してもよい。
【０１６２】
（変形例９）
上記実施形態では、第１計算式Σ１の右辺において、着座時点以外の時点の当該対象の部分の表面温度から成る項の係数がゼロである。しかし、第１計算式Σ１の右辺において、着座時点以外の時点の当該対象の部分の表面温度から成る項の係数がゼロでなくてもよい。ただしその場合でも、第１計算式Σ１に対する初期温度依存要素の寄与度は、第１計算式Σ１に対する着座時点以外の時点の当該対象の部分の表面温度から成る要素（例えば項）の寄与度よりも高い。
【０１６３】
（変形例１０）
上記実施形態では、ヒータクーラ５、６、７、８は、取り付け先の座席に着座する人を、あるときには暖め、またあるときには冷やすようになっている。しかし、ヒータクーラ５、６、７、８のそれぞれは、取り付け先の座席に着座する人を、あるときには暖め、またあるときには暖めない電気ヒータ（例えば電熱線）に置き換えられてもよい。また、ヒータクーラ５、６、７、８のそれぞれは、取り付け先の座席に着座する人にあるときには送風し、またあるときには送風しない電動ファンに置き換えられてもよい。この場合、電気ヒータは置き換え対象のヒータクーラと同じ位置に配置されてもよい。
【０１６４】
電動ファンは、吸い込みタイプの電動ファンである。すなわち、電動ファンは、取り付け先の座席２、３に座る人を基準として当該ファンとは反対側から当該ファン側に風を流す。つまり、電動ファンは、座部に取り付けられた場合は、上から下に風を流し、背もたれ部に取り付けられた場合は、前方から後方に風を流す。電動ファンは、取り付け対象の座席を基準として、当該座席に着座する人とは反対側に取り付けられる。例えば、電動ファンは、座部の下側、または背もたれ部の後ろ側に配置される。
【０１６５】
また・BR>Aヒータクーラ５、６、７、８の各々は、上述の電気ヒータおよび電動ファンを備えたヒータファンに置き換えられてもよい。以下、ヒータクーラ５、６、７、８の各々がヒータファンに置き換えられた場合の例について説明する。
【０１６６】
この例では、取り付け先の座席２、３に着座する人を暖める場合は、電気ヒータに通電され、電動ファンへの通電が停止される。これにより、電気ヒータからの熱伝導および電気ヒータからの輻射により、取り付け先の座席２、３に着座する人が暖められる。
【０１６７】
取り付け先の座席２、３に着座する人に送風する場合は、電動ファンに通電され、電気ヒータへの通電が停止される。電動ファンが通電されて作動すると、取り付け先の座席２、３に着座する人を基準として当該座席とは反対側から当該座席側に空気が流れる。つまり、
取り付け先の座席２、３に座る人に当該座席の反対側から送風されると共に、座席２、３内の空気が電動ファンに吸い込まれる。これにより、取り付け先の座席が長時間日射に晒された後に人が当該座席に着座した場合に、当該座席から当該着座者への熱の移動を抑え、熱による着座者の不快感を低減することができる。このように、電動ファンは、座席から熱を奪うことで取り付け先の座席に座る人を冷やすこともできる。
【０１６８】
このような例において、各ヒータファンの作動状態に応じて決定される温調依存要素εは、上述の通り、対象のヒータクーラの現時点の作動状態に応じて決定される。具体的には、上記実施形態と同様、図６のように決定されてもよい。
【０１６９】
なお、低出力冷却モードでは、電動ファンに通電され、電気ヒータに通電されない。また、高出力冷却モードでは、電動ファンに低出力冷却モードよりも高い電力で通電され、電気ヒータに通電されない。また、低出力加熱モードでは、電気ヒータに通電され、電動ファンに通電されない。また、高出力加熱モードでは、電気ヒータに低出力加熱モードよりも高い電力で通電され、電気ヒータに通電されない。
【０１７０】
このようにした場合において、図１０に、運転席２における座部と人の間の熱流束の実測値１０１の推移と、当該熱流束の推移をエアコンＥＣＵ１４が熱流速推定処理１４Ａにおいて推定した推定値１０２と、運転席２の座席近傍温度Ｔｂの実測値１０３の推移について、１つの実験結果を示す。座席近傍温度Ｔｂは、対象の座席（本実験では運転席２）の座部よりも上方かつ対象の座席２の背もたれ部よりも前方であると共に、当該座部からも当該背もたれ部からも３０ｃｍ以上離れた領域における、空気の温度である。当該座部からも当該背もたれ部からも３０ｃｍ以上離れた領域とするのは、着座した人がいない領域とするためである。座席近傍温度Ｔｂは、熱流束推定部位の近辺を流れる空気の温度でもある。
【０１７１】
この実験においては、時点ｔ０において、車室内に送風される前の空気がエバポレータ１１で冷却される冷房運転が始まる。そして、その後の時点ｔ１において運転席２の座部の電動ファンが高出力冷却モードで作動を開始する。そして、その後の時点ｔ２において当該冷房運転が停止する。時点ｔ０以降において、運転席２の座部の電気ヒータは作動していない。
【０１７２】
この実験では、図１０に示すように、熱流束の推定値１０２は、概ね実測値１０１に近い値を算出している。しかし、本発明者は、冷房が停止した時点ｔ２以降において、座席近傍温度Ｔｂの実測値１０３の上昇と共に、熱流束の実測値１０１が上昇しているにもかかわらず、熱流束の推定値１０２が低下していることに着目した。
【０１７３】
そして、本発明者は、座席近傍温度Ｔｂの実測値１０３と熱流束の実測値１０１との相関について検討した結果、電動ファンが吸い込む空気の温度に熱流束が影響を受けることを見出した。
【０１７４】
そこで、本発明者は、各ヒータファンの作動状態に応じて決定される温調依存要素εを、図６ではなく図１１のように決定してもよいことを着想した。図１１と図６の違いは、低出力冷却モードと高出力冷却モードにおけるεの式のみである。
【０１７５】
具体的には、各ヒータファンの作動状態に応じて決定される温調依存要素εは、冷却モードでは、対象のヒータファンの現時点の作動状態および取り付け先の座席の近傍の温度に応じて、図１１に示すように−Ｃ１とｆ（Ｔｂ）の乗算結果または−Ｃ２とｆ（Ｔｂ）の乗算結果となる。Ｃ１、Ｃ２、Ｈ１、Ｈ２は、図６と同じ量である。ｆ（Ｔｂ）は、温度Ｔｂの関数であり、温度Ｔｂによらず正の値を取る。ｆ（Ｔｂ）は、座席近傍温度Ｔｂの増加関数である。例えばｆ（Ｔｂ）はＴｂそのものであってもよい。ｆ（Ｔｂ）が座席近傍温度Ｔｂの増加関数かつ正の値であるのは、当該領域の温度が高いほど、対象の座席の温度低下が遅くなるからである。
【０１７６】
エアコンＥＣＵ１４が座席近傍温度Ｔｂを検出する方法としては、例えば、赤外線カメラ４で撮影された当該領域の画素値によって特定されてもよいし、当該領域にあらかじめ設けられた温度センサによって特定されてもよい。あるいは、内気温センサ１７が検出した温度を座席近傍温度Ｔｂとして用いてもよい。内気温センサ１７は当該領域に配置されてはいないものの、当該領域の温度と強い正の相関があるので、内気温センサ１７の検出値を用いても大きな問題はない。
【０１７７】
このように、図１１のように温調依存要素εを決定することで、電動ファンが吸い込む空気の温度に応じた熱流束の変化も、熱流束の推定に組み込むことができる。
【０１７８】
運転席背部ヒータクーラ５、運転席座部ヒータクーラ６、後席背部ヒータクーラ７、後席座部ヒータクーラ８の各々は、取り付け先の座席２、３に着座する人を、あるときには暖め、またあるときには冷やす温度調整装置である。
【０１７９】
例えば、ヒータクーラ５、６、７、８の各々は、ペルチェ素子から構成されていてもよい。取り付け先の座席２、３に着座する人を暖める場合は、ペルチェ素子の両面のうち、当該人側の面が放熱面となるよう、ペルチェ素子に電流が印加される。取り付け先の座席２、３に着座する人を冷やす場合は、ペルチェ素子の両面のうち、当該人側の面が吸熱面となるよう、ペルチェ素子に電流が印加される。
【０１８０】
（まとめ）
上記各実施形態の一部または全部で示された第１の観点によれば、熱流算出部は、着座時点を起点とする経過時間に依存して変化する時間依存要素と初期温度に依存する初期温度依存要素とを含む第１計算式により、着座時点より後の期間における前記熱流束の推移を算出する。
【０１８１】
また、第２の観点によれば、前記初期温度依存要素の前記第１計算式に対する寄与度は、前記着座時点以外の時点における前記座席の温度の推定値を含む要素の、前記第１計算式に対する寄与度よりも、高い。このように、当該対象の部分の着座後の温度の寄与が小さいので、着座後の測定による誤差が発生する可能性が低い。したがって、熱流束の推定が安定化する。
【０１８２】
また、第３の観点によれば、前記熱流算出部は、前記着座時点における前記座席と前記人との間の熱流束である着座時熱流束を前記初期温度に基づいて推定する。前記時間依存要素は、前記経過時間が長くなるほど前記熱流束が前記着座時熱流束からゼロ側に大きく離れるよう、前記第１計算式に寄与する。
【０１８３】
このようになっていることで、第１計算式によって推定された熱流束の、経過時間の
変化に応じた推移が、実際の熱流束に近くなる。
【０１８４】
また、第４の観点によれば、前記時間依存要素は、前記着座時点を起点とする経過時間が長くなるほど前記熱流束の変化速度の絶対値が小さくなるよう、前記第１計算式に寄与する。このようになっていることで、第１計算式によって推定された熱流束の、経過時間の変化に応じた推移が、実際の熱流束に近くなる。
【０１８５】
また、第５の観点によれば、前記初期温度依存要素は、前記初期温度と基準温度（ａ１）との差の絶対値が大きいほど、前記第１計算式に対する前記時間依存要素の寄与を増幅させる。
【０１８６】
このようになっていることで、第１計算式によって推定された熱流束の、経過時間の変化に応じた変化速度が、初期温度に応じて異なることで、算出された熱流束が実際の熱流束に近くなる。
【０１８７】
また、第６の観点によれば、前記車両の前記座席には、温度調整装置が配置されており、前記温度調整装置は、前記座席に着座する人を暖める、前記座席に着座する人を冷やす、前記座席に着座する人に送風する、のうち少なくとも１つを行い、前記第１計算式は、さらに前記温度調整装置の作動状態に依存する温調依存要素を含む。このように、ヒータクーラの作動状態に依存する温調依存要素を第１計算式が含むことで、温度調整装置の作動状態に起因する熱流束の変化が第１計算式に反映される。
【０１８８】
また、第７の観点によれば、前記熱流算出部は、前記第１計算式により前記熱流束の推移を算出すると共に、前記第１計算式とは異なる第２計算式により前記熱流束の推移を算出し、前記第２計算式は、前記座席の温度を含まず、時間経過と共に前記熱流束の経時変化量の絶対値を減少させ、前記熱流算出部は、前記着座時点の後の切替タイミングよりも前においては、前記第１計算式で算出された前記熱流束の推移を採用し、前記切替タイミングの後においては、前記第２計算式で算出された前記熱流束の推移を採用し、前記切替タイミングは、前記第１計算式で算出された前記熱流束と前記第２計算式で算出された前記熱流束との差の絶対値が所定の許容値より小さくなるタイミングである。
【０１８９】
このようになっていることで、第１計算式のみで熱流束を算出する場合や、第２計算式のみで熱流束を算出する場合に比べて、精度の高い熱流束算出を行うことができる。すなわち、前期は第１計算式によって初期温度に応じた熱流束を算出でき、後期は座席の温度にかかわらず時間経過と共に経時変化量の絶対値が低下していく熱流束の特性を第２計算式で再現することができる。
【０１９０】
また、第８の観点によれば、前記センサは、前記座席を前記座席の前側から撮影する赤外線カメラである。通常、車両に搭載される赤外線カメラは、座席に座った乗員の前側の表面の温度を特定するために用いられる。しかし、ここでは、乗員の背中または臀部と座席との間の、すなわち隠れた部分の、熱流束を推定するために用いられる。つまり、このような赤外線カメラの用途は、従来の用途とは観点が大きく異なっている。
【符号の説明】
【０１９１】
１                車両
２                運転席
３                後席
４                赤外線カメラ
５、６、７、８    ヒータクーラ
１４              エアコンＥＣＵ
【図１】
2020063029.tif 000003
【図２】
2020063029.tif 000004
【図３】
2020063029.tif 000005
【図４】
2020063029.tif 000006
【図５】
2020063029.tif 000007
【図６】
2020063029.tif 000008
【図７】
2020063029.tif 000009
【図８】
2020063029.tif 000010
【図９】
2020063029.tif 000011
【図１０】
2020063029.tif 000012
【図１１】
2020063029.tif 000013