Top > Search of Japanese Patents > ANALYSIS METHOD OF BUBBLE FLOW, BUBBLE FLOW SIMULATION PROGRAM, STORAGE MEDIUM WITH IT STORED AND BUBBLE FLOW SIMULATION DEVICE

ANALYSIS METHOD OF BUBBLE FLOW, BUBBLE FLOW SIMULATION PROGRAM, STORAGE MEDIUM WITH IT STORED AND BUBBLE FLOW SIMULATION DEVICE commons

Patent code P06P003530
File No. NU-0018
Posted date Apr 28, 2006
Application number P2004-272011
Publication number P2006-085603A
Patent number P4032123
Date of filing Sep 17, 2004
Date of publication of application Mar 30, 2006
Date of registration Nov 2, 2007
Inventor
  • (In Japanese)内山 知実
  • (In Japanese)出川 智啓
Applicant
  • (In Japanese)学校法人名古屋大学
Title ANALYSIS METHOD OF BUBBLE FLOW, BUBBLE FLOW SIMULATION PROGRAM, STORAGE MEDIUM WITH IT STORED AND BUBBLE FLOW SIMULATION DEVICE commons
Abstract

PROBLEM TO BE SOLVED: To provide an analysis method of a bubble flow capable of accurately analyzing the behavior of the bubble flow where a gas and a liquid are mixed; to provide a bubble flow simulation program, a storage medium with it stored, and a bubble flow simulation device.

SOLUTION: In this analysis method of a bubble flow, the bubble flow is discretized with a plurality of vortex elements 16 and 17 by focusing attention on a vorticity field of the bubble flow where bubbles G and a liquid L are mixed; a vorticity transportation equation is formulated from a momentum conservation relationship and a mass conservation relationship considering a gas-liquid volume ratio (a bubble volume fraction αg, a liquid volume fraction αl) between the bubbles G and the liquid L in the bubble flow; and the behavior of the discrete vortex element is obtained from Lagrange calculation of the vorticity transportation equation.

Outline of related art and contending technology (In Japanese)


微小な気泡と液体とが混在して相互作用を及ぼし合いながら流れる気泡流は、例えば反応装置、熱交換器及び廃水処理装置など、様々な工業装置において観察される。このため、この気泡流の挙動について、合理的な解析方法の確立が望まれている。



従来の解析方法では、流れの変数として速度、圧力及び気泡体積率が用いられ、それらを支配する複数の偏微分方程式を様々な数値解法(差分法、有限体積法、有限要素法)によって計算するようにしていた。これらの偏微分方程式は、微分を差分で近似表現することにより、代数方程式に置き換えられるが、非線形項が数値不安定性を有するため、気泡流のレイノルズ数が高い場合には計算過程で発散が生じるおそれがある。特に、工業装置で見受けられる気泡流は、レイノルズ数が高いものが多く、精度の高い気泡流解析ができないという問題があった。



ところで、従来から、液体又は気体のみの単相流の解析については渦法が有効に利用されている。ここで、渦法とは、渦度をもつ渦要素を追跡して渦度場の時間変化を求めるLagrange 型解法であり,大規模渦の形成や変形など渦構造の発展過程を良好に計算できるものである。そして、下記非特許文献1には、この渦法を利用して、気泡流の挙動を解析する技術が開示されている。

【非特許文献1】Sene, K. J., et al., The role of coherent structures in bubble transport by turbulent shear flows,J. Fluid Mech., 259(1994), 219-240.

【特許文献1】特開2004-126925公報

【特許文献2】特開2004-118394公報

【特許文献3】特開2001-99748公報

【特許文献4】特開平9-329109号公報

【特許文献5】特開2003-331208公報

【特許文献6】特開2004-503005公報

【特許文献7】特開2001-132700公報

Field of industrial application (In Japanese)


本発明は、液中に分散した気泡が液体と相互作用を及ぼし合いながら流れる気泡流(二相流)の挙動の解析方法に関する。

Scope of claims (In Japanese)
【請求項1】
 
気泡と液体とが混在した気泡流の解析するための気泡流シミュレーションプログラムであって、コンピュータに、次のステップ(a)~(h)を含むステップを繰り返し実行させて、解析領域内に想定した評価点について求まる前記液体のスカラーポテンシャル及びベクトルポテンシャルの時間的変化に基づいて、時間的に変化する前記気泡流の挙動を解析する気泡流シミュレーションプログラム。
前記コンピュータは、前記ステップ(c)、(d)、(e)、(f)用の記憶手段を備え、前記気泡流シミュレーションプログラムは、
(a)前記解析領域内に複数の離散渦要素を順次導入し、そのときに前記解析領域内に存在する各離散渦要素の各位置データを算出するステップ。
(b)前記解析領域内に複数の気泡を導入し、そのときに前記解析領域内に存在する各気泡の位置を算出するステップ。
(c)前記解析領域を複数に区画した各解析要素内における気泡と液体との気液体積割合を、前記ステップ(b)で算出された各気泡の位置に基づき求めて記憶手段に保存するステップ。
(d)前記ステップ(c)によって保存された気液体積割合データを読み出して、前記各解析要素内における循環強度変化を、当該気液体積割合の気泡流についての運動量保存関係から導出された気泡流の渦度輸送方程式と、レイノルズの輸送定理とに基づき求めて記憶手段に保存するステップ。
前記気泡流の渦度輸送方程式は、次の数式[数1]である。
【数式1】
 



Γ :解析要素内における循環強度変化
αl:気泡と液体との気液体積割合
g :重力加速度
ul:液体の速度
A :解析要素の面積
(e)前記ステップ(d)によって保存された前記循環強度変化データを読み出して、前記各解析要素内に存在する各離散渦要素の循環強度及び位置を、当該循環強度変化と、前記ステップ(a)で算出された各離散渦要素の位置に基づきラグランジェ的に求めて記憶手段に保存するステップ。
(f)前記ステップ(e)によって保存された各離散渦要素の循環強度及び位置データを読み出して、前記解析領域内の前記評価点における渦度を、当該各離散渦要素の循環強度及び位置に基づき求めて記憶手段に保存するステップ。
(g)前記ステップ(f)によって保存された前記評価点における渦度データを読み出して当該渦度から前記液体のベクトルポテンシャルを求め、前記評価点における前記液体のスカラーポテンシャルを、当該液体のベクトルポテンシャルと、前記体積割合の気泡流の質量保存関係及びヘルムホルツの定理によって導出された液体のスカラーポテンシャル式とに基づき求めるステップ。
前記液体のスカラーポテンシャル式は、次の数式[数2]である。
【数式2】
 



αl:気泡と液体との気液体積割合
φ :液体のスカラーポテンシャル
ψ :液体のベクトルポテンシャル
(h) 前記ステップ(g)によって求められた前記評価点における前記液体の速度ポテンシャルおよびスカラーポテンシャルを用いて、前記評価点における前記液体の速度を計算するステップ。

【請求項2】
 
前記請求項1に記載の気泡流シミュレーションプログラムを記憶した記憶媒体。

【請求項3】
 
気泡と液体とが混在した気泡流の解析するための気泡流シミュレーション装置であって、次の手段(a)~(h)を備えるとともに、
前記手段(a)で算出される前記離散渦要素の位置データが保存される渦要素情報記憶手段と、 前記手段(b)で算出される前記気泡の位置データが保存される気泡情報記憶手段と、前記手段(c)、(d)、(e)、(f)用の記憶手段とを備えており、
前記手段(a)~(h)を繰り返し実行させて、解析領域内に想定した評価点について求まる前記液体のスカラーポテンシャル及びベクトルポテンシャルの時間的変化に基づいて、時間的に変化する前記気泡流の挙動を解析する気泡流シミュレーション装置。
(a)前記解析領域内に複数の離散渦要素を順次導入し、そのときに前記解析領域内に存在する各離散渦要素の各位置データを算出する手段。
(b)前記解析領域内に複数の気泡を導入し、そのときに前記解析領域内に存在する各気泡の位置を算出する手段。
(c)前記解析領域を複数に区画した各解析要素内における気泡と液体との気液体積割合を、前記手段(b)で算出された各気泡の位置に基づき求めて記憶手段に保存する手段。
(d)前記手段(c)によって保存された気液体積割合データを読み出して、前記各解析要素内における循環強度変化を、当該気液体積割合の気泡流についての運動量保存関係から導出された気泡流の渦度輸送方程式と、レイノルズの輸送定理とに基づき求めて記憶手段に保存する手段。
前記気泡流の渦度輸送方程式は、次の数式[数1]である。
【数式1】
 



Γ :解析要素内における循環強度変化
αl:気泡と液体との気液体積割合
g :重力加速度
ul:液体の速度
A :解析要素の面積

(e)前記手段(d)によって保存された前記循環強度変化データを読み出して、前記各解析要素内に存在する各離散渦要素の循環強度及び位置を、当該循環強度変化と、前記手段(a)で算出された各離散渦要素の位置に基づきラグランジェ的に求めて記憶手段に保存する手段。
(f)前記手段(e)によって保存された各離散渦要素の循環強度及び位置データを読み出して、前記解析領域内の前記評価点における渦度を、当該各離散渦要素の循環強度及び位置に基づき求めて記憶手段に保存する手段。
(g)前記手段(f)によって保存された前記評価点における渦度データを読み出して当該渦度から前記液体のベクトルポテンシャルを求め、前記評価点における前記液体のスカラーポテンシャルを、当該液体のベクトルポテンシャルと、前記体積割合の気泡流の質量保存関係及びヘルムホルツの定理によって導出された液体のスカラーポテンシャル式とに基づき求める手段。
液体のスカラーポテンシャル式は、次の数式[数2]である。
【数式2】
 



αl:気泡と液体との気液体積割合
φ :液体のスカラーポテンシャル
ψ :液体のベクトルポテンシャル
(h) 前記手段(g)によって求められた前記評価点における前記液体の速度ポテンシャルおよびスカラーポテンシャルを用いて、前記評価点における前記液体の速度を計算する手段。

【請求項4】
 
前記手段(a)で算出される前記離散渦要素の位置データが保存される渦要素情報記憶手段と、
前記手段(b)で算出される前記気泡の位置データ及び前記手段(c)で算出される気液体積割合データが保存される気泡情報記憶手段とを備え、
前記手段(a)による処理と、前記手段(b)及び前記手段(c)による処理とを並列的に実行することを特徴とする請求項3に記載の気泡流シミュレーション装置。
Industrial division
  • Application of computer
IPC(International Patent Classification)
Drawing

※Click image to enlarge.

JP2004272011thum.jpg
State of application right Right is in force
(In Japanese)名古屋大学の公開特許情報を掲載しています。ご関心のある案件がございましたら、下記まで電子メールでご連絡ください。


PAGE TOP

close
close
close
close
close
close
close