【手続補正書】
【提出日】平成17年8月19日(2005.8.19)
【手続補正１】
【補正対象書類名】明細書
【補正対象項目名】全文
【補正方法】変更
【補正の内容】
【書類名】          明細書
【発明の名称】空力音源探査システム及び空力音源探査方法
【特許請求の範囲】
【請求項１】
流れ場を可視化するためのトレーサ粒子を発生するトレーサ粒子発生手段と、
流れ場をシート状照明光により照明する光シート照明手段と、
前記光シート照明手段により照明された所定の測定領域を撮影するための撮影手段と、
前記撮影手段により撮影された複数の撮影画像を格納する記録手段と、
前記記録手段に記録されている微小時間間隔で撮影された撮影画像から、流れ場の速度ベクトルの分布を解析する第１の演算手段と、
前記第１の演算装置により解析された速度ベクトルの分布から、空力音の二重極音の音源分布を解析する第２の演算手段と、
を備えたことを特徴とする空力音源探査システム。
【請求項２】
前記第２の演算手段は、Howeの渦音の式にグリーン関数を適用した式を用いて音源分布を解析することを特徴とする請求項１記載の空力音源探査システム。
【請求項３】
前記グリーン関数は、コンパクトグリーン関数であることを特徴とする請求項２記載の空力音源探査システム。
【請求項４】
前記第２の演算手段が用いる式は、前記速度ベクトルをｕ、渦度ベクトルをω（ω＝rotｕ）、観測点の座標をｘ、前記速度ベクトルの座標をｙとして、下記式で表されることを特徴とする請求項３記載の空力音源探査システム。
【数１】
2004156968.tif 000001
【請求項５】
前記第１の演算手段は、前記撮影画像を複数の微小領域セルに分割し、各セルごとに前記トレーサ粒子の局所変位を求めることで、前記速度ベクトルの分布を解析することを特徴とする請求項１乃至請求項４いずれか１項に記載の空力音源探査システム。
【請求項６】
流れ場の中の空力音の二重極音源を解析する空力音源探査方法であって、
流れ場を可視化するためにトレーサ粒子を流れ場に混入する可視化工程と、
流れ場をシート状照明光により照明する照明工程と、
シート状に照明された流れ場の測定領域を撮影する撮影工程と、
前記撮影工程により撮影された画像を記録装置に格納する格納工程と、
前記記録装置に格納された画像であって、微小時間間隔で撮影された画像から、流れ場の速度ベクトル分布を解析する第１の解析工程と、
前記第１の解析工程により解析された速度ベクトルの分布から、二重極音の音源分布を解析する第２の解析工程と、
を備えたことを特徴とする空力音源探査方法。
【請求項７】
前記第２の解析工程は、Howeの渦音の式にグリーン関数を適用した式を用いて音源分布を解析する工程であることを特徴とする請求項６記載の空力音源探査方法。
【請求項８】
前記グリーン関数は、コンパクトグリーン関数であることを特徴とする請求項７記載の空力音源探査方法。
【請求項９】
前記第２の解析工程が用いる式は、前記速度ベクトルをｕ、渦度ベクトルをω（ω＝rotｕ）、観測点の座標をｘ、前記速度ベクトルの座標をｙとして、下記式で表されることを特徴とする請求項８記載の空力音源探査方法。
【数１】
2004156968.tif 000002
【請求項１０】
前記第１の解析工程は、前記撮影画像を複数の微小領域セルに分割し、各セルごとに前記トレーサ粒子の局所変位を求めることで、前記速度ベクトルの分布を解析する工程であることを特徴とする請求項６乃至請求項９いずれか１項に記載の空力音源探査方法。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
本発明は、空力音の音源分布を解析するシステム及び方法に関し、特に、ＰＩＶ（Particle Image Velocimetry：粒子画像流速測定法）を用いて、空力音の音源分布を解析するシステム及び方法に関する。
【０００２】
【従来の技術】
近年、流れの速度分布を計測する方法として、ＰＩＶと呼ばれる手法が知られるようになってきている。
【０００３】
通常、流れは目に見えないものであるが、トレーサ粒子（１μｍの油滴等）と呼ばれるマーカーを挿入することによって見えるようにすることが流れの可視化である。この流れの可視化技術にデジタル画像処理技術を加え、流れ場の瞬時・多点の速度情報を抽出する方法がＰＩＶである。ＰＩＶは、流れの速度分布を調べる方法として最も進化した計測法の一つであり、有力な流れ解析手法として威力を発揮している。
【０００４】
図６は、ＰＩＶシステム１の概略構成を示す構成図である。ＰＩＶシステム１は、照明装置２、光シート光学系３、撮影装置４、トレーサ粒子発生装置５から構成される。
【０００５】
測定に際しては、まず、このシステム１自体を測定対象である流れ場の周囲に載置する。続いて、トレーサ粒子発生装置５からトレーサ粒子を発生させ、流れの中に混入させる。そして、照明装置２からレーザ光を発生させ、光シート光学系３を通過させることで、トレーサ粒子を含む流れの中に光シートを照射する。このレーザ光は、短い時間間隔（ｄｔ）で２回照射する。そして、撮影装置４により、上記２回の照射の際のトレーサ粒子を撮影して、２枚の画像データを得る。
【０００６】
このようにして撮影された２枚の画像内の各トレーサ粒子の移動距離が、図示しない解析装置により、２次元のＦＦＴ等の統計的な手法を用いて求められる。各トレーサ粒子の移動距離を、照射時間間隔ｄｔで割ると、速度ベクトルが得られるので、流れの速度分布を求めることができる。
【０００７】
なお、ＰＩＶシステムに関しては、例えば、下記非特許文献１で詳細に説明されている。
【０００８】
【非特許文献１】
可視化情報学会編，「ＰＩＶハンドブック」，森北出版株式会社，２００２年７月
【０００９】
【発明が解決しようとする課題】
上述した従来のＰＩＶシステムは、まだ普及の初期段階であり、基本技術と応用技術の研究・開発が並行して行われている状況である。
【００１０】
その応用としては、自動車エンジン・燃料噴霧・タービンの流れ場の測定、風洞実験の流れ場の測定、構造物周囲・内部の気体の流れ場の測定、水・オイルなど液体の流れ場の測定、ファン周辺の流れ場の測定等が挙げられる。しかし、このＰＩＶシステムによる流れ場の測定値を用いて後述する空力音の解析を行う、総合的な解析システムについては、報告がなされていない。
【００１１】
他方、近年の新幹線の高速化に伴い、沿線騒音に占める空力音の寄与率が大きくなってきている。従来の空力音低減対策法は、実際の車両で問題となっている現象について、条件を変えた模型を用いた風洞実験により、手探り的に行っていた。
【００１２】
この方法は、現実の問題を解決するためには優れた方法である。しかし、反面、実験を行う者の経験によるところが大きい。よって、どういった物理現象が起こっているかを正しく捉えないと、改良の方向性を誤り、有効な低減効果を見込めない。また、近年の新幹線の高速化に伴い、経験則による手探り的な空力音低減対策法の限界も見えてきている。
【００１３】
このような事情から、流体騒音の理論に基づいた空力音低減対策法に対する期待が高まっている。
【００１４】
本発明は、このような課題に鑑みてなされたものであり、ＰＩＶ手法を具体的に空力音低減対策に適用し、空力音の音源分布を解析する空力音源探査システム及び空力音源探査方法を提供することを目的とする。
【００１５】
【課題を解決するための手段】
上記課題を解決するために、本発明に係る空力音源探査システムは、流れ場を可視化するためのトレーサ粒子を発生するトレーサ粒子発生手段と、流れ場をシート状照明光により照明する光シート照明手段と、前記光シート照明手段により照明された所定の測定領域を撮影するための撮影手段と、前記撮影手段により撮影された複数の撮影画像を格納する記録手段と、前記記録手段に記録されている微小時間間隔で撮影された撮影画像から、流れ場の速度ベクトルの分布を解析する第１の演算手段と、前記第１の演算装置により解析された速度ベクトルの分布から、空力音の二重極音の音源分布を解析する第２の演算手段と、を備えたことを特徴とする。
【００１６】
また、本発明に係る空力音源探査システムは、Howeの渦音の式にコンパクトグリーン関数を適用した式を用いて音源分布を解析することを特徴としている。音源がコンパクトであり、観測点が遠方の場合には、Howeの渦音の式にコンパクトグリーン関数を適用した式を用いて解析することになる。
【００１７】
また、上記課題を解決するために、本発明に係る空力音源探査方法は、流れ場を可視化するためにトレーサ粒子を流れ場に混入する可視化工程と、流れ場をシート状照明光により照明する照明工程と、シート状に照明された流れ場の測定領域を撮影する撮影工程と、前記撮影工程により撮影された画像を記録装置に格納する格納工程と、前記記録装置に格納された画像であって、微小時間間隔で撮影された画像から、流れ場の速度ベクトル分布を解析する第１の解析工程と、前記第１の解析工程により解析された速度ベクトルの分布から、二重極音の音源分布を解析する第２の解析工程と、を備えたことを特徴とする。
【００１８】
【発明の実施の形態】
以下、図面を参照しながら本発明の実施の形態について説明する。
図１は、本実施の形態に係る空力音源探査システム１０の全体構成を概略的に示す図である。
【００１９】
図１に示すように、本実施の形態に係る空力音源探査システム１０は、煙発生装置１１、シート照明装置１２、ＣＣＤカメラ１３、及び音源解析装置２０から構成される。音源解析装置２０は、パーソナルコンピュータ（以下、「ＰＣ」とする）２１とディスプレイ２２を含み、ＰＣ２１は図示しない演算装置と記録装置を含む。
【００２０】
また、図１に示すように本実施の形態に係る空力音源探査システムの測定対象は、円柱模型１４である。円柱模型１４は上下面を地面に平行に配置されており、本実施の形態でのサイズは、直径５０ｍｍ、長さ６００ｍｍである。
【００２１】
また、空力音源探査システム１０は、図示しない風洞装置の中に配置されている。そして、風洞装置を作動させると、図中左側から右側に向けて風が吹き、円柱模型１４の周囲には流れ場が構成されるようになっている。本実施の形態に係る風洞装置の風速は、１５０ｋｍ／ｈ（４１．７ｍ／ｓ）である。この流れ場により、円柱模型１４の周囲に空力音が発生しており、本実施の形態に係る空力音源探査システム１０によって、この空力音源の空間分布を解析することになる。
【００２２】
ＰＣ２１は、ＣＣＤカメラ１３に接続されており、ＣＣＤカメラ１３で撮影された画像は、ＰＣ２１内の記録装置に格納される。ＰＣ２１内の演算装置は、この記録装置に格納された画像データに対して、処理・解析を行うよう構成されている。また、その解析結果は、ディスプレイ２２に表示されることになる。
【００２３】
シート照明装置１２は、流れ場の所定の二次元測定領域（図１参照）を、４Ｈｚ周期でシート状の照明光により照明しており、ＣＣＤカメラ１３は、この照明に同期させて撮影を行っている。シート照明装置１２とＣＣＤカメラ１３の同期は、これらに接続された図示しない同期装置により取られている。よって、ＰＣ２１内の記録装置には、微小時間間隔で撮影された時系列の画像が複数格納されることになる。
【００２４】
図２は、ＣＣＤカメラ１３と円柱模型１４の位置関係を示す図であり、図２（ａ）は図１の手前からみた位置関係を示す図、図２（ｂ）は図１の上方からみた位置関係を示す図である。
【００２５】
図２に示すように、ＣＣＤカメラ１３は、測定領域面に対して垂直方向の位置であって、円柱模型１４周囲の所定の測定領域を撮影できる位置に配置されている。本実施の形態では、図中、円柱模型１４の右側（下流側）を撮影できる位置に配置されている。また、ｘ軸、ｙ軸、ｚ軸の原点は、円柱模型１４の重心付近の位置としている。
【００２６】
次に、この空力音源探査システム１０により、空力音源の空間分布を解析する際の動作について説明する。
【００２７】
まず、風洞装置を作動させてから、煙発生装置１１により煙を発生させる。この煙発生装置１１は、図中、円柱模型１４よりも流れの上流側である左側に設置されており、発生した煙は、円柱模型１４周囲の流れ場に混入されることになる。これにより、流れ場が可視化される。この煙は、流れに追従できるような粒子を入れる必要がある。粒子の径が大き過ぎたり、粒子の質量が重過ぎたりすると、流れに追従できなくなるからである。
【００２８】
煙が混入された後、シート照明装置１２によりシート状照明光が円柱模型１４付近の測定領域に照射される。
【００２９】
続いて、照明された瞬間の流れ場の画像をＣＣＤカメラ１３により撮影し、ＰＣ２１内の記録装置に記録する。ＣＣＤカメラ１３により撮影された画像は、音源解析装置２０のＰＣ２１に転送され、記録装置に保存される。
【００３０】
続いて、音源解析装置２０において、記録装置に保存された画像を解析する処理について、図３を参照しながら説明する。図３は、音源解析装置２０での解析処理について説明するためのフロー図である。
【００３１】
まず、ＰＣ２１は、ステップ１（以下、「Ｓ１」とする）において、記録装置に保存された時系列的に隣り合う２枚の画像間での粒子の移動を解析することで、速度ベクトルの分布（流速分布）を算出する。この画像解析では、各画像を３５×３８＝１３３０個のセルに分割し、各セルごとに速度ベクトルｕを算出している。この処理は、従来のＰＩＶシステムにおける処理と同様であるので、詳細な説明は省略する。
【００３２】
得られる速度ベクトルとしては、例えば、ｘ＝１００ｍｍ、ｙ＝１５０ｍｍ、ｚ＝０ｍｍの点で、ｘ方向の瞬時流速＝１０ｍ／ｓ、ｙ方向の瞬時流速＝−５ｍ／ｓ、また、ｘ＝１００ｍｍ、ｙ＝１００ｍｍ、ｚ＝０ｍｍの点で、ｘ方向の瞬時流速＝２０ｍ／ｓ、ｙ方向の瞬時流速＝１０ｍ／ｓ、といった速度ベクトルを得ることができる。
【００３３】
図４は、Ｓ１の画像解析により求められた速度ベクトルｕの分布を示す図である。図４（ａ）は、時刻ｔ1における瞬時の速度分布を示す図であり、図４（ｂ）は、時刻ｔ2における瞬時の速度分布を示す図である。図４においては、矢印が速度ベクトルを示しており、矢印の長さが速度ベクトルの大きさを表す。図中、上下端に位置する水平を向いている速度ベクトルの大きさが、上述した風洞装置の風速１５０ｋｍ／ｈとなる。
【００３４】
図４に示す速度ベクトルの分布から、円柱模型１４の下流側にカルマン渦が発生していることが分かる。
【００３５】
次に、Ｓ２において、Ｓ１で求めた速度ベクトルｕから、渦度ベクトルω（ω＝rot ｕ）を求める。
【００３６】
次に、Ｓ３において、下記（１）式に、速度ベクトルｕ及び渦度ベクトルωを代入することで、空力音の二重極音源（ｍ／ｓ3）を求める。
【００３７】
【数２】
2004156968.tif 000003
また、前記速度ベクトルｕの位置ベクトルをｙ、観測点の位置ベクトルをｘとして、Ｙj（ｙ）は、無限遠方でｊ方向に単位速度を持ち、
【数３】
2004156968.tif 000004
を満たす、仮想的な非圧縮流れの速度ポテンシャルである。さらに、Ｙj（ｙ）は、ラプラス方程式∇2Ｙj（ｙ）＝０を満たす。
【００３８】
また、ｊは、制約のない空間では、ｊ＝１，２，３の値をとり、それぞれｘ，ｙ，ｚ成分を示す。ただし、本実施の形態では、図１に示すように、ｘ−ｙ平面における二次元の音源分布を求めているので、ｊ＝１，２となる。
【００３９】
具体的には、まず、Ｓ１及びＳ２で求めた位置ベクトルｙのセルでの速度ベクトルｕと渦度ベクトルωを、（１）式にそれぞれ代入する。次に、別途求めてあるgradＹj（ｙ）を掛け合わせ、さらに時間微分することで、音源を求めることができる。そして、１３３０個全てのセルに対して、この処理を実行すると、音源分布を求めることができる。
【００４０】
図５は、Ｓ３により求められた、二重極音源のｙ成分の分布を示す図である（Δｔ＝ｔ2−ｔ1）。図中、色の濃い部分が絶対値の大きな音源が存在する部分である。このように、二重極音源が解析できれば、空力音低減のために有効な対策を施すことが可能となる。
【００４１】
ここで、（１）式の算出方法について説明する。下記（２）式に示すHoweの渦音の式と、物体形状に適したコンパクトグリーン関数を組み合わせることで、下記（３）式が算出される。そして、（３）式の被積分関数を時間微分したものが、（１）式である。
【００４２】
【数４】
2004156968.tif 000005
【数５】
2004156968.tif 000006
【００４３】
ここで、Ｐaは位置ベクトルｘの観測点における音圧（Ｐａ）であり、Ｇ（ｘ，ｙ，ｔ−τ）はグリーン関数である。
そして、Howeの渦音の式である（２）式は、以下のようにして求められる。
Lighthillは、連続の式（４）とNavier-Stokes方程式（５）から、Lighthillの流体音響の基礎方程式（６）を導いた。これは、音速ｃ0と密度ρ0が一定である均一な流体の無限の広がりの中で、その一部を占める比較的小さな不均一領域から発生する音場の非同次波動方程式である。
【００４４】
【数６】
2004156968.tif 000007
【数７】
2004156968.tif 000008
【数８】
2004156968.tif 000009
ここで、ｐは瞬時圧力変動、ρは流体密度である。
【００４５】
また、τijは、粘性応力テンソルであり、粘性係数をμとすると、下記式（７）で表わされる。
【数９】
2004156968.tif 000010
但し、δijはディラック関数である。
【００４６】
また、ＴijはLighthillの応力テンソルと呼ばれ、式（８）で表わされる。
【数１０】
2004156968.tif 000011
ここで、ｐ0は均一な流体の基準圧力である。
【００４７】
そして、Howeは、マッハ数が小さい場合について、上記Lighthillの式（６）を近似して、次式（９）を導いた。
【数１１】
2004156968.tif 000012
ここで、Ｂは次式（１０）で表わされる。
【数１２】
2004156968.tif 000013
【００４８】
式（９）の右辺は、流れの中の渦が音源となることを示している。静止した固体面Ｓ（本実施の形態では、円柱模型の表面が該当）が存在するとき、式（９）の解は、グリーンの定理を用いて次式（１１）により求めることができる。
【数１３】
2004156968.tif 000014
ただし、ｎjは固体面Ｓの法線ベクトルを表わし、流体域に向かう方向を正とする。グリーン関数Ｇ（ｘ，ｙ，ｔ−τ）は、次式（１２）を満たし、また固体面Ｓ上で法線方向の微分が０、すなわち、式（１３）を満たして物体形状に適したものを選ぶものとする。
【００４９】
【数１４】
2004156968.tif 000015
【数１５】
2004156968.tif 000016
そして、発散定理を用いて、次式（１４）となる。
【数１６】
2004156968.tif 000017
【００５０】
一方、上記Navie-Stokes方程式（５）を、渦度を用いて表わし、さらに高レイノルズ数流れとして、粘性項を無視すると、式（１５）となる。
【数１７】
2004156968.tif 000018
式（１４）と式（１５）を、式（１２）に入れて整理すると、次式（１６）になる。
【数１８】
2004156968.tif 000019
【００５１】
このとき、式（１６）の第２項は、固体面Ｓが剛体であるならば０となる。
次に、観測点ｘが音源領域から十分離れているとして、式（１０）を次式（１７）のように近似すると、渦と固体面Ｓの干渉による圧力変動（空力音）ｐaは、式（２）で表わされる。この式（２）が、Howeの渦音の式である。
【数１９】
2004156968.tif 000020
【００５２】
次に、本実施の形態のように、音源がコンパクトであり、また、観測点が音源から十分遠方にある場合には、グリーン関数Ｇを次式（１８）で表わされるコンパクトグリーン関数Ｇcに置き換えることができる。
【数２０】
2004156968.tif 000021
【００５３】
ここで、音源がコンパクトとは、音波の波長の範囲内に音源が集まっていることをいい、観測点が音源から十分遠方であるとは、音波の波長の１／４以上離れていることをいう。
【００５４】
このコンパクトグリーン関数（１８）を、Howeの渦音の式（２）に代入することで、上述した式（３）を得ることができる。式（３）によって抽出される音は、物体形状によって四重極音が散乱されたことに伴う二重極音に相当する。よって、式（３）の被積分関数から二重極音源の分布を求めることができる。
【００５５】
以上、詳細に説明したが、本実施の形態によれば、円柱模型の周囲に発生する空力音の二重極音源の分布を求めることができる。このように物理現象を正確に把握すれば、空力音低減対策において、流体騒音の理論に基づいた低減対策が可能となり、従来の経験則に基づく騒音対策に比較して有効な騒音対策を施すことができる。
【００５６】
また、本実施の形態によれば、ＰＩＶシステムの測定値を用いて二重極音源分布を解析しており、剥離、乱流遷移、壁面近くでの境界層の挙動など、流体の数値計算等の理論的な方法では精度良く捉えることの困難なものであっても、十分な精度で解析することが可能である。また、ＰＩＶシステムの測定値を用いることで、キャビティ音のように、発生した圧縮波が流体の挙動に影響を与えるような現象についても、精度良く解析することができる。
【００５７】
なお、本発明の実施の形態は、上述した実施の形態に限定されるものではなく、本発明の主旨を逸脱しない範囲で種々の変形が可能である。
【００５８】
例えば、本実施の形態においては、二次元の空力音源分布を求めているが、ＰＩＶシステムとして、三次元の速度分布を計測可能な公知のステレオＰＩＶシステムを使用すれば、三次元の空力音源分布を求めることも可能である。
【００５９】
また、説明を簡単にするため、円柱模型周りの空力音源分布を求めているが、円柱模型に限られるものでなく、列車模型、自動車模型等の、空力音の発生を解析する対象であれば、どのような形状の模型であっても良い。この場合、ラプラス方程式を満たす式（１）のＹを形状に応じて数値的に求めれば良い。
【００６０】
また、本実施の形態に係るＰＩＶシステムは、４Ｈｚ周期で流れ場を撮影しているが、リアルタイムで音源分布を解析するために、リアルタイムＰＩＶシステムを用いても良い。例えば、撮影の周期を１ｋＨｚ以上とすることで、リアルタイムＰＩＶシステムを実現することができる。
【００６１】
【発明の効果】
以上、本発明によれば、ＰＩＶシステムを用いて空力音の音源分布を解析する、空力音源探査システム及び空力音源探査方法を提供することができる。
【図面の簡単な説明】
【図１】本発明の実施の形態に係る空力音源探査システムの概略構成を示す図である。
【図２】本発明の実施の形態に係るＣＣＤカメラと円柱模型との位置関係を示す図である。
【図３】本発明の実施の形態に係る音源解析装置での画像解析処理を説明するためのフロー図である。
【図４】本発明の実施の形態に係る速度ベクトルの分布を示す図である。
【図５】本発明の実施の形態に係る二重極音源の分布を示す図である。
【図６】従来のＰＩＶシステムの概略構成を示す図である。
【符号の説明】
１０      空力音源探査システム
１１      煙発生装置
１２      シート照明装置
１３      ＣＣＤカメラ
１４      円柱模型
２０      音源解析装置
２１      ＰＣ
２２      ディスプレイ