JPH063365A - 二相流計測装置 - Google Patents
二相流計測装置Info
- Publication number
- JPH063365A JPH063365A JP16573692A JP16573692A JPH063365A JP H063365 A JPH063365 A JP H063365A JP 16573692 A JP16573692 A JP 16573692A JP 16573692 A JP16573692 A JP 16573692A JP H063365 A JPH063365 A JP H063365A
- Authority
- JP
- Japan
- Prior art keywords
- probe
- distribution
- bubble
- velocity
- bubbles
- 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.)
- Pending
Links
Landscapes
- Indicating Or Recording The Presence, Absence, Or Direction Of Movement (AREA)
Abstract
(57)【要約】
【構成】計測対象である配管５の中を流れる水の中の気
泡６の速度分布と気泡径分布を測定するため、電気抵抗
式二探針プローブ３，計測器１０，演算器２０，表示器
３０と入力器３１を構成する。データ処理部２２では、
上流側プローブ１と下流側プローブ２の信号から、気泡
がプローブに接触している時間と、気泡が上流側と下流
側のプローブに接触する時刻の差を算出し、それらの頻
度分布を求める。一方、モンテカルロ演算部２３では、
気泡の速度分布，気泡径分布を仮定して、頻度分布を求
め、二つの頻度分布を比較計算部２４で比較し、仮定を
修正して、頻度分布が十分一致したとみなせた時の気泡
の速度分布，気泡径分布を表示器３０に表示する。 【効果】気液二相流中の気泡や液滴の実際の速度分布と
直径分布とを、精度良く検出できる。
泡６の速度分布と気泡径分布を測定するため、電気抵抗
式二探針プローブ３，計測器１０，演算器２０，表示器
３０と入力器３１を構成する。データ処理部２２では、
上流側プローブ１と下流側プローブ２の信号から、気泡
がプローブに接触している時間と、気泡が上流側と下流
側のプローブに接触する時刻の差を算出し、それらの頻
度分布を求める。一方、モンテカルロ演算部２３では、
気泡の速度分布，気泡径分布を仮定して、頻度分布を求
め、二つの頻度分布を比較計算部２４で比較し、仮定を
修正して、頻度分布が十分一致したとみなせた時の気泡
の速度分布，気泡径分布を表示器３０に表示する。 【効果】気液二相流中の気泡や液滴の実際の速度分布と
直径分布とを、精度良く検出できる。
Description
【０００１】
【産業上の利用分野】本発明は、複数の探針を用いて、
気液二相流中の気泡あるいは液滴の速度分布と、気泡径
あるいは液滴径分布を計測する二相流計測装置に関す
る。
気液二相流中の気泡あるいは液滴の速度分布と、気泡径
あるいは液滴径分布を計測する二相流計測装置に関す
る。
【０００２】
【従来の技術】気液二相流は、化学プラント，発電プラ
ント，各種流体機械など、多くの装置で発生する。これ
らの装置の制御や性能向上のためには、流れの特性を把
握する必要がある。二相流は、条件によって、気泡流，
スラグ流，環状流，噴霧流など、流動様式が複雑に変化
し、測定が比較的容易な圧力損失や平均ボイド率といっ
たグロスな量では、必ずしも現象を正しく把握できない
場合がある。
ント，各種流体機械など、多くの装置で発生する。これ
らの装置の制御や性能向上のためには、流れの特性を把
握する必要がある。二相流は、条件によって、気泡流，
スラグ流，環状流，噴霧流など、流動様式が複雑に変化
し、測定が比較的容易な圧力損失や平均ボイド率といっ
たグロスな量では、必ずしも現象を正しく把握できない
場合がある。
【０００３】気泡流とスラグ流で重要な測定対象は、気
泡の寸法と速度である。環状流，噴霧流で重要な測定対
象は、液膜及び液滴の寸法と速度である。分散相といわ
れる気泡と液滴の寸法や速度を局所的に計測するには、
光プロ−ブ，電気抵抗式プローブ，ホットワイヤなど探
針部を持つ接触測定法が主に用いられている。
泡の寸法と速度である。環状流，噴霧流で重要な測定対
象は、液膜及び液滴の寸法と速度である。分散相といわ
れる気泡と液滴の寸法や速度を局所的に計測するには、
光プロ−ブ，電気抵抗式プローブ，ホットワイヤなど探
針部を持つ接触測定法が主に用いられている。
【０００４】接触測定法による一般的な速度の測定法に
は、例えば世古口らが“ベロシティメジャメント ウィ
ズ エレクトリカル ダブル−センシング デバイス
イントウ−フェイズ フロー”メジャリング テクニッ
ク イン ガス−リクィドトウ−フェイズ フローズ
（“Velocity Measurement with Electorical Double−
Sensing Devices in Two−Phase Flow”，Measuring Te
chniques in Gas− Liquid Two−Phase Flows, IUTAM
Symposium, France（１９８３））で述べているよう
に、上流側と下流側に設置した一対のプローブの出力信
号の相関や出力信号の時間遅れから、平均速度を求める
手法がある。この場合、平均速度は比較的精度良く求ま
るが、計測された速度分布は実際の速度分布とは一致し
ない。
は、例えば世古口らが“ベロシティメジャメント ウィ
ズ エレクトリカル ダブル−センシング デバイス
イントウ−フェイズ フロー”メジャリング テクニッ
ク イン ガス−リクィドトウ−フェイズ フローズ
（“Velocity Measurement with Electorical Double−
Sensing Devices in Two−Phase Flow”，Measuring Te
chniques in Gas− Liquid Two−Phase Flows, IUTAM
Symposium, France（１９８３））で述べているよう
に、上流側と下流側に設置した一対のプローブの出力信
号の相関や出力信号の時間遅れから、平均速度を求める
手法がある。この場合、平均速度は比較的精度良く求ま
るが、計測された速度分布は実際の速度分布とは一致し
ない。
【０００５】気泡や液滴の寸法の測定法には、Hinataが
“ア スタディ オン ザ メジャメント オブ バブ
ル サイズ ディストリビューション イン ザ アプ
ワード バブリィ フロー”（“A Study on the Measu
rement of Bubble Size Distribution in the Upward B
ubbly Flow”，Bulletin of ＪＳＭＥ，２７，ｐ１１２
７(１９８４)）で述べているように、流れの上流側と下
流側、あるいは流れに対して直角な方向に設置した一対
のプローブの出力信号を処理して得られる確率密度と、
分布と気泡径を仮定したときに解析的に求まる確率密度
とを比較して求める手法がある。この手法の欠点は、解
析的に扱うために、気泡の形状を球とするなどの多数の
仮定が必要なことであり、仮定と実際との差が精度の低
下をもたらし、またその補正量を見積もることが困難な
ことである。
“ア スタディ オン ザ メジャメント オブ バブ
ル サイズ ディストリビューション イン ザ アプ
ワード バブリィ フロー”（“A Study on the Measu
rement of Bubble Size Distribution in the Upward B
ubbly Flow”，Bulletin of ＪＳＭＥ，２７，ｐ１１２
７(１９８４)）で述べているように、流れの上流側と下
流側、あるいは流れに対して直角な方向に設置した一対
のプローブの出力信号を処理して得られる確率密度と、
分布と気泡径を仮定したときに解析的に求まる確率密度
とを比較して求める手法がある。この手法の欠点は、解
析的に扱うために、気泡の形状を球とするなどの多数の
仮定が必要なことであり、仮定と実際との差が精度の低
下をもたらし、またその補正量を見積もることが困難な
ことである。
【０００６】
【発明が解決しようとする課題】気液二相流を利用した
装置の制御や性能向上のために、従来にも増して速度分
布や気泡径分布といったミクロの情報が重要になってい
る。従来の手法では、主に気泡や液滴の平均速度や平均
寸法に着目していたため、分布については考慮していな
かったり、低い精度で満足していた。
装置の制御や性能向上のために、従来にも増して速度分
布や気泡径分布といったミクロの情報が重要になってい
る。従来の手法では、主に気泡や液滴の平均速度や平均
寸法に着目していたため、分布については考慮していな
かったり、低い精度で満足していた。
【０００７】本発明の目的は、気泡や液滴の実際の速度
分布と直径分布とを精度良く検出する方法を提供するこ
とにある。
分布と直径分布とを精度良く検出する方法を提供するこ
とにある。
【０００８】
【課題を解決するための手段】本発明は、上記目的を達
成するために、気泡や液滴が二個のプローブに接触する
時刻の差の分布と、気泡や液滴がプローブに接触してい
る時間の分布を測定し、一方で、気泡や液滴の形状，速
度や寸法の分布を仮定したモンテカルロ計算により、プ
ローブから得られる信号をシミュレートし、両者を比較
して、気泡や液滴の実際の速度分布と気泡径分布を算出
する。
成するために、気泡や液滴が二個のプローブに接触する
時刻の差の分布と、気泡や液滴がプローブに接触してい
る時間の分布を測定し、一方で、気泡や液滴の形状，速
度や寸法の分布を仮定したモンテカルロ計算により、プ
ローブから得られる信号をシミュレートし、両者を比較
して、気泡や液滴の実際の速度分布と気泡径分布を算出
する。
【０００９】
【作用】モンテカルロ法により気泡や液滴がプローブに
接触した信号を模擬できる原理を説明する。
接触した信号を模擬できる原理を説明する。
【００１０】流れの上流側と下流側に設置した一対のプ
ローブを考える。プローブの間隔を十分小さく（３mm程
度）すれば、気泡や液滴は上流側プローブと下流側プロ
ーブの間を直線的に移動するとみなせる。以下、簡単の
ため、球形の気泡の場合を例に説明する。
ローブを考える。プローブの間隔を十分小さく（３mm程
度）すれば、気泡や液滴は上流側プローブと下流側プロ
ーブの間を直線的に移動するとみなせる。以下、簡単の
ため、球形の気泡の場合を例に説明する。
【００１１】図２に示すように、上流側プローブ，下流
側プローブの位置を各々(ｘ1,ｙ1,ｚ1)，(ｘ2,ｙ2,
ｚ2)、気泡の中心位置を(ｘ′,ｙ′,ｚ′）、気泡の半
径をｒとする。流れの中の任意の点(ｘ,ｙ,ｚ)につい
て、関数ｆ(ｘ,ｙ,ｚ)を以下のように定義する。
側プローブの位置を各々(ｘ1,ｙ1,ｚ1)，(ｘ2,ｙ2,
ｚ2)、気泡の中心位置を(ｘ′,ｙ′,ｚ′）、気泡の半
径をｒとする。流れの中の任意の点(ｘ,ｙ,ｚ)につい
て、関数ｆ(ｘ,ｙ,ｚ)を以下のように定義する。
【００１２】
【数１】
【００１３】気泡の表面では、ｆ(ｘ,ｙ,ｚ)＝０であ
る。したがって、気泡表面が上流側プローブに接触した
ときには
る。したがって、気泡表面が上流側プローブに接触した
ときには
【００１４】
【数２】
【００１５】が成り立つ。
【００１６】図３に示すように角度θ，ωを定義し気泡
速度をＵb とすると、気泡中心位置(ｘ′,ｙ′,ｚ′）
は、初期の中心位置を(ｘ0,ｙ0,ｚ0)、時刻をｔとして
次式で表せる。
速度をＵb とすると、気泡中心位置(ｘ′,ｙ′,ｚ′）
は、初期の中心位置を(ｘ0,ｙ0,ｚ0)、時刻をｔとして
次式で表せる。
【００１７】
【数３】
【００１８】
【数４】
【００１９】
【数５】
【００２０】以上の式からｘ′，ｙ′，ｚ′を消去して
次式を得る。
次式を得る。
【００２１】
【数６】
【００２２】この式は、時刻ｔに関する二次式であり、
ａｔ2＋ｂｔ＋ｃ＝０ と表せる。したがって、ｂ2−４
ａｃ≧０ であれば、気泡が上流側プローブに接触する
時刻ｔ1 及び、離脱する時刻ｔ2 は、それぞれ
ａｔ2＋ｂｔ＋ｃ＝０ と表せる。したがって、ｂ2−４
ａｃ≧０ であれば、気泡が上流側プローブに接触する
時刻ｔ1 及び、離脱する時刻ｔ2 は、それぞれ
【００２３】
【数７】
【００２４】
【数８】
【００２５】となる。気泡が下流側プローブに接触する
時刻ｔ3 ，離脱する時刻ｔ4 を同様にして求めれば、気
泡が上流側プローブに接触している時間ｔtip ，気泡が
下流側プローブに接触している時間ｔrear，気泡が上流
側プローブに接触してから下流側プローブに接触するま
での時間遅れｔtrは、
時刻ｔ3 ，離脱する時刻ｔ4 を同様にして求めれば、気
泡が上流側プローブに接触している時間ｔtip ，気泡が
下流側プローブに接触している時間ｔrear，気泡が上流
側プローブに接触してから下流側プローブに接触するま
での時間遅れｔtrは、
【００２６】
【数９】
【００２７】
【数１０】
【００２８】
【数１１】
【００２９】となる。
【００３０】乱数を用いて流速Ｕb や角度θ，ω、気泡
径ｒの分布を発生させ、プローブの信号と比較すること
で、二探針プローブの計測信号のシミュレーションがで
きる。なお、ここではプローブの先端は点であるとした
が、実際と同じようにある有限の面積を有しているとし
て評価することもできる。また、気泡形状を球形とした
が、実際の気泡形状に近い楕円体として評価することも
できる。
径ｒの分布を発生させ、プローブの信号と比較すること
で、二探針プローブの計測信号のシミュレーションがで
きる。なお、ここではプローブの先端は点であるとした
が、実際と同じようにある有限の面積を有しているとし
て評価することもできる。また、気泡形状を球形とした
が、実際の気泡形状に近い楕円体として評価することも
できる。
【００３１】
【実施例】以下、本発明の装置の一実施例を図１により
説明する。図１は、計測対象である配管５の中を流れる
液体（水）中の気泡６の速度分布と気泡径分布を測定す
る本発明の装置の一実施例である。装置は、大別して、
電気抵抗式二探針プローブ３，計測器１０，演算器２
０，表示器３０と入力器３１の五つから構成する。電気
抵抗式二探針プローブ３は腐食に強いステンレス鋼や白
金でできた上流側プローブ１，下流側プローブ２からな
り、それぞれ、定圧電源１１ａ，１１ｂと抵抗１２ａ，
１２ｂを有する計測器１０に接続している。気泡は電気
抵抗が大きく、液相（水）は溶解しているイオンにより
電気抵抗が小さいため、プローブが気泡に接したときの
出力端子１３の電位は高く、液相に接したときの出力端
子１３の電位は低くなる。出力端子１３ａ，１３ｂの信
号は、Ａ／Ｄ変換器２１を内蔵した演算器２０に取り込
まれ、データ処理部２２に送られる。データ処理部２２
では、上流側プローブ１と下流側プローブ２の信号か
ら、気泡がプローブに接触している時間と、気泡が上流
側と下流側のプローブに接触する時刻の差を算出し、そ
れらの頻度分布を求める。
説明する。図１は、計測対象である配管５の中を流れる
液体（水）中の気泡６の速度分布と気泡径分布を測定す
る本発明の装置の一実施例である。装置は、大別して、
電気抵抗式二探針プローブ３，計測器１０，演算器２
０，表示器３０と入力器３１の五つから構成する。電気
抵抗式二探針プローブ３は腐食に強いステンレス鋼や白
金でできた上流側プローブ１，下流側プローブ２からな
り、それぞれ、定圧電源１１ａ，１１ｂと抵抗１２ａ，
１２ｂを有する計測器１０に接続している。気泡は電気
抵抗が大きく、液相（水）は溶解しているイオンにより
電気抵抗が小さいため、プローブが気泡に接したときの
出力端子１３の電位は高く、液相に接したときの出力端
子１３の電位は低くなる。出力端子１３ａ，１３ｂの信
号は、Ａ／Ｄ変換器２１を内蔵した演算器２０に取り込
まれ、データ処理部２２に送られる。データ処理部２２
では、上流側プローブ１と下流側プローブ２の信号か
ら、気泡がプローブに接触している時間と、気泡が上流
側と下流側のプローブに接触する時刻の差を算出し、そ
れらの頻度分布を求める。
【００３２】一方、モンテカルロ演算部２３では、気泡
の速度分布，気泡径分布を仮定して、作用の項で説明し
た手法を利用して、気泡がプローブに接触している時間
を数９と数１０で、気泡が上流側と下流側のプローブに
接触する時刻の差を数１１で算出し、それらの頻度分布
を求める。比較計算部２４では、データ処理部２２とモ
ンテカルロ演算部２３で得た頻度分布を比較し、頻度分
布あるいは比較結果をＣＲＴ，プリンタなどの表示器３
０に表示すると同時に、頻度分布が大きく異なる場合に
は、計算条件変更部２５で仮定した気泡の速度分布，気
泡径分布などの計算条件を変更し、再び、モンテカルロ
演算部２３で頻度分布を計算する。データ処理部２２と
モンテカルロ演算部２３で得た頻度分布が十分一致した
とみなせたときには、その時の気泡の速度分布，気泡径
分布を表示器３０に表示する。キーボードなどの入力器
３１からは、測定の開始信号をデータ処理部２２に送っ
たり、計算条件変更部２５の条件変更に対する制約条件
を送ることができる。また、データ処理部２２から、計
算条件変更部２５の条件変更に対する情報を送ることが
できる。なお、この実施例では、気泡がプローブに接触
している時間と、気泡が上流側と下流側のプローブに接
触する時刻の差の頻度分布を、測定値とシミュレーショ
ンで比較したが、これらの時間から気泡の速度分布と気
泡が探針を横切る長さ（コード長さ）の分布を計算し、
これら二つの頻度分布をモンテカルロ計算と比較しても
よい。
の速度分布，気泡径分布を仮定して、作用の項で説明し
た手法を利用して、気泡がプローブに接触している時間
を数９と数１０で、気泡が上流側と下流側のプローブに
接触する時刻の差を数１１で算出し、それらの頻度分布
を求める。比較計算部２４では、データ処理部２２とモ
ンテカルロ演算部２３で得た頻度分布を比較し、頻度分
布あるいは比較結果をＣＲＴ，プリンタなどの表示器３
０に表示すると同時に、頻度分布が大きく異なる場合に
は、計算条件変更部２５で仮定した気泡の速度分布，気
泡径分布などの計算条件を変更し、再び、モンテカルロ
演算部２３で頻度分布を計算する。データ処理部２２と
モンテカルロ演算部２３で得た頻度分布が十分一致した
とみなせたときには、その時の気泡の速度分布，気泡径
分布を表示器３０に表示する。キーボードなどの入力器
３１からは、測定の開始信号をデータ処理部２２に送っ
たり、計算条件変更部２５の条件変更に対する制約条件
を送ることができる。また、データ処理部２２から、計
算条件変更部２５の条件変更に対する情報を送ることが
できる。なお、この実施例では、気泡がプローブに接触
している時間と、気泡が上流側と下流側のプローブに接
触する時刻の差の頻度分布を、測定値とシミュレーショ
ンで比較したが、これらの時間から気泡の速度分布と気
泡が探針を横切る長さ（コード長さ）の分布を計算し、
これら二つの頻度分布をモンテカルロ計算と比較しても
よい。
【００３３】図４には、本発明によるモンテカルロ演算
部の構成の一実施例を表わす。まず、気泡速度の平均値
と標準偏差，気泡径の平均値と標準偏差といったよう
な、分布を特徴づけるパラメータを設定する。乱数発生
ルーチンを使って、気泡径，気泡速度，流れの方向に対
する気泡の運動方向を表わす角度を次々に設定し、気泡
が上流側プローブに接触している時間ｔtip ，気泡が下
流側プローブに接触している時間ｔrear，気泡が上流側
と下流側のプローブに接触する時刻の差ｔtrを、所望の
個数の気泡について計算する。最後に、ｔtip ，ｔrea
r，ｔtrの頻度分布を求める。
部の構成の一実施例を表わす。まず、気泡速度の平均値
と標準偏差，気泡径の平均値と標準偏差といったよう
な、分布を特徴づけるパラメータを設定する。乱数発生
ルーチンを使って、気泡径，気泡速度，流れの方向に対
する気泡の運動方向を表わす角度を次々に設定し、気泡
が上流側プローブに接触している時間ｔtip ，気泡が下
流側プローブに接触している時間ｔrear，気泡が上流側
と下流側のプローブに接触する時刻の差ｔtrを、所望の
個数の気泡について計算する。最後に、ｔtip ，ｔrea
r，ｔtrの頻度分布を求める。
【００３４】計算条件変更部２５における、仮定した気
泡の速度分布，気泡径分布などの計算条件の変更方法
は、最小二乗法で評価した実測値とシミュレーション結
果の差の標準偏差を評価関数として、この値が小さくな
るようにパラメータを変更する。また、入力器３１か
ら、パラメータの取り得る範囲を指定し、この制約条件
のもとで最適化することもできる。
泡の速度分布，気泡径分布などの計算条件の変更方法
は、最小二乗法で評価した実測値とシミュレーション結
果の差の標準偏差を評価関数として、この値が小さくな
るようにパラメータを変更する。また、入力器３１か
ら、パラメータの取り得る範囲を指定し、この制約条件
のもとで最適化することもできる。
【００３５】図１の実施例では電気抵抗式プローブを用
いているが、二探針の光プローブを用いても、計測器１
０をそれに適したものに置き換えれば、同じ構成で気泡
の速度分布，気泡径分布を測定できる。また、同じ装置
で、環状流や噴霧流中の液滴の速度分布，液滴径分布を
測定できる。
いているが、二探針の光プローブを用いても、計測器１
０をそれに適したものに置き換えれば、同じ構成で気泡
の速度分布，気泡径分布を測定できる。また、同じ装置
で、環状流や噴霧流中の液滴の速度分布，液滴径分布を
測定できる。
【００３６】さらに、プローブの数が三個以上の複数探
針プローブの場合でも、同様の処理をして、気泡の速度
分布，気泡径分布を測定できる。この場合、二探針プロ
ーブに比べて情報量が増加しているので、気泡の形状
が、楕円体になっている影響も評価でき、気泡の寸法を
一層正確に測定できる。
針プローブの場合でも、同様の処理をして、気泡の速度
分布，気泡径分布を測定できる。この場合、二探針プロ
ーブに比べて情報量が増加しているので、気泡の形状
が、楕円体になっている影響も評価でき、気泡の寸法を
一層正確に測定できる。
【００３７】
【発明の効果】本発明によれば、気泡や液滴の実際の速
度分布や直径の分布を精度良く算出することができる。
度分布や直径の分布を精度良く算出することができる。
【図１】本発明の二相流計測装置の第一の実施例を表わ
すブロック図。
すブロック図。
【図２】角度θ，ωの定義を表わす説明図。
【図３】モンテカルロ法による解析モデルを表わす説明
図。
図。
【図４】モンテカルロ演算部の構成の一実施例を表わす
説明図。
説明図。
１…上流側プローブ、２…下流側プローブ、３…電気抵
抗式二探針プローブ、５…計測配管、６…気泡、１０…
計測器、１１ａ，１１ｂ…定圧電源、１２ａ，１２ｂ…
抵抗、１３ａ，１３ｂ…出力端子、２０…演算器、２１
…Ａ／Ｄ変換器、２２…データ処理部、２３…モンテカ
ルロ演算部、２４…比較計算部、２５…計算条件変更
部、３０…表示器、３１…入力器。
抗式二探針プローブ、５…計測配管、６…気泡、１０…
計測器、１１ａ，１１ｂ…定圧電源、１２ａ，１２ｂ…
抵抗、１３ａ，１３ｂ…出力端子、２０…演算器、２１
…Ａ／Ｄ変換器、２２…データ処理部、２３…モンテカ
ルロ演算部、２４…比較計算部、２５…計算条件変更
部、３０…表示器、３１…入力器。
Claims (2)
- 【請求項１】複数の探針を用いて気液二相流の分散相、
すなわち、液相中の気泡、あるいは気相中の液滴の速度
分布を計測する二相流計測装置において、前記分散相が
二個の前記探針に接触する時刻の差の分布と、前記分散
相が探針に接触している時間の分布を測定し、前記分散
相の速度分布と気泡径分布を仮定したモンテカルロ計算
による前記二項目の分布とを比較することにより、前記
分散相の実際の速度分布と気泡径分布を検出することを
特徴とする二相流計測装置。 - 【請求項２】複数の探針を用いて気液二相流の分散相、
すなわち、液相中の気泡、あるいは気相中の液滴の速度
分布を計測する二相流計測装置において、前記分散相が
二個の前記探針に接触する時刻の差から前記分散相の見
かけの速度分布を算出し、前記分散相の見かけの速度分
布と前記分散相が前記探針に接触している時間の分布か
ら、前記探針が前記分散相を横切る長さの分布を測定
し、前記分散相の速度分布と気泡径分布を仮定したモン
テカルロ計算による前記二項目の分布とを比較すること
により、前記分散相の実際の速度分布と気泡径分布を検
出することを特徴とする二相流計測装置。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP16573692A JPH063365A (ja) | 1992-06-24 | 1992-06-24 | 二相流計測装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP16573692A JPH063365A (ja) | 1992-06-24 | 1992-06-24 | 二相流計測装置 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| JPH063365A true JPH063365A (ja) | 1994-01-11 |
Family
ID=15818106
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP16573692A Pending JPH063365A (ja) | 1992-06-24 | 1992-06-24 | 二相流計測装置 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JPH063365A (ja) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2013238490A (ja) * | 2012-05-15 | 2013-11-28 | Kyoto Univ | 気液２相流パラメータ測定装置及びコンピュータプログラム |
| JP2013246124A (ja) * | 2012-05-29 | 2013-12-09 | Central Research Institute Of Electric Power Industry | 三次元速度の計測システム |
-
1992
- 1992-06-24 JP JP16573692A patent/JPH063365A/ja active Pending
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2013238490A (ja) * | 2012-05-15 | 2013-11-28 | Kyoto Univ | 気液２相流パラメータ測定装置及びコンピュータプログラム |
| JP2013246124A (ja) * | 2012-05-29 | 2013-12-09 | Central Research Institute Of Electric Power Industry | 三次元速度の計測システム |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Stutz et al. | Measurements within unsteady cavitation | |
| Hibiki et al. | Local measurement of interfacial area, interfacial velocity and liquid turbulence in two-phase flow | |
| US20140200836A1 (en) | Flow rate determination method and apparatus | |
| McClure et al. | Experimental investigation into the drag volume fraction correction term for gas-liquid bubbly flows | |
| Rodrigues et al. | Statistical features of the flow evolution in horizontal liquid-gas slug flow | |
| Pietri et al. | Simultaneous measurements of temperature and velocity fluctuations in a slightly heated jet combining a cold wire and Laser Doppler Anemometry | |
| JP2012132797A (ja) | 超音波流量計の検証システム及び検証方法 | |
| CN113188613B (zh) | 一种基于不确定度分析的多相流量测量方法和系统 | |
| Hayakawa et al. | Turbulence measurements in a compressible reattaching shear layer | |
| WO2019045703A1 (en) | DETECTION AND IDENTIFICATION OF CHANGE IN A VIBRATORY MEASURING DEVICE | |
| Olbrich et al. | Statistical characterization of horizontal slug flow using snapshot proper orthogonal decomposition | |
| Geršl et al. | Flow rate measurement in stacks with cyclonic flow–Error estimations using CFD modelling | |
| Dayev et al. | Modeling of coefficient of contraction of differential pressure flowmeters | |
| Miller et al. | Evaluation of hot-wire spatial filtering corrections for wall turbulence and correction for end-conduction effects | |
| Mishra et al. | A model for obtaining the velocity vectors of spherical droplets in multiphase flows from measurements using an orthogonal four-sensor probe | |
| JPH063365A (ja) | 二相流計測装置 | |
| JP3586671B2 (ja) | 大気汚染物質の濃度予測方法及び装置 | |
| Suzanne et al. | Local measurements in turbulent bubbly flows | |
| Gurau et al. | Measurement of gas and liquid velocities in an air–water two-phase flow using cross-correlation of signals from a double sensor hot-film probe | |
| Freire et al. | Bubble plumes and the Coanda effect | |
| US9488509B2 (en) | Method or determining an absolute flow rate of a volume or mass flow | |
| JP2003228646A (ja) | 河川状態シミュレーション方法 | |
| Loureiro et al. | Water tank and numerical model studies of flow over steep smooth two-dimensional hills | |
| Poole et al. | Turbulent flow of a viscoelastic shear-thinning liquid through a plane sudden expansion of modest aspect ratio | |
| Lapin et al. | Investigation of conversion function for vortex sonic flowmeter using monte-carlo method |