(57)【要約】【課題】本発明は、多段に配置された熱電モジュール全体としての総発電電力の取得量を増大させることができる、熱発電装置の提供を目的とする。【解決手段】熱エネルギーを持った排ガスの流れに沿って配置される２つの熱電モジュール４と、熱電モジュール４に流れる発電電流をそれぞれ調整する電力変換器１２を備え、排ガスの流れの上流側に配置された前段熱電モジュール４Ａに流れる発電電流の目標値をその２つの熱電モジュール４のそれぞれが発電した発電電力を合計した総発電電力が最大となるように設定し、前段熱電モジュール４Ａより下流側に配置された後段熱電モジュール４Ｂに流れる発電電流の目標値を後段熱電モジュール４Ｂ自体の発電電力が最大となるように設定する、ことを特徴とする熱発電装置。【選択図】図５
2008167626.tif 000002
【特許請求の範囲】
【請求項１】
熱エネルギーを持った流体の流れに沿って配置される複数の熱電モジュールと、
前記熱電モジュールに流れる発電電流をそれぞれ制御する発電電流制御手段とを備え、
前記発電電流制御手段は、
流体の流れの上流側に配置された前段熱電モジュールに流れる発電電流の目標値を前記複数の熱電モジュールのそれぞれが発電した発電電力を合計した総発電電力が最大となるように設定し、
前段熱電モジュールより下流側に配置された後段熱電モジュールに流れる発電電流の目標値を後段熱電モジュール自体の発電電力が最大となるように設定する、ことを特徴とする熱発電装置。
【請求項２】
前記発電電流手段は、前記熱電モジュールに流れる発電電流の流通に伴って当該発電電流を減少させる、請求項１に記載の熱発電装置。
【請求項３】
前記発電電流制御手段は、前段熱電モジュールの発電電流の流通／非流通を所定周期で繰り返す、請求項１又は２に記載の熱発電装置。
【請求項４】
前記所定周期のうち流通期間の比率は、前段熱電モジュールの前記複数の熱電モジュール全体としての発電電力を最大にする発電電流を前段熱電モジュールのそれ自体の発電電力を最大にする発電電流で除算した値であるとした、請求項３に記載の熱発電装置。
【請求項５】
前記発電電流手段は、前段熱電モジュールの発電電流の流通期間における平均電流値を前段熱電モジュールのそれ自体の発電電力を最大にする発電電流値に設定し、前段熱電モジュールの前記所定周期期間における平均電流値が前段熱電モジュールの前記複数の熱電モジュール全体としての発電電力を最大にする発電電流値となるように制御する、請求項３又は４に記載の熱発電装置。
【請求項６】
前記発電電流制御手段は、流体の熱容量に応じて、前段熱電モジュールに流れる発電電流の目標値を変化させる、請求項１から５のいずれかに記載の熱発電装置。
【請求項７】
前記発電電流制御手段は、流体の熱容量が少なくなるにつれて、前段熱電モジュールに流れる発電電流の目標値を前段熱電モジュール自体の発電電力を最大にする電流値よりも前段熱電モジュール自体の発電効率を最大にする電流値に近づくように変化させる、請求項６に記載の熱発電装置。
【請求項８】
前記発電電流制御手段は、流体の熱容量が多くなるにつれて、前段熱電モジュールに流れる発電電流の目標値を前段熱電モジュール自体の発電効率を最大にする電流値よりも前段熱電モジュール自体の発電電力を最大にする電流値に近づくように変化させる、請求項６に記載の熱発電装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、熱電効果を利用して発電を行う熱発電装置に関する。
【背景技術】
【０００２】
従来、熱電素子材料で構成された発電モジュールに温度差を生じさせて発電を行う自動車用排熱発電装置が知られている（例えば、特許文献１参照）。この排熱発電装置は、エンジンからの排気を導出する排気管の外表面に一側面が接触するように取付けられる発電モジュールと、当該発電モジュールの他側面に接触するように取付けられる冷却部とによって構成されている。そして、この発電モジュールにおいて、排気管と接触する高温端と冷却部と接触する低温端との間に温度差を生じさせることによって、排温エネルギーを電気エネルギーとして回収している。
【０００３】
しかしながら、エンジンの出力変化などによって排ガスの流量や温度は大きく変動するため、構造的な改良を行って熱電モジュールに温度差を与えて発電する上述の従来技術のままでは、そのような出力変化の対応には限界があり、その出力変化の態様によっては熱の回収効率が低下することが考えられ得る。
【０００４】
そこで、熱電モジュールに流れる電流を積極的に調整制御することによって、熱エネルギーの回収効率を向上させる技術の開発が進んでいる。例えば、熱電モジュールにはその発電電力が最大となる最適電流が存在するため、その最適電流になるように熱電モジュールに流す電流を熱電モジュール毎に調整制御する技術がある。このような調整制御を行うことによって、複数の熱電モジュール全体として発電できる総発電電力の最大値が得られるものと考えられ得る。
【特許文献１】特開２００１−１２２４０号公報
【発明の開示】
【発明が解決しようとする課題】
【０００５】
しかしながら、実際には熱電モジュールに電流を流すと、ペルチェ効果で熱の流れが発生することで熱電素子の熱抵抗が変化するため、発電電力が最大となる電流値と熱電素子の熱抵抗が最適となる電流値は一致しなくなる。すなわち、熱電モジュールに電流を流したことによるペルチェ効果によって熱電素子の熱抵抗が変化して熱電素子の温度差が小さくなるので、熱電モジュールに同じ電流を流したままだと実際に得られる発電電力は少なくなる。したがって、各熱電モジュールに流れる電流を適切に制御しなければ、実際に得られた熱電モジュール全体としての総発電電力が想定より少ない場合がある。
【０００６】
また、熱エネルギーを持った排ガス等の流体に沿って複数の熱電モジュールを多段に配置している場合、流体の上流側に配置された前段熱電モジュールで熱エネルギーが回収されることによって流体の下流側になるほど回収可能な熱エネルギーは減る。そのため、流体の下流側に配置された後段熱電モジュールでの高温端と低温端との温度差は必然的に小さくなることによって、後段熱電モジュール単体での発電電力は減少するとともにその発電効率も低下することになる。したがって、各熱電モジュールに流れる電流を適切に制御しなければ、実際に得られた熱電モジュール全体としての総発電電力が想定より少ない場合がある。
【０００７】
そこで、本発明は、多段に配置された熱電モジュール全体としての総発電電力の取得量を増大させることができる、熱発電装置の提供を目的とする。
【課題を解決するための手段】
【０００８】
上記目的を達成するため、本発明の熱発電装置は、
熱エネルギーを持った流体の流れに沿って配置される複数の熱電モジュールと、
前記熱電モジュールに流れる発電電流をそれぞれ制御する発電電流制御手段とを備え、
御手段とを備え、
前記発電電流制御手段は、
流体の流れの上流側に配置された前段熱電モジュールに流れる発電電流の目標値を前記複数の熱電モジュールのそれぞれが発電した発電電力を合計した総発電電力が最大となるように設定し、
前段熱電モジュールより下流側に配置された後段熱電モジュールに流れる発電電流の目標値を後段熱電モジュール自体の発電電力が最大となるように設定する、ことを特徴としている。これにより、前段熱電モジュール単体の発電電力を最大にするのではなく前後段合わせた熱電モジュール全体の総発電電力を最大にするように設定された発電電流が前段熱電モジュールに流れるので、前段熱電モジュールに吸収されずに後段熱電モジュールに到達する熱エネルギーを増やしつつ、前後段合わせた熱電モジュール全体で得られる総発電電力を増大させることができる。そして、後段熱電モジュール単体の発電電力を最大にするように設定された発電電流が後段熱電モジュールに流れるので、後段熱電モジュール単体で得られる発電電力も増大させることができる。すなわち、多段に配置された熱電モジュール全体としての総発電電力の取得量を増大させることができる。
【０００９】
また、前記発電電流手段は、前記熱電モジュールに流れる発電電流の流通に伴って当該発電電流を減少させることが好ましい。このように発電電流の一定のままではなく発電電流を減少させることによって、ペルチェ効果で熱電素子の温度差が小さくなることによって実際に得られる発電電力が減ることを抑えることができる。
【００１０】
また、前記発電電流制御手段は、前段熱電モジュールの発電電流の流通／非流通を所定周期で繰り返すようにしてもよい。このように流通／非流通を繰り返せば、常時流通する場合よりも、前段熱電モジュールに吸収されずに後段熱電モジュールに到達する熱エネルギーを増やすことができる。なお、前記所定周期のうち流通期間の比率は、例えば、前段熱電モジュールの前記複数の熱電モジュール全体としての発電電力を最大にする発電電流を前段熱電モジュールのそれ自体の発電電力を最大にする発電電流で除算した値であるとする。また、前記発電電流手段は、前段熱電モジュールの発電電流の流通期間における平均電流値を前段熱電モジュールのそれ自体の発電電力を最大にする発電電流値に設定し、前段熱電モジュールの前記所定周期期間における平均電流値が前段熱電モジュールの前記複数の熱電モジュール全体としての発電電力を最大にする発電電流値となるように制御してもよい。
【００１１】
一方、前記発電電流制御手段は、流体の熱容量に応じて、前段熱電モジュールに流れる発電電流の目標値を変化させると好適である。つまり、熱電モジュールの能力に対して、流量と流体の比熱の積が大きいか小さいかによって変化させる。例えば、前記発電電流制御手段は、流体の熱容量が少なくなるにつれて、前段熱電モジュールに流れる発電電流の目標値を前段熱電モジュール自体の発電電力を最大にする電流値よりも前段熱電モジュール自体の発電効率を最大にする電流値に近づくように変化させてもよいし、流体の熱容量が多くなるにつれて、前段熱電モジュールに流れる発電電流の目標値を前段熱電モジュール自体の発電効率を最大にする電流値よりも前段熱電モジュール自体の発電電力を最大にする電流値に近づくように変化させてもよい。このように変化させることにより、前段熱電モジュールに流れる発電電流の目標値を前後段合わせた熱電モジュール全体の総発電電力を最大にする電流値に近づけることができる。
【発明の効果】
【００１２】
本発明によれば、多段に配置された熱電モジュール全体としての総発電電力の取得量を増大させることができる。
【発明を実施するための最良の形態】
【００１３】
以下、図面を参照しながら、本発明を実施するための最良の形態の説明を行う。図１は、本実施形態の熱発電装置に用いる複数の熱電モジュールの配置関係を示した図である。熱電モジュール４は、その一側面がエンジン（図示せず）からの排ガスを導出する排気管２の外側に接触して取付けられ、その一側面と異なる他側面が冷却水を流す冷却装置５の外側に接触して取付けられている。図１の場合、排気管２の外側に沿って５個の熱電モジュール４ａ〜４ｅが配置されている。排気管２の所定の基準位置に対してエンジン側に配置された（すなわち、排ガスの上流側に配置された）熱電モジュールを前段熱電モジュールといい、前段熱モジュール以外の熱電モジュールであって排気管２の所定の基準位置に対して外気への排気口側に配置された（すなわち、排ガスの下流側に配置された）熱電モジュールを後段熱電モジュールという。例えば、熱電モジュール４ａ，４ｂ，４ｃを前段熱電モジュールといい、熱電モジュール４ｄ，４ｅを後段熱電モジュールという。詳細は後述するが、前段熱電モジュールと後段熱電モジュールでは、熱電モジュールに流れる出力電流の調整制御の方法が異なる。
【００１４】
熱電モジュール４は、熱電半導体等の熱電素子を構成する。熱電素子は、熱エネルギーから電気エネルギーへの変換効率ができるだけ高いほうがよい。また、単位変化温度当たりに発生する起電力を表すゼーベック係数Ｓは、できるだけ大きいほうがよい。さらに、電流が流れたときに電気抵抗が大きいと、ジュール熱によってエネルギーが失われるので、熱電素子の電気抵抗は、できるだけ小さいほうがよい。
【００１５】
図２は、熱電モジュール４に生ずる温度差△Ｔを説明するための図である。熱電モジュール４の排気管２との接触面（高温端）と熱電モジュール４の冷却装置５との接触面（低温端）の間に温度差△Ｔが生ずることにより熱電モジュール４に電圧が発生することを利用して（ゼーベック効果）、排気管２を流れる排ガスの熱エネルギーを電気エネルギーに変換して回収することができる。そして、排気管２に沿って配置された熱電モジュール４ａ〜４ｅのそれぞれが、排気管２を流れる排ガスの熱エネルギーを電気エネルギーに変換して回収する。
【００１６】
図３は、本実施形態の熱発電装置の基本構成を示した図である。コンピュータ１４は、排ガス流量などの情報に基づいて熱電モジュール４全体としての発電効率を最適化する制御プログラムやその制御データを記憶するＲＯＭ、制御プログラムの処理データなどを一時的に記憶するＲＡＭ、制御プログラムを処理するＣＰＵ、外部と情報をやり取りするための入出力インターフェースなど、複数の回路要素によって構成されたものである。また、コンピュータ１４は一つの電子制御ユニットとは限らず、制御が分担されるように複数の電子制御ユニットで構成されてもよい。
【００１７】
コンピュータ１４は、熱電モジュール４の熱特性に関する情報（例えば、ゼーベック係数、電気抵抗、性能指数など）を外部から取得したり、予め有していたりする。コンピュータ１４は、エンジンを制御するエンジンコンピュータなどの他のコンピュータや計測センサから、排気管２を流れる排ガスの流量や温度などの情報や、冷却装置５の冷却水の温度などの情報を取得する。コンピュータ１４は、それらの取得情報に基づいて前段熱電モジュールの出力電流Ｉｆ及び後段熱電モジュールの出力電流Ｉｒの最適値などを演算し、その演算結果になるように電力変換器１２に対して指令する。
【００１８】
電力変換器１２は、シリーズレギュレータやスイッチングレギュレータなどの内部回路を備える、入出力間の電力を変換する変換手段である。電力変換器１２の入力側に熱電モジュール４の出力電流Ｉｆ，Ｉｒが入力され、電力変換器１２の出力側からその内部回路によって電力変換された電流Ｉｏｕｔが出力される。したがって、電力変換器１２がコンピュータ１４の指令信号に従って熱電モジュール４の出力電流Ｉｆ，Ｉｒを調整することによって、電力変換器１２の出力側から所望の電力を取り出すことができる。例えば、電力変換器１２から取り出せる電力が最大となる出力電流Ｉｆ，Ｉｒなどの演算値がコンピュータ１４によって指令されるによって、電力変換器１２の出力側から最大電力を取り出すことができる。
【００１９】
ところが、熱電モジュール４に流す電流を変えることによって、熱電モジュール４から取り出せる電力が変わるとともに、熱電モジュール４の発電効率も変わってしまう。発電効率とは、熱電モジュール４において、吸収された熱エネルギーに対して発電できた電力の比率である。図１１は、熱電モジュール４の出力電流Ｉと発電電力Ｐ又は発電効率ηとの関係を示した図である。図１１に示されるように、それらの関係は或る出力電流Ｉで極大点を有する放物線のような特性を有している。また、熱電モジュール４の発電電力を最大にする出力電流と熱電モジュール４の発電効率を最大にする出力電流は異なり、最大発電効率になる出力電流は最大発電電力になる出力電流より小さい。
【００２０】
電力を取り出すために熱電モジュール４に電流を流すと、ペルチェ効果で熱の流れが発生することで熱電モジュール４内の熱電素子の熱抵抗が変化するため、電力が最大になる出力電流値と熱電素子の熱抵抗が最適となる電流値が一致せず、発電効率は低下してしまる。すなわち、図１１に示されるように、熱電モジュール４の発電電力が最大となる出力電流Ｉｏに設定したとしても（点Ａ１）、ペルチェ効果により熱が吸収されることにより熱電素子の温度差が小さくなるため、出力電流Ｉｏのままでは発電電力は減少してしまう（点Ｃ１）とともに発電効率も低下してしまう（点Ｄ１）。したがって、各熱電モジュールに流れる電流を適切に制御しなければ、実際に得られた熱電モジュール全体としての総発電電力が想定より少ない場合がある。
【００２１】
また、エンジンの出力が低い場合など排ガスの熱エネルギーが小さい場合には、排ガスが通る排気管２に沿って配置された熱電モジュール４を通過する熱量が小さくなるため、熱電モジュール４に生ずる温度差が不十分となりその発電量は低下する。とりわけ、複数の熱電モジュール４を排気管２に沿って多段に並べている場合、排気管２の上流側に配置された前段熱電モジュールで熱エネルギーが回収されることによって排気管２の下流側になるほど回収可能な熱エネルギーは減る。そのため、排気管２の下流側に配置された後段熱電モジュールでの高温端と低温端との温度差は必然的に小さくなることによって、後段熱電モジュール単体での発電電力は減少するとともにその発電効率も低下することになる。したがって、各熱電モジュールに流れる電流を適切に制御しなければ、実際に得られた熱電モジュール全体としての総発電電力が想定より少ない場合がある。
【００２２】
図４は、排ガスの熱エネルギーと回収エネルギーとの関係を示した模式図である。図４（ａ）のように、エンジンの出力が大きいこと等により集熱フィン３を介して得られる排ガスの熱量が大きい場合には、排気管２に沿って多段に配置された各熱電モジュール４の発電量が最大となるようにそれらの出力電流をそれぞれ設定したとしても、下流側の熱電モジュール（例えば、４ｅ）からでも十分な発電量を得ることができる。しかしながら、図４（ｂ）のように、エンジンの出力が小さいこと等により集熱フィン３を介して得られる排ガスの熱量が小さい場合には、各熱電モジュール４の発電量が最大となるようにそれらの出力電流をそれぞれ設定すると、上流側の前段熱電モジュール（例えば４ａ）単独で熱がほとんど奪われ、その下流側に配置された後段熱電モジュールの熱の取り分が少なくなる。その結果、発電量が最大となるように出力電流を設定したとしてもそのときの発電効率は必ずしも最適ではないため、発電モジュール全体としての発電効率は低下し、取り出せるべき総発電量も減っていることになる。
【００２３】
そこで、コンピュータ１４は、熱電モジュール全体としての総発電電力の取得量を増大させるために、以下の各実施例に記載の理論計算に基づいて導出された出力電流Ｉｆ，Ｉｒ等の算出式に従って、熱電モジュール４の出力電流Ｉｆ，Ｉｒを電力変換器１５によって制御する。なお、計算の簡単化のため、図５に示されるように、熱電モジュール４が前段と後段の２段から構成されている場合について考える。
【実施例１】
【００２４】
実施例１では、熱電モジュール全体としての総発電電力の取得量を増大させるために、コンピュータ１４は、前段熱電モジュール４Ａと後段熱電モジュール４Ｂを合わせた総発電電力が最大となるように前段熱電モジュール４Ａの出力電流Ｉｆをデューティ比Ｄでオン／オフ制御し、後段熱電モジュール４Ｂ自体の発電電力が最大となるように後段熱電モジュール４Ｂの出力電流Ｉｒを制御する。さらに、コンピュータ１４は、エンジンの出力が小さいこと等により集熱フィン３Ａを介して得られる排ガスの熱量や温度が低下するにつれて、前段熱電モジュール４Ａの出力電流Ｉｆをオフする時間を長くするように制御する。このように制御することで、前段熱電モジュール４Ａで回収される排ガスの熱量を減らすとともに後段熱電モジュール４Ｂで回収される排ガスの熱量を増やすことができるので、前段熱電モジュール４Ａと後段熱電モジュール４Ｂを合わせた全体としての最大の発電電力を取り出すことができる。
【００２５】
図６は、実施例１における、前段熱電モジュール４Ａの出力電流が間欠に流れているときの電流波形である。前段熱電モジュール４Ａの出力電圧Ｖｆと後段熱電モジュール４Ｂの出力電圧Ｖｒは、熱電モジュールに電流が流れることによる電圧降下を考慮して、式（１），（２）によって表すことができる。なお、Ｓは熱電モジュールのゼーベック係数（定数）、Ｒは熱電モジュールの電気抵抗（定数）、ΔＴｆは前段熱電モジュール４Ａにおける温度差、Ｉｆは前段熱電モジュール４Ａの出力電流、ΔＴｒは後段熱電モジュール４Ｂにおける温度差、Ｉｒは後段熱電モジュール４Ｂの出力電流である。また、実施例１内の演算式中のＩｆ，Ｖｆのそれぞれは、図６に示されるように、前段発電モジュール４Ａの出力電流がオンしている状態での（通電状態での）電流値と電圧値として定義する。
【００２６】
【数１】
2008167626.tif 000003
ここで、出力電流Ｉｆが零のときの前段熱電モジュール４Ａの温度差をΔＴｆｏと定義するとともに、ペルチェ効果で熱電モジュールに生ずる温度差が減少する係数をαと定義するならば、前段熱電モジュール４Ａに発生する温度差ΔＴｆは、式（３）によってモデル化することができる。また、前段の出力電流Ｉｆ及び後段の出力電流Ｉｒが零のときの後段熱電モジュールの温度差をΔＴｒｏと定義するとともに、前段熱電モジュール４Ａの出力電流Ｉｆによるペルチェ効果で後段熱電モジュール４Ｂに生ずる温度差が減少する係数をβと定義するならば、後段熱電モジュール４Ｂに発生する温度差ΔＴｒは、式（４）によってモデル化することができる。
【００２７】
【数２】
2008167626.tif 000004
ここで、式（４）の第３項『β・Ｉｆ・Ｄ』の意味は、前段の出力電流Ｉｆを流すとペルチェ効果によって排ガス温度が低下して後段熱電モジュール４Ｂに発生する温度差△Ｔｒが小さくなるということを示している。そして、βの設定値を大きくするにつれて後段熱電モジュール４Ｂに発生する温度差△Ｔｒが小さくなりやすいモデルとなり、βの設定値を小さくするにつれて後段熱電モジュール４Ｂに発生する温度差△Ｔｒが大きくなりやすいモデルとなる。したがって、βの設定値を変更することによって後段熱電モジュール４Ｂに発生する温度差△Ｔｒを大きくするか小さくするかの優先度合いを決めることができる。
【００２８】
△Ｔｆｏ，△Ｔｒｏ，α，βは、例えば排ガスの温度と流量によって決まる定数とする。排ガスの温度や流量が上昇すると排ガスの熱容量が増えるので、βを排ガスの温度や流量が上昇するにつれて小さくなるような可変値に設定することによって、排ガスの温度が高かったり流量が多かったりするときには前段の出力電流Ｉｆを流しても後段熱電モジュール４Ｂに発生する温度差△Ｔｒは小さくなりにくいという現象をモデル式（４）に反映させることができる。αについても同様の考えの下、αを排ガスの温度や流量が上昇するにつれて小さくなるような可変値に設定することによって、排ガスの温度が高かったり流量が多かったりするときには前段の出力電流Ｉｆを流しても前段熱電モジュール４Ａに発生する温度差△Ｔｆは小さくなりにくいという現象や後段の出力電流Ｉｒを流しても後段熱電モジュール４Ｂに発生する温度差△Ｔｒは小さくなりにくいという現象をモデル式（４）（５）に反映させることができる。
【００２９】
一方、前段熱電モジュール４Ａから取り出せる電力をＰｆとし、後段熱電モジュール４Ｂからの取り出せる電力をＰｒとすると、前後段合わせた熱電モジュール全体から取り出せる総電力Ｐは、式（５）で表せる。式（５）に式（１）〜（４）を代入すると、総電力Ｐは式（６）のように変形することができる。
【００３０】
【数３】
2008167626.tif 000005
電力変換器１２Ａ，Ｂは、出力電流Ｉｆ，Ｉｒを調整することによって、総電力Ｐが最大となるように制御する。総電力Ｐを最大にする前段熱電モジュール４Ａの出力電流Ｉｆｍａｘは、総電力ＰをＩｆについて偏微分した偏微分方程式（７）を満たすＩｆに相当する。したがって、総電力Ｐを最大にする前段熱電モジュール４Ａの出力電流Ｉｆｍａｘは、式（８）と導出することができる。
【００３１】
【数４】
2008167626.tif 000006
また、総電力Ｐを最大にするデューティ比Ｄｍａｘは、総電力Ｐを最大にする前段熱電モジュール４Ａの出力電流Ｉｆｍａｘと同様に、総電力ＰをＤについて偏微分した偏微分方程式（９）を満たすＤに相当する。したがって、総電力Ｐを最大にするデューティ比Ｄｍａｘは、式（１０）と導出することができる。
【００３２】
【数５】
2008167626.tif 000007
また、後段熱電モジュール４Ｂ自体の発電電力が最大となる後段熱電モジュール４Ｂの出力電流Ｉｒｍａｘは、図１１に示されるように、検出された排ガスの温度や流量のときの後段熱電モジュール４Ｂの発電電力を最大にする出力電流に設定される。後段熱電モジュール４Ｂ自体の発電電力が最大となる後段熱電モジュール４Ｂの出力電流Ｉｒｍａｘについては、排ガスの温度や流量及び冷却水の温度などに応じて変化する。そこで、予めシミュレーション等によって、排ガスの温度や流量及び冷却水の温度と後段熱電モジュール４Ｂ自体の発電電力が最大となる後段熱電モジュール４Ｂの出力電流Ｉｒｍａｘとの関係をメモリに記憶させておく。計測センサ等によって実際に計測された排ガスの温度や流量及び冷却水の温度に基づいて、そのメモリに記憶された対応関係を用いて、後段熱電モジュール４Ｂの発電電力を最大にする出力電流を設定する。
【００３３】
したがって、式（８）と式（１０）の連立解が求められることによって、コンピュータ１４（電力変換器１２Ａ）が制御すべき前段熱電モジュール４Ａの出力電流Ｉｆと前段熱電モジュール４Ａのデューティ比Ｄが決定される。すなわち、コンピュータ１４は、式（８）と式（１０）に従って、前段熱電モジュール４Ａの出力電流Ｉｆと前段熱電モジュール４Ａのデューティ比Ｄを制御するとともに、上述のＩｒｍａｘに従って後段熱電モジュール４Ｂの出力電流Ｉｒを制御する。ただし、デューティ比Ｄｍａｘの計算値が１を超える場合にはデューティ比Ｄを１に固定する。また、式（８）と式（１０）内の未知変数については、シミュレーション等によって適当な値を予め算出しておく。
【００３４】
そして、エンジンの出力が小さいこと等により排ガスの温度や流量が低下するにつれてα，βは大きくなるように連続的に変化するので、コンピュータ１４によって制御される前段熱電モジュール４Ａの出力電流Ｉｆとそのデューティ比Ｄは、式（８）（１０）に従いながら、排ガスの温度や流量が低下するにつれて小さくなるように補正される。その結果、前段熱電モジュール４Ａにおける熱エネルギーの取り分を減らして、後段熱電モジュール４Ｂにおける熱エネルギーの取り分を増やすことができる。
【００３５】
したがって、前段熱電モジュール４Ａに吸収されずに後段熱電モジュール４Ｂに到達する熱エネルギーを増やしつつ、前後段合わせた熱電モジュール全体で得られる総発電電力を増大させることができる。そして、後段熱電モジュール４Ｂ単体の発電電力を最大にするように設定された発電電流が後段熱電モジュール４Ｂに流れるので、後段熱電モジュール４Ｂ単体で得られる発電電力も増大させることができる。すなわち、多段に配置された熱電モジュール全体としての総発電電力の取得量を増大させることができる。
【実施例２】
【００３６】
実施例２では、実施例１と同様に、熱電モジュール全体としての総発電電力の取得量を増大させるために、コンピュータ１４は、前段熱電モジュール４Ａと後段熱電モジュール４Ｂを合わせた発電電力が最大となるように前段熱電モジュール４Ａの出力電流Ｉｆをデューティ比Ｄでオン／オフ制御し、後段熱電モジュール４Ｂ自体の発電電力が最大となるように後段熱電モジュール４Ｂの出力電流Ｉｒを制御する。さらに、コンピュータ１４は、エンジンの出力が小さいこと等により集熱フィン３Ａを介して得られる排ガスの熱量や温度が低下するにつれて、前段熱電モジュール４Ａの出力電流Ｉｆをオフする時間を長くするように制御する。実施例２では、実施例１と異なる導出方法で、総電力Ｐを最大にする前段熱電モジュール４Ａの出力電流Ｉｆｍａｘ及びデューティ比Ｄｍａｘを導出する。
【００３７】
図７は、実施例２における、前段熱電モジュール４Ａの出力電流が間欠に流れているときの電流波形である。前段熱電モジュール４Ａの出力電圧Ｖｆと後段熱電モジュール４Ｂの出力電圧Ｖｒは、熱電モジュールに電流が流れることによる電圧降下を考慮して、式（１１），（１２）によって表すことができる。なお、Ｓは熱電モジュールのゼーベック係数（定数）、Ｒは熱電モジュールの電気抵抗（定数）、ΔＴｆは前段熱電モジュール４Ａにおける温度差、Ｉｆは前段熱電モジュール４Ａの出力電流、ΔＴｒは後段熱電モジュール４Ｂにおける温度差、Ｉｒは後段熱伝モジュール４Ｂの出力電流である。
【００３８】
【数６】
2008167626.tif 000008
ここで、出力電流Ｉｆが零のときの前段熱電モジュール４Ａの温度差をΔＴｆｏと定義するとともに、ペルチェ効果で熱電モジュールに生ずる温度差が減少する係数をαと定義するならば、前段熱電モジュール４Ａに発生する温度差ΔＴｆは、式（１３）によってモデル化することができる。また、前段の出力電流Ｉｆ及び後段の出力電流Ｉｒが零のときの後段熱電モジュールの温度差をΔＴｒｏと定義するとともに、前段熱電モジュール４Ａの出力電流Ｉｆによるペルチェ効果で後段熱電モジュール４Ｂに生ずる温度差が減少する係数をβと定義するならば、後段熱電モジュール４Ｂに発生する温度差ΔＴｒは、式（１４）によってモデル化することができる。
【００３９】
【数７】
2008167626.tif 000009
ここで、△Ｔｆｏ，△Ｔｒｏ，α，βは、排ガスの温度と流量によって決まる定数とする。実施例１と同様の考えの下、α，βは排ガスの温度や流量が増加するにつれて小さくなるような可変値に設定する。
【００４０】
一方、前段熱電モジュール４Ａから取り出せる電力をＰｆとし、後段熱電モジュール４Ｂからの取り出せる電力をＰｒとすると、前後段合わせた熱電モジュール全体から取り出せる総電力Ｐは、式（１５）で表せる。式（１５）に式（１１）〜（１４）を代入すると、総電力Ｐは式（１６）のように変形することができる。
【００４１】
【数８】
2008167626.tif 000010
電力変換器１２Ａ，Ｂは、出力電流Ｉｆ，Ｉｒを調整することによって、総電力Ｐが最大となるように制御する。総電力Ｐを最大にする前段熱電モジュール４Ａの出力電流Ｉｆｍａｘ０は、総電力ＰのＩｆについて偏微分した偏微分方程式（１７）を満たすＩｆに相当する。したがって、総電力Ｐを最大にする前段熱電モジュール４Ａの出力電流Ｉｆｍａｘ０は、式（１８）として導出することができる。
【００４２】
【数９】
2008167626.tif 000011
また、前段熱電モジュール４Ａから取り出せる電力Ｐｆを最大にする前段熱電モジュールの出力電流Ｉｆｍａｘ１は、式（１８）においてβを零にした場合に相当するので、式（１９）として導出することができる。
【００４３】
【数１０】
2008167626.tif 000012
そこで、前段熱電モジュール４Ａの出力電流のオン期間における平均電流値をＩｆｍａｘ１に設定し、オフ期間も含めた平均電流値がＩｆｍａｘ０となるようにオン／オフ制御を行えば、総電力Ｐが最大となるようにしつつ前段熱電モジュール４Ａの電力Ｐｆが最大となるような最大の発電効率が得られることになる。このときのオン／オフ制御におけるオン期間ｔｏｎと周期ｔｐの比（デューティ比Ｄｍａｘに相当）は、式（２０）として導出することができる。
【００４４】
【数１１】
2008167626.tif 000013
また、後段熱電モジュール４Ｂ自体の発電電力が最大となる後段熱電モジュール４Ｂの出力電流Ｉｒｍａｘは、実施例１と同様に、検出された排ガスの温度や流量のときの後段熱電モジュール４Ｂの発電電力を最大にする出力電流に設定される。
【００４５】
したがって、式（１８）と式（２０）の連立解が求められることによって、コンピュータ１４（電力変換器１２Ａ）が制御すべき前段熱電モジュール４Ａの出力電流Ｉｆと前段熱電モジュール４Ａのデューティ比Ｄが決定される。すなわち、コンピュータ１４は、式（１８）と式（２０）に従って、前段熱電モジュール４Ａの出力電流Ｉｆと前段熱電モジュール４Ａのデューティ比Ｄを制御するとともに、上述のＩｒｍａｘに従って後段熱電モジュール４Ｂの出力電流Ｉｒを制御する。ただし、デューティ比Ｄｍａｘの計算値が１を超える場合にはデューティ比Ｄを１に固定する。また、式（１８）と式（２０）内の未知変数については、シミュレーション等によって適当な値を予め算出しておく。
【００４６】
そして、エンジンの出力が小さいこと等により排ガスの温度や流量が低下するにつれてα，βは大きくなるように連続的に変化するので、コンピュータ１４によって制御される前段熱電モジュール４Ａの出力電流Ｉｆとそのデューティ比Ｄは、式（８）（１０）に従いながら、排ガスの温度や流量が低下するにつれて小さくなるように補正される。その結果、前段熱電モジュール４Ａにおける熱エネルギーの取り分を減らして、後段熱電モジュール４Ｂにおける熱エネルギーの取り分を増やすことができる。
【００４７】
したがって、前段熱電モジュール４Ａに吸収されずに後段熱電モジュール４Ｂに到達する熱エネルギーを増やしつつ、前後段合わせた熱電モジュール全体で得られる総発電電力を増大させることができる。そして、後段熱電モジュール４Ｂ単体の発電電力を最大にするように設定された発電電流が後段熱電モジュール４Ｂに流れるので、後段熱電モジュール４Ｂ単体で得られる発電電力も増大させることができる。すなわち、多段に配置された熱電モジュール全体としての総発電電力の取得量を増大させることができる。
【００４８】
このように、上述の実施例１，２によると、出力電流がオフ期間のときは電力を取り出していないので、熱電モジュール全体としての総発電電力は電流をオンし続けた場合より小さくなるように思われる。しかしながら、実際には、エンジンの出力が小さい場合には排ガスの熱量が小さく、後段の熱電モジュールでは著しく出力電力が小さくなっている。そこで、出力電流を間欠的に制御することによって前段から取り出される熱エネルギーが小さくなるため、後段の熱電モジュール部分での排ガス温度が出力電流のオンを継続した場合よりも高くなり、後段の熱電モジュールで取り出せる電力は向上する。一方、前段の熱電モジュールでは熱電素子の温度差ΔＴを高く保っていることで発電効率が高く保たれる。したがって、前後段合わせた発電モジュール全体としての総発電電力を一層増大させることができる。
【００４９】
なお、上述の実施例１，２において、前段熱電モジュール４Ａの出力電流がオン期間のとき、その電流が流れることによってペルチェ効果が生じ、熱電モジュールにおける温度差ΔＴは減少する。△Ｔの減少につれて最大電力が得られる電流値も減少してしまう。そこで、コンピュータ１４は、図８に示されるように、最大電力が得られるように、前段熱電モジュール４Ａの出力電流のオン期間における電流値を徐々に減少させる補正制御を行う。
【００５０】
図９は、熱電モジュール４の出力電圧と出力電流の関係を示した図である。直線１は温度差ΔＴが大きい場合の電流−電圧特性であり、直線Ｂは温度差ΔＴが直線１の条件より小さい場合の電流−電圧特性である。これらの電流−電圧特性は、熱電モジュール４内の熱電材料の電気抵抗によって決まる。図９のＡ点〜Ｄ点は、図８のＡ点〜Ｄ点にそれぞれ対応する。Ａ点は直線１のＤ点とＥ点の中点であり、Ｂ点は直線２のＣ点とＦ点の中点である。すなわち、Ａ点、Ｂ点は各温度差△Ｔにおいて得られるべき最大電力点である。
【００５１】
出力電流をオンした直後では温度差ΔＴは大きいので電流−電圧特性は直線１上にあり、この状態での最大電力点に電流を制御するとＡ点に到達する。そして、ペルチェ効果によって温度差ΔＴが徐々に減少するので、この間に最大電力が得られるように電流を制御することで、Ａ点からＢ点に徐々に移動させる。次に、温度差ΔＴの減少を抑えるために電流をオフすることでＣ点に到達する。ペルチェ効果が無くなったことで温度差ΔＴが高い値に回復してＤ点に到達する。再び、電流をオンしてＡ点に到達する。以上を繰り返す。以上の繰り返し周期は数秒程度以下の時間周期で行うことが望ましい。
【００５２】
なお、Ａ点からＢ点までの間に最大電力が得られるように電流をＡ点からＢ点に徐々に移動させているが、出力電流を多少変化させてＡ点からＢ点までの間の最大電力点を探索しながらＡ点からＢ点まで移動させると望ましい。図１０は、最大電力点を探索した場合の電流変化を示した図である。図７と図１０に示される出力電流のデューティ比は同じでよい。
【００５３】
このような最大電力点の探索方法は、山登り法などが知られている。図１２は、出力電流Ｉを変化させることによって最大電力点を探索する山登り法を表現した図である。具体的には、現状の動作点で得られる電力と移動させた動作点で得られる電力の増減により動作点への方向判断を行い、この２点間の各動作点での電力増減判定操作を繰り返し行うことにより、最適点へと追従させる方法である。例えば、動作点を△Ｉ移動させた時に電力が増加している場合には、△Ｉ移動させた方向と同じ方向に次回の動作点を移動させ、動作点を△Ｉ移動させた時に電力が減少している場合には、最適点は今回動作させた逆方向にあるとして、その逆方向に次回の動作点を移動させ、この２点間の各動作点での電力増減判定操作を繰り返し行うことにより最適点へと追従させることができる。
【実施例３】
【００５４】
実施例３では、熱電モジュール全体としての総発電電力の取得量を増大させるために、コンピュータ１４は、エンジンの出力が小さくなること等により排ガスの熱エネルギーが少なくなるにつれて（例えば、排ガスの流量が少なくなるにつれて）、前段熱電モジュール４Ａの出力電流が前段熱電モジュール４Ａ自体の発電電力を最大にする電流よりも前段熱電モジュール４Ａ自体の効率を最大にする電流に近づくように制御する。そして、コンピュータ１４は、エンジンの出力が大きくなること等により排ガスの熱エネルギーが多くなるにつれて（例えば、排ガスの流量が多くなるにつれて）、前段熱電モジュール４Ａの出力電流が前段熱電モジュール４Ａ自体の効率を最大にする電流よりも前段熱電モジュール４Ａ自体の発電電力を最大にする電流に近づくように制御する。
【００５５】
実施例２と同様の導出方法によって、総電力Ｐを最大にする前段熱電モジュール４Ａの出力電流Ｉｆｍａｘは、式（１８）として導出することができる。
【００５６】
式（１８）によって表されるＩｆｍａｘを演算するため、式（１８）内のαとβを算出する必要がある。αとβの算出にあたり、まず、出力電流Ｉｆが零のときの前段熱電モジュール４Ａの温度差ΔＴｆｏを算出する。前段熱電モジュール４Ａと排気管２との接触部の排ガス温度をＴｇｆ、冷却装置５の冷却水の温度をＴｗ、熱電モジュール４の熱抵抗をγｂと定義すると、△Ｔｆｏは式（２１）によって表すことができる。
【００５７】
また、前段熱電モジュール４ＡでＱ［Ｗ］の熱が取られるとすると、排ガス流量をＧａ［ｋｇ／ｓｅｃ］と定義するとともに、排ガスの比熱をＣｇ［Ｊ／（ｋｇ℃）］と定義するならば、後段熱電モジュール４Ｂと排気管２との接触部の排ガス温度Ｔｇｒは、式（２２）によって表すことができる。
【００５８】
また、前段熱電モジュール４Ａで取られる熱Ｑは、ペルチェ効果と熱伝導によって、式（２３）のように表すことができる。式（２３）の第１項がペルチェ効果を表し、式（２３）の第２項が熱伝導を表す。ここで、Ｔは前段熱電モジュール４Ａの代表温度とし、ｔｍは前段熱電モジュール４Ａの厚さ、λｍは前段熱電モジュール４Ａの熱伝導率、Ａは前段熱電モジュール４Ａの断面積である。
【００５９】
【数１２】
2008167626.tif 000014
式（２２），（２３）より、後段熱電モジュール４Ｂと排気管２との接触部の排ガス温度Ｔｇｒは式（２４）によって表すことができ、後段熱電モジュール４Ｂに発生する温度差△Ｔｒは式（２５）によって表すことができる。なお、Ｋは式（２６）とする。
【００６０】
【数１３】
2008167626.tif 000015
したがって、熱電モジュールの熱流れの考察より、熱電モジュールの発電性能を示す性能指数をＺＴとすると、αは式（２７）によって表すことができるとともに、式（２５）により、βは式（２８）によって表すことができる。また、前段の出力電流Ｉｆ及び後段の出力電流Ｉｒが零のときの後段熱電モジュールの温度差△Ｔｒｏは式（２９）によって表すことができる。
【００６１】
【数１４】
2008167626.tif 000016
よって、式（１８）に示される、総電力Ｐを最大にする前段熱電モジュール４Ａの出力電流Ｉｆｍａｘは、式（３０）のように算出することができる。
【００６２】
【数１５】
2008167626.tif 000017
さらに、後段熱電モジュール４Ｂの出力電流Ｉｒを最適化する。すなわち、上述と同様に式（１６）で示される総電力ＰをＩｒについての偏微分が零である偏微分方程式を解くことによって、総電力Ｐを最大にする後段熱電モジュール４Ｂの出力電流Ｉｒｍａｘが式（３１）（３２）のように導出することができる。
【００６３】
【数１６】
2008167626.tif 000018
したがって、式（３０）内のＩｒに式（３２）のＩｒｍａｘを代入することによって、最終的な前段熱電モジュールの出力電流の解Ｉｆを求めると、式（３３）のように表すことができる。
【００６４】
【数１７】
2008167626.tif 000019
また、後段熱電モジュール４Ｂ自体の発電電力が最大となる後段熱電モジュール４Ｂの出力電流Ｉｒｍａｘは、実施例１と同様に、検出された排ガスの温度や流量のときの後段熱電モジュール４Ｂの発電電力を最大にする出力電流に設定される。
【００６５】
したがって、コンピュータ１４は、式（３３）に示されるＩｆｍａｘと上述のＩｒｍａｘに従って、前段熱電モジュール４Ａの出力電流Ｉｆと後段熱電モジュール４Ｂの出力電流Ｉｒを制御する。また、式（３３）内の未知変数については、シミュレーション等によって適当な値を予め算出しておく。
【００６６】
したがって、前段熱電モジュール４Ａに吸収されずに後段熱電モジュール４Ｂに到達する熱エネルギーを増やしつつ、前後段合わせた熱電モジュール全体で得られる総発電電力を増大させることができる。そして、後段熱電モジュール４Ｂ単体の発電電力を最大にするように設定された発電電流が後段熱電モジュール４Ｂに流れるので、後段熱電モジュール４Ｂ単体で得られる発電電力も増大させることができる。すなわち、多段に配置された熱電モジュール全体としての総発電電力の取得量を増大させることができる。
【００６７】
そして、式（３３）内の変数Ｇａの位置からわかるように、排ガスの流量Ｇａが小さくなるにつれて電流最適値Ｉｆｍａｘを小さな値に補正することができ、排ガスの流量Ｇａが大きくなるにつれて電流最適値Ｉｆｍａｘを大きな値に補正することができる。
【００６８】
図１３（ａ）は、排ガス流量が少ない場合（Ｇａ＝０．０２［ｋｇ／ｓｅｃ］）における、前段発電モジュール４Ａの出力電流と前段発電モジュール４Ａの又は前後段合わせた発電モジュール全体としての発電電力との関係、並びに前段発電モジュール４Ａの出力電流と前段発電モジュール４Ａ自体の効率との関係を、上述のモデル計算に従って得られた結果を示した図である。図１３（ｂ）は、排ガス流量が多い場合（Ｇａ＝０．０６［ｋｇ／ｓｅｃ］）における、それらの関係を上述のモデル計算に従って得られた結果を示した図である。モデル計算時の各パラメータは、図１３に記載の数値とする。
【００６９】
図１３において、前後段合わせた熱電モジュール全体としての総発電電力を最大にする前段の出力電流は、前段熱電モジュール自身の発電電力を最大にする前段の出力電流より少ないことがわかる。これは、前段で奪う熱エネルギーを少なくする電流値であって後段の発電電力を電力が多く取り出せる電流値に設定されることによって、前後段合わせた熱電モジュール全体としての総発電電力を最大にしていることにほかならない。
【００７０】
また、図１３（ｂ）の排ガス流量Ｇａが多い場合における前後段合わせた熱電モジュール全体としての総発電電力を最大にする電流値は、図１３（ａ）の排ガス流量Ｇａが少ない場合における前後段合わせた熱電モジュール全体としての総発電電力を最大にする電流値に比べ、前段熱電モジュール自身の発電電力を最大にする電流値に近くなっている。これは、排ガスの熱エネルギーに余裕があるため前段で奪う熱エネルギーの後段への影響が小さいためである。
【００７１】
以上、本発明の好ましい実施例について詳説したが、本発明は、上述した実施例に制限されることはなく、本発明の範囲を逸脱することなく、上述した実施例に種々の変形及び置換を加えることができる。
【００７２】
例えば、上述の実施形態は、熱電モジュール４が前段と後段の２段から構成されている場合であるか、熱電モジュール４が３段以上配置された場合でも上記同様の観点で考えればよい。
【図面の簡単な説明】
【００７３】
【図１】本実施形態の熱発電装置に用いる複数の熱電モジュールの配置関係を示した図である。
【図２】熱電モジュール４に生ずる温度差△Ｔを説明するための図である。
【図３】本実施形態の熱発電装置の基本構成を示した図である。
【図４】排ガスの熱エネルギーと回収エネルギーとの関係を示した模式図である。
【図５】熱電モジュール４を２段とする計算モデルを示した図である。
【図６】前段の発電モジュールの出力電流Ｉｆについての電流波形である。
【図７】前段の発電モジュールの出力電流Ｉｆについての電流波形である。
【図８】発電モジュール４の出力電流と温度差との関係を示した図である。
【図９】熱電モジュール４の出力電圧と出力電流の関係を示した図である。
【図１０】最大電力点を探索した場合の電流変化を示した図である。
【図１１】熱電モジュール４の出力電流Ｉと発電電力Ｐ又は発電効率ηとの関係を示した図である。
【図１２】出力電流Ｉを変化させることによって最大電力点を探索する山登り法を表現した図である。
【図１３】前段発電モジュール４Ａの出力電流と前段発電モジュール４Ａの又は前後段合わせた発電モジュール全体としての発電電力との関係、並びに前段発電モジュール４Ａの出力電流と前段発電モジュール４Ａ自体の効率との関係を、本発明に係るモデル計算に従って得られた結果を示した図である。
【符号の説明】
【００７４】
２  排気管
３  集熱フィン
４  熱電モジュール
５  冷却装置
１２  電力変換器
１４  コンピュータ
【図１】
2008167626.tif 000020
【図２】
2008167626.tif 000021
【図３】
2008167626.tif 000022
【図４】
2008167626.tif 000023
【図５】
2008167626.tif 000024
【図６】
2008167626.tif 000025
【図７】
2008167626.tif 000026
【図８】
2008167626.tif 000027
【図９】
2008167626.tif 000028
【図１０】
2008167626.tif 000029
【図１１】
2008167626.tif 000030
【図１２】
2008167626.tif 000031
【図１３】
2008167626.tif 000032