JPH0670638B2 - 流れ場データの格子点補間方法 - Google Patents
流れ場データの格子点補間方法Info
- Publication number
- JPH0670638B2 JPH0670638B2 JP30047390A JP30047390A JPH0670638B2 JP H0670638 B2 JPH0670638 B2 JP H0670638B2 JP 30047390 A JP30047390 A JP 30047390A JP 30047390 A JP30047390 A JP 30047390A JP H0670638 B2 JPH0670638 B2 JP H0670638B2
- Authority
- JP
- Japan
- Prior art keywords
- flow velocity
- data
- velocity data
- extracted
- points
- Prior art date
- Legal status (The legal status 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 status listed.)
- Expired - Lifetime
Links
Landscapes
- Indicating Or Recording The Presence, Absence, Or Direction Of Movement (AREA)
- Aerodynamic Tests, Hydrodynamic Tests, Wind Tunnels, And Water Tanks (AREA)
Description
次元断面内の任意の場所で測定した2次元流速データか
ら等間隔(又は任意な間隔:非等間隔)の格子点上の2
次元流速データを算出する方法に関する。
計等を用いて測定を行う場合、予め測定点を測定空間の
2次元断面内で測定者が決定している。測定点は、等間
隔又は非等間隔の格子点として決定し、各格子点におい
て流速の測定を行うことで流れ場の測定を行う。これら
格子点で測定された2次元流速データを基にして例えば
楕円形ポアソン方程式を解析して流れ関数を算出し、流
れ関数に基いて各種の物理量の把握を行う。
が膨大となる欠点がある。そこで、可視粒子を流れ場に
混入して、その軌跡の画像データから流れ場を算出する
画像処理方式も知られている。この画像処理方式は、例
えば論文「日本機械学会論文(B編)」第55巻、509(1
989−1)、107〜114頁に開示されていて、流れ場にト
レーサ粒子を混入し、そのトレーサ粒子に連続光又はス
トロボ光を当てて、その軌跡を画像処理する。
流速データが等間隔(又は非等間隔)の格子点上に並ぶ
ことがなく、処理画面内の不特定な位置に流速データが
生成されている。従って、測定断面内の各種物理量(例
えば流れ関数等)を例えば楕円形ポアソン方程式を解析
して算出するに当たっては、測定された流速データを等
間隔(又は非等間隔)の格子点配列に補間演算すること
が不可欠である。
関係な点において得られている2次元流速データから、
等間隔(又は非等間隔)な格子点における2次元流速デ
ータを比較的簡単に算出して、流れ場の流れ関数の算出
を短時間に行えるようにすることを目的とする。
タを表す不規則な分布の測定点における2次元流速デー
タから、規則性のある格子点における2次元流速データ
を生成するデータ補間方法であって、補間すべき格子点
の近傍の3個以上の測定点における流速データを抽出
し、抽出された流速データから3個の流速データを取出
す全ての組合せについて下式に基づく2次元流速ベクト
ルを算出し、算出された各組の流速ベクトルを統計
処理して補間すべき格子点における流速データを近似し
て得ることを特徴とする流れ場データの格子点補間方
法: K、L、Mは抽出された流速データから取り出された3
個の流速データに係る測定点、(dx)K、(dy)K、
(dx)L、(dy)L、(dx)M、(dy)MはXY座標系で
表現したときの各測定点K、L、Mと補間すべき格子点
との間の夫々X軸方向及びY軸方向の距離、∂/∂
x、∂/∂yは補間すべき格子点における夫々X軸方
向及びY軸方向の流速勾配、K、L、Mは測定点
K、L、Mにおける流速データ。
補間ができる。扱うデータ量は多くなるが、高速計算機
を利用すれば比較的短時間に処理することができる。抽
出された流速データを10個とした場合、そのうちの3つ
を取り出す組合せは10C3通りである。全ての組合せにつ
いての補間すべき格子点における流速ベクトルの統計処
理は、例えば相加平均の計算や分布関数を設定して近似
計算により行う。
を示す図であり、第2図は画像処理装置の構成例を示す
ブロック図である。
を実施するための装置は、アルゴンレーザ1、光ファイ
バー2、ピンホール3、光学レンズ4、肉厚10mmのアク
リル製容器5、CCDカメラ8、録画再生装置9(12)、
録画再生装置9(12)に接続されたモニター10(11)、
計算機15、及び、計算機15に接続されたフレームメモリ
14、モニター16及び外部記憶装置17から構成されてい
る。
ーザ光は一度光ファイバー2に入射し、この光ファイバ
ー2の先端に取り付けられた光学レンズ4にて厚み2mm
以下のレーザ光シートを作り、アクリル製容器5内の流
体に照射される。本発明の実施例で用いた液体は水であ
る。
光ファイバー内での光の損失を少なくし、先端の光学レ
ンズ4と光ファイバー2との間に、光ファイバー2の先
端から30mmの位置に内径1mmのピンホール3を設けてレ
ーザ光シート厚みを調整した。
り、このアクリル製容器5の図中手前の面(従って、
上)が大気開放になっている。このアクリル製容器5に
は図中上面より連続的に水が供給されており、図中下面
から連続的に排水されている。
に、水の流れに追従する粒子をアクリル製容器5の上面
から連続的に入れ、容器5内の可視化したい断面に光学
レンズ4からアルゴンレーザを光シート状に連続的に照
射した。そして、レーザシート面1Aと直交る方向にCCD
カメラ8をおき、上記粒子がレーザシート面1Aを通過す
る時に散乱若しくは蛍光した光をそのシート面と直交す
る方向からCCDカメラ8にて撮像し、録画再生装置9に
連続的に録画することでシート面内の映像を録画した。
なお、録画時には、録画再生装置9に接続したモニター
10で録画状態を確認した。
面をメチレンブルー等の染料で表面改質して散乱若しく
は蛍光強度を強くしたものを用いるのが望ましい。本例
では、表面改質をしない粒子径30μm、密度1g/cm3の
粒子を用いた。
秒単位のフレーム情報を、計算機15に接続したフレーム
メモリ14を通し、外部記憶装置17に記録した。
録された例えば4時刻分の画像情報は、計算機15内ソフ
トウェアで1/60秒のフィールド情報に分割され、連続
した4枚の例えば奇数フィールドのみの画像情報が抽出
される。抽出された画像情報は、一旦、フレームメモリ
14に格納される。
により画像処理がなされる。このステップS2の画像処理
は、ノイズ処理、2値化及び粒子の重心の算出である。
ルド間で移動している粒子の対応付けが行われる。この
対応付けは、例えば前記論文に開示された粒子トレース
方法、即ち流速の連続性に基いて4時刻の粒子の対応付
けを行う方法を適用することができる。
次元流速データが算出され、ベクトル表示が行われる。
第4図に、算出された流速データのベクトル表示を示
す。
を基に測定平面上の格子点における2次元流速データを
補間により生成する。生成された流速データは、例えば
第5図に示すように等間隔の格子点上に位置する。この
ようにして生成された流速データは、例えば楕円形ポア
ソン方程式に適用され、流れ関数が算出される。
ルゴリズムを説明するフローチャートである。
作る。第7図は格子点の一例を示し、この例では等間隔
で互いに直交する直線を仮定し、その交点を格子点とす
る。この例では、各直線は等間隔であるが、非等間隔で
もよい。例えば、着目している部分とその周囲とで格子
間隔を変えてもよく、或いは、同心円と放射直線から成
る極座標を想定してもよい。
き格子点Gから一定範囲(半径lの円)内のn個の流速
データの測定点(夫々をK、L、M、N……とする)に
夫々対応したn個の2次元流速データを抽出する。lは
例えば格子点間距離の6倍とする。nは補間すべき格子
点Gに近い測定点から順次取り出した一定個数p(3個
以上で、例えば10個)とする。
個に達しているか否かを判定する。ここで、抽出された
流速データ数nがp個に足りないと判定された場合に
は、ステップS13で、抽出範囲の半径lを大きくしてス
テップS11を再度行う。なお、ステップS13でpの値を例
えば8個、6個…のように順次低減してもよい。
していると判定された場合には、ステップS14で、n個
の流速データから3個の流速データを取り出して得られ
る全部でnC3通りの組合せの1組(例えば3個の測定点
K、L、Mにおける流速データ)を用いて、格子点Gに
おける流速ベクトルを補間により算出する。補間は、半
径lの抽出円内での夫々X軸方向及びY軸方向の流速勾
配を一定と仮定し、この円内での流速勾配が格子点Gに
おける流速勾配に等しいとした線形補間で行う。即ち、
第8図において、補間すべき格子点Gの流速ベクトル
(2次元ベクトル)を=(ux,uy)とし(なお、本明
細書中、頭に矢印がついた文字は全て2次元ベクトルを
表す。)、測定点Kの流速データ(2次元ベクトル)を
K=(uKx,uKy)とし、(dx)K、(dy)Kを測定点
Kと格子点Gとの間の夫々X軸方向及びY軸方向の距離
とし、∂/∂x、∂/∂yを格子点Gでの夫々X軸
方向及びY軸方向の流速勾配、即ち、∂/∂x=(∂
ux/∂x,∂uy/∂x)、∂/∂y=(∂ux/∂y,∂uy
/∂y)とすると、 という3つの未知量(夫々が2次元ベクトルである、
∂/∂x、∂/∂y)を含む方程式が成り絶つ。他
の2つの測定点L、M(各測定点での流速データを夫々
L、Mとし、格子点GとのX軸方向及びY軸方向の
距離を夫々(dx)L、(dx)M、(dy)L、(dy)Mと
する)についても夫々同様の線形式が成り立つので、下
記の連立方程式(2)が得られる。この連立方程式
(2)を解くことにより、求める流速ベクトル(2次
元ベクトル)が得られる。
1,2…nC3)が得られた否かを判定し、nC3個の流速ベク
トルiが得られるまでnC3個の全ての組合せの夫々に
ついてステップS14を繰り返して実行させる。
1〜nC3)の平均値を統計処理して、補間すべき格子点
Gにおける流速データを近似して求める。統計処理は、
例えば相加平均計算Σi/nC3でよい。或いは、第9
図に示すように、流速ベクトルiの度数Nを求めてガ
ウス近似又は二項近似により平均値を求める統計計算を
行ってもよい。尚、第9図では、図面を簡略にするため
に、2次元流速ベクトルiを1次元の量として表示し
てある。ここまでの処理によって、格子点Gにおける2
次元流速データを生成する補間計算が終了する。
る。
補間計算の終了を検出するまで、ステップS11〜S17の処
理を繰り返す。
おける2次元流速データを補間計算することができる。
この補間された格子点での流速データを用いて、流れ関
数を例えば楕円形ポアソン方程式を解析して算出するこ
とができ、更に流れ関数を基にして流れ場の各種物理量
を求めることができる。第10図は、一例として隣接する
4つの格子点(I、J)1-4での流体の流量洩れΔQを
計算して表示した結果である。
おける2次元流速データから、規則性のある等間隔又は
非等間隔の格子点上の2次元流速データを補間計算によ
り高精度に生成することができる。従って、生成された
流速データから例えば楕円形ポアソン方程式を解析する
ことにより流れ関数を算出でき、流れ場の物理的把握が
短時間でできるようになる。特に流れ場を表す不規則な
分布の測定点における流速データとして、複数の時刻に
おける画像データから画像処理によって得たデータを用
いることにより、流れ場の物理量の把握が短時間に高精
度でできるようになる。
像処理装置の構成例を示すブロック回路図、第3図は画
像処理手順のフローチャート、第4図は画像処理によっ
て得られた測定点における流速データの分布図、第5図
は補間計算によって得られた格子点における流速データ
の分布図、第6図は補間処理の手順を示すフローチャー
ト、第7図は補間すべき格子点を示す補間領域の図、第
8図は補間アルゴリズムを説明するための格子点及び測
定点の図、第9図は補間計算中の統計処理の一例を示す
グラフ、第10図は格子点における流速データを基に計算
した流量洩れのグラフである。 なお、図面に用いた符号において、 1……アルゴンレーザ 2……光ファイバー 4……光学レンズ 5……アクリル容器 8……CCDカメラ 9、12……録画再生装置 14……フレームメモリ 15……計算機 16……モニター 17……外部記憶装置 である。
Claims (6)
- 【請求項1】流れ場データを表す不規則な分布の測定点
における2次元流速データから、規則性のある格子点に
おける2次元流速データを生成するデータ補間方法であ
って、 補間すべき格子点の近傍の3個以上の測定点における流
速データを抽出し、抽出された流速データから3個の流
速データを取出す全ての組合せについて下式に基づく2
次元流速ベクトルを算出し、算出された各組の流速ベ
クトルを統計処理して補間すべき格子点における流速
データを近似して得ることを特徴とする流れ場データの
格子点補間方法: K、L、Mは抽出された流速データから取り出された3
個の流速データに係る測定点、(dx)K、(dy)K、
(dx)L、(dy)L、(dx)M、(dy)MはXY座標系で
表現したときの各測定点K、L、Mと補間すべき格子点
との間の夫々X軸方向及びY軸方向の距離、∂/∂
x、∂/∂yは補間すべき格子点における夫々X軸方
向及びY軸方向の流速勾配、K、L、Mは測定点
K、L、Mにおける流速データ。 - 【請求項2】上記格子点が、互いに直交する等間隔又は
非等間隔の複数の直線の交点として設定されていること
を特徴とする請求項1の補間方法。 - 【請求項3】上記格子点が、複数の同心円及び放射能直
線の交点として設定されていることを特徴とする請求項
1の補間方法。 - 【請求項4】上記格子点から所定の距離の範囲内に存在
する3個以上の測定点での流速データを抽出することを
特徴とする請求項1の補間方法。 - 【請求項5】上記統計処理が相加平均計算であることを
特徴とする請求項1の補間方法。 - 【請求項6】上記統計処理が分布関数に基づいた近似計
算であることを特徴とする請求項1の補間方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP30047390A JPH0670638B2 (ja) | 1990-11-06 | 1990-11-06 | 流れ場データの格子点補間方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP30047390A JPH0670638B2 (ja) | 1990-11-06 | 1990-11-06 | 流れ場データの格子点補間方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
JPH04174362A JPH04174362A (ja) | 1992-06-22 |
JPH0670638B2 true JPH0670638B2 (ja) | 1994-09-07 |
Family
ID=17885221
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP30047390A Expired - Lifetime JPH0670638B2 (ja) | 1990-11-06 | 1990-11-06 | 流れ場データの格子点補間方法 |
Country Status (1)
Country | Link |
---|---|
JP (1) | JPH0670638B2 (ja) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP5027293B2 (ja) * | 2007-04-04 | 2012-09-19 | テレフオンアクチーボラゲット エル エム エリクソン(パブル) | 電気通信における無線ネットワーク計画、シミュレーションおよび解析を改善するための方法および装置 |
-
1990
- 1990-11-06 JP JP30047390A patent/JPH0670638B2/ja not_active Expired - Lifetime
Also Published As
Publication number | Publication date |
---|---|
JPH04174362A (ja) | 1992-06-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Pereira et al. | Defocusing digital particle image velocimetry and the three-dimensional characterization of two-phase flows | |
JP2013088898A (ja) | 3dデータ解析のための装置、方法及びプログラムと、微小粒子解析システム | |
CN103376367A (zh) | 放射信号可视化装置 | |
JP6136018B2 (ja) | 可視化流体の流速計測方法及び流速計測システム | |
CN107247926B (zh) | 一种人体检测方法及装置 | |
CN106062824A (zh) | 边缘检测装置、边缘检测方法和程序 | |
Bourgoin et al. | Using ray-traversal for 3D particle matching in the context of particle tracking velocimetry in fluid mechanics | |
JPWO2017104662A1 (ja) | 粒子分析装置及び粒子分析方法 | |
US10324021B2 (en) | Magnetophorisis measuring system for determining motion status of object and quantifying amount of magnetic particles contained therein | |
Pu et al. | Off-axis holographic particle image velocimetry for diagnosing particulate flows | |
US3621130A (en) | System for quantizing graphic displays | |
Zhu et al. | Volumetric resolution of light field imaging and its effect on the reconstruction of light field PIV | |
Smith et al. | A technique for evaluation of three-dimensional behavior in turbulent boundary layers using computer augmented hydrogen bubble-wire flow visualization | |
Sheu et al. | A three-dimensional measurement technique for turbulent flows | |
JPH0670638B2 (ja) | 流れ場データの格子点補間方法 | |
JP2013088167A (ja) | 3dデータ解析のための装置、方法及びプログラムと、微小粒子解析システム | |
JP5691452B2 (ja) | 3dデータ解析装置および3dデータ解析方法ならびに3dデータ解析プログラム | |
JPH0663905B2 (ja) | 流れ場データの格子点補間方法 | |
JP2681745B2 (ja) | レーザ光を利用したスペックルパターンによる被計測物の上下および横移動量の測定方法。 | |
JPH0635944B2 (ja) | 画像処理による流れの可視化方法 | |
JP2879501B2 (ja) | 可視化画像のトレーサ追跡方法 | |
Koenderink et al. | MARVIN: high speed 3d imaging for seedling classification | |
Bawden | Applications of aerial photography in land system mapping | |
Kawasue et al. | 3-D measurement of moving particles by circular image shifting | |
JPH0663906B2 (ja) | 画像処理による流れの可視化方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
FPAY | Renewal fee payment (prs date is renewal date of database) |
Year of fee payment: 13 Free format text: PAYMENT UNTIL: 20070907 |
|
FPAY | Renewal fee payment (prs date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20080907 Year of fee payment: 14 |
|
FPAY | Renewal fee payment (prs date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090907 Year of fee payment: 15 |
|
FPAY | Renewal fee payment (prs date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100907 Year of fee payment: 16 |
|
FPAY | Renewal fee payment (prs date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100907 Year of fee payment: 16 |
|
FPAY | Renewal fee payment (prs date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110907 Year of fee payment: 17 |
|
EXPY | Cancellation because of completion of term | ||
FPAY | Renewal fee payment (prs date is renewal date of database) |
Year of fee payment: 17 Free format text: PAYMENT UNTIL: 20110907 |