(57)【要約】【課題】液滴の物性や初期分布から液滴の粒径や空間分布の時間変化を精度良く推定することができる液滴成長シミュレーション装置、液滴成長シミュレーション方法、液滴成長シミュレーションプログラムおよび、液滴成長シミュレーションプログラムを記録したコンピュータ読み取り可能な記録媒体を提供する。【解決手段】入力手段１１により、液滴の粒径、物性、分布と、液滴に作用する外力と、流体の物性とを含む初期情報を入力し、液滴計算手段１２により、所定の経過時間後の液滴の粒径と分布とを計算し、出力手段１３により、計算された液滴の粒径と分布とを出力する。液滴計算手段１２は、運動方程式に基づいて液滴の分布を計算する運動計算手段１２ａと、その液滴の分布から液滴が合体成長する確率を求め、液滴の合体処理を行って液滴の粒径を求める液滴合体手段１２ｂとを有している。【選択図】図２
2014049054.tif 000002
【特許請求の範囲】
【請求項１】
流体中の液滴の粒径や分布の時間変化を推定するための液滴成長シミュレーション装置であって、
前記液滴の粒径、物性、分布と、前記液滴に作用する外力と、前記流体の物性とを含む初期情報を入力する入力手段と、
前記入力手段で入力された前記初期情報に基づいて、所定の経過時間後の前記液滴の粒径と分布とを計算する液滴計算手段と、
前記液滴計算手段で計算された前記液滴の粒径と分布とを出力する出力手段とを有し、
前記液滴計算手段は、下記の（１）式の運動方程式に基づいて、前記液滴の分布を計算する運動計算手段と、
【数１】
2014049054.tif 000009
前記運動計算手段で計算された前記液滴の分布から、下記の（２）式の表面自由エネルギーＧに基づいて、前記液滴が合体成長する確率を求め、その確率に基づいて、前記液滴の合体処理を行って前記液滴の粒径を求める液滴合体手段とを有し、
【数２】
2014049054.tif 000010
前記所定の経過時間になるまで、所定の時間間隔で前記運動計算手段と前記液滴合体手段とを繰り返すことを
特徴とする液滴成長シミュレーション装置。
【請求項２】
前記初期情報は、前記液滴に作用する外力が重力のみであり、前記流体の流速がゼロであることを特徴とする請求項１記載の液滴成長シミュレーション装置。
【請求項３】
前記運動計算手段は、前記液滴間に働く相互作用として、ファンデアワールス相互作用と静電相互作用とを考慮することを特徴とする請求項１または２記載の液滴成長シミュレーション装置。
【請求項４】
前記液滴合体手段は、前記（２）式の表面自由エネルギーＧを、下記の（３）式に基づいて、単純な表面積の変化によるエネルギー変化のうち、減少する表面積分の項Ｇ１、増加する表面積分の項Ｇ２、曲がっていることによるエネルギー変化Ｇ３、ガウス曲率Ｇｇに分けて算出することを特徴とする請求項１乃至３のいずれか１項に記載の液滴成長シミュレーション装置。
【数３】
2014049054.tif 000011
【請求項５】
前記液滴合体手段は、自発曲率Ｈ０および表面エネルギーγが、下記の（４）式および（５）式で表されることを特徴とする請求項４記載の液滴成長シミュレーション装置。
【数４】
2014049054.tif 000012
【請求項６】
流体中の液滴の粒径や分布の時間変化を推定するための液滴成長シミュレーション方法であって、
前記液滴の粒径、物性、分布と、前記液滴に作用する外力と、前記流体の物性とを含む初期情報を入力する入力工程と、
前記入力工程で入力された前記初期情報に基づいて、所定の経過時間後の前記液滴の粒径と分布とを計算する液滴計算工程と、
前記液滴計算工程で計算された前記液滴の粒径と分布とを出力する出力工程とを有し、
前記液滴計算工程は、下記の（１）式の運動方程式に基づいて、前記液滴の分布を計算する運動計算工程と、
【数１】
2014049054.tif 000013
前記運動計算工程で計算された前記液滴の分布から、下記の（２）式の表面自由エネルギーＧに基づいて、前記液滴が合体成長する確率を求め、その確率に基づいて、前記液滴の合体処理を行って前記液滴の粒径を求める液滴合体工程とを有し、
【数２】
2014049054.tif 000014
前記所定の経過時間になるまで、所定の時間間隔で前記運動計算工程と前記液滴合体工程とを繰り返すことを
特徴とする液滴成長シミュレーション方法。
【請求項７】
流体中の液滴の粒径や分布の時間変化を推定するための液滴成長シミュレーションプログラムであって、
コンピュータを、
前記液滴の粒径、物性、分布と、前記液滴に作用する外力と、前記流体の物性とを含む初期情報を入力する入力手段、
前記入力手段で入力された前記初期情報に基づいて、所定の経過時間後の前記液滴の粒径と分布とを計算する液滴計算手段、
前記液滴計算手段で計算された前記液滴の粒径と分布とを出力する出力手段として機能させ、
前記液滴計算手段は、下記の（１）式の運動方程式に基づいて、前記液滴の分布を計算する運動計算手段と、
【数１】
2014049054.tif 000015
前記運動計算手段で計算された前記液滴の分布から、下記の（２）式の表面自由エネルギーＧに基づいて、前記液滴が合体成長する確率を求め、その確率に基づいて、前記液滴の合体処理を行って前記液滴の粒径を求める液滴合体手段とを有し、
【数２】
2014049054.tif 000016
前記所定の経過時間になるまで、所定の時間間隔で前記運動計算手段と前記液滴合体手段とを繰り返すことを
特徴とする液滴成長シミュレーションプログラム。
【請求項８】
流体中の液滴の粒径や分布の時間変化を推定するための液滴成長シミュレーションプログラムを記録したコンピュータ読み取り可能な記録媒体であって、
コンピュータを、
前記液滴の粒径、物性、分布と、前記液滴に作用する外力と、前記流体の物性とを含む初期情報を入力する入力手段、
前記入力手段で入力された前記初期情報に基づいて、所定の経過時間後の前記液滴の粒径と分布とを計算する液滴計算手段、
前記液滴計算手段で計算された前記液滴の粒径と分布とを出力する出力手段として機能させ、
前記液滴計算手段は、下記の（１）式の運動方程式に基づいて、前記液滴の分布を計算する運動計算手段と、
【数１】
2014049054.tif 000017
前記運動計算手段で計算された前記液滴の分布から、下記の（２）式の表面自由エネルギーＧに基づいて、前記液滴が合体成長する確率を求め、その確率に基づいて、前記液滴の合体処理を行って前記液滴の粒径を求める液滴合体手段とを有し、
【数２】
2014049054.tif 000018
前記所定の経過時間になるまで、所定の時間間隔で前記運動計算手段と前記液滴合体手段とを繰り返すことを
特徴とする液滴成長シミュレーションプログラムを記録したコンピュータ読み取り可能な記録媒体。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、流体中の液滴の粒径や分布の時間変化を推定するための液滴成長シミュレーション装置、液滴成長シミュレーション方法、液滴成長シミュレーションプログラムおよび、液滴成長シミュレーションプログラムを記録したコンピュータ読み取り可能な記録媒体に関する。
【背景技術】
【０００２】
各種材料の溶液を用いた製造プロセスの設計を行うにあたり、各種液滴に関する粒成長すなわち液滴の粒径や分布の時間変化を推定する必要がある。従来、液滴に関する時間変化を推定するための液滴成長シミュレーションとして、自動車のフロントガラス上の水滴シミュレーションを行って水滴マップを生成することにより、雨滴による視界の悪化をリアルに表現することができるシミュレーション装置（例えば、特許文献１参照）や、燃料電池のガスチャンネルを流れる反応ガスの流動解析を行うにあたり、反応ガスを気液二相流とみなして水滴の挙動を解析し、数値演算が発散することなく確実に収束解を得ることができるシミュレーション方法（例えば、特許文献２参照）、燃料電池のガスチャンネルを流れる反応ガスの流動状態がミスト流またはやや大きめの液滴であることに着目し、水滴をすべて球体として取り扱ってラグランジュ的なシミュレーション解析を行うことにより、数値演算が発散することなく確実に収束解を得ることができるシミュレーション方法（例えば、特許文献３参照）がある。
【０００３】
なお、液滴の合体を検討する方法として、界面曲率を考慮したモデルを構築し、そのモデルに基づいて液滴の合体成長を検討する方法が提供されている（例えば、非特許文献１または２参照）。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特開２００９−０７０２１５号公報
【特許文献２】特開２００５−１３５８１３号公報
【特許文献３】特開２００５−１３５８１４号公報
【非特許文献】
【０００５】
【非特許文献１】A. Kabalnov, H. Wennerstrom, “Macroemulsion Stability: TheOriented Wedge Theory Revisited”, Langmuir, 1996, 12, p.276-292
【非特許文献２】H. Wennerstrom, U. Olsson, “Microemulations as modelsystems”, C. R. Chimie, 2009, 12, p.4-17
【発明の概要】
【発明が解決しようとする課題】
【０００６】
特許文献１に記載のシミュレーション装置は、ドライブシミュレータにおける画像処理に関する技術であり、特許文献２および３に記載のシミュレーション方法は、燃料電池シミュレータにおけるガスチャネルでの気液二相流の収束計算に関するモデリング技術に関する技術である。このように、これらの従来の装置や方法は、液滴の物性や初期分布から液滴の粒径や空間分布の時間変化を推定するものではなく、液滴の粒径や空間分布の時間変化を推定する技術の開発が望まれていた。
【０００７】
本発明は、このような課題に着目してなされたもので、液滴の物性や初期分布から液滴の粒径や空間分布の時間変化を精度良く推定することができる液滴成長シミュレーション装置、液滴成長シミュレーション方法、液滴成長シミュレーションプログラムおよび、液滴成長シミュレーションプログラムを記録したコンピュータ読み取り可能な記録媒体を提供することを目的としている。
【課題を解決するための手段】
【０００８】
本発明に係る液滴成長シミュレーション装置は、流体中の液滴の粒径や分布の時間変化を推定するための液滴成長シミュレーション装置であって、前記液滴の粒径、物性、分布と、前記液滴に作用する外力と、前記流体の物性とを含む初期情報を入力する入力手段と、前記入力手段で入力された前記初期情報に基づいて、所定の経過時間後の前記液滴の粒径と分布とを計算する液滴計算手段と、前記液滴計算手段で計算された前記液滴の粒径と分布とを出力する出力手段とを有し、前記液滴計算手段は、下記の（１）式の運動方程式に基づいて、前記液滴の分布を計算する運動計算手段と、
【数１】
2014049054.tif 000003
前記運動計算手段で計算された前記液滴の分布から、下記の（２）式の表面自由エネルギーＧに基づいて、前記液滴が合体成長する確率を求め、その確率に基づいて、前記液滴の合体処理を行って前記液滴の粒径を求める液滴合体手段とを有し、
【数２】
2014049054.tif 000004
前記所定の経過時間になるまで、所定の時間間隔で前記運動計算手段と前記液滴合体手段とを繰り返すことを特徴とする。
【０００９】
本発明に係る液滴成長シミュレーション方法は、流体中の液滴の粒径や分布の時間変化を推定するための液滴成長シミュレーション方法であって、前記液滴の粒径、物性、分布と、前記液滴に作用する外力と、前記流体の物性とを含む初期情報を入力する入力工程と、前記入力工程で入力された前記初期情報に基づいて、所定の経過時間後の前記液滴の粒径と分布とを計算する液滴計算工程と、前記液滴計算工程で計算された前記液滴の粒径と分布とを出力する出力工程とを有し、前記液滴計算工程は、（１）式の運動方程式に基づいて、前記液滴の分布を計算する運動計算工程と、前記運動計算工程で計算された前記液滴の分布から、（２）式の表面自由エネルギーＧに基づいて、前記液滴が合体成長する確率を求め、その確率に基づいて、前記液滴の合体処理を行って前記液滴の粒径を求める液滴合体工程とを有し、前記所定の経過時間になるまで、所定の時間間隔で前記運動計算工程と前記液滴合体工程とを繰り返すことを特徴とする。
【００１０】
本発明に係る液滴成長シミュレーションプログラムは、流体中の液滴の粒径や分布の時間変化を推定するための液滴成長シミュレーションプログラムであって、コンピュータを、前記液滴の粒径、物性、分布と、前記液滴に作用する外力と、前記流体の物性とを含む初期情報を入力する入力手段、前記入力手段で入力された前記初期情報に基づいて、所定の経過時間後の前記液滴の粒径と分布とを計算する液滴計算手段、前記液滴計算手段で計算された前記液滴の粒径と分布とを出力する出力手段として機能させ、前記液滴計算手段は、（１）式の運動方程式に基づいて、前記液滴の分布を計算する運動計算手段と、前記運動計算手段で計算された前記液滴の分布から、（２）式の表面自由エネルギーＧに基づいて、前記液滴が合体成長する確率を求め、その確率に基づいて、前記液滴の合体処理を行って前記液滴の粒径を求める液滴合体手段とを有し、前記所定の経過時間になるまで、所定の時間間隔で前記運動計算手段と前記液滴合体手段とを繰り返すことを特徴とする。
【００１１】
本発明に係るコンピュータ読み取り可能な記録媒体は、流体中の液滴の粒径や分布の時間変化を推定するための液滴成長シミュレーションプログラムを記録したコンピュータ読み取り可能な記録媒体であって、コンピュータを、前記液滴の粒径、物性、分布と、前記液滴に作用する外力と、前記流体の物性とを含む初期情報を入力する入力手段、前記入力手段で入力された前記初期情報に基づいて、所定の経過時間後の前記液滴の粒径と分布とを計算する液滴計算手段、前記液滴計算手段で計算された前記液滴の粒径と分布とを出力する出力手段として機能させ、前記液滴計算手段は、（１）式の運動方程式に基づいて、前記液滴の分布を計算する運動計算手段と、前記運動計算手段で計算された前記液滴の分布から、（２）式の表面自由エネルギーＧに基づいて、前記液滴が合体成長する確率を求め、その確率に基づいて、前記液滴の合体処理を行って前記液滴の粒径を求める液滴合体手段とを有し、前記所定の経過時間になるまで、所定の時間間隔で前記運動計算手段と前記液滴合体手段とを繰り返すことを特徴とする。
【００１２】
本発明に係る液滴成長シミュレーション装置、液滴成長シミュレーション方法、液滴成長シミュレーションプログラムおよび、液滴成長シミュレーションプログラムを記録したコンピュータ読み取り可能な記録媒体は、ブラウン運動を考慮した（１）式のＳｔｏｋｅｓの抵抗法則を用いて液滴の運動を算出して液滴の分布を計算するだけでなく、（２）式のＨｏｏｋｅａｎの方程式を用いて液滴が合体成長する確率を求めることにより、液滴の物性や初期分布から、液滴の移動および合体による粒成長の挙動を考慮して、液滴の粒径や空間分布の時間変化を精度良く推定することができる。
【００１３】
なお、（２）式のＨｏｏｋｅａｎの方程式は、非特許文献１および２に記載のように、界面曲率を考慮したモデルに基づいて液滴の合体成長を検討することにより求められた式であり、本来、運動を伴う液滴を対象とするものではない。しかし、本発明者等は、鋭意努力の結果、（２）式を（１）式の運動方程式と組み合わせることにより、液滴の移動および合体による粒成長の挙動を取り入れたシミュレーションを構築することに初めて成功し、本発明に想到した。
【００１４】
本発明に係る液滴成長シミュレーション装置、液滴成長シミュレーション方法、液滴成長シミュレーションプログラムおよび、液滴成長シミュレーションプログラムを記録したコンピュータ読み取り可能な記録媒体で、前記初期情報は、前記液滴に作用する外力が重力のみであり、前記流体の流速がゼロであってもよい。この場合、重力以外の外力が作用せず、流体の流速がゼロの条件下でのシミュレーションを高速で行うことができる。
【００１５】
本発明に係る液滴成長シミュレーション装置、液滴成長シミュレーション方法、液滴成長シミュレーションプログラムおよび、液滴成長シミュレーションプログラムを記録したコンピュータ読み取り可能な記録媒体で、前記運動計算手段または前記運動計算工程は、前記液滴間に働く相互作用として、ファンデアワールス相互作用と静電相互作用とを考慮することが好ましい。この場合、（１）式の運動方程式での計算過程で、ファンデアワールス相互作用および静電相互作用以外の相互作用を考慮しないため、高速でシミュレーションを行うことができる。また、ファンデアワールス相互作用のみ、または静電相互作用のみを考慮した場合に比べて、高精度でシミュレーションを行うことができる。比較的影響が大きいファンデアワールス相互作用および静電相互作用のみを考慮するため、高速で高精度なシミュレーションを行うことができる。
【００１６】
本発明に係る液滴成長シミュレーション装置、液滴成長シミュレーション方法、液滴成長シミュレーションプログラムおよび、液滴成長シミュレーションプログラムを記録したコンピュータ読み取り可能な記録媒体で、前記液滴合体手段または前記液低合体工程は、前記（２）式の表面自由エネルギーＧを、下記の（３）式に基づいて、単純な表面積の変化によるエネルギー変化のうち、減少する表面積分の項Ｇ１、増加する表面積分の項Ｇ２、曲がっていることによるエネルギー変化Ｇ３、ガウス曲率Ｇｇに分けて算出することが好ましい。
【数３】
2014049054.tif 000005
【００１７】
（３）式は非特許文献１および２で得られた式であり、これを使用する場合、精度良く液滴の合体成長の確率を計算することができる。また、この場合、高速でシミュレーションを行うために、ガウス曲率Ｇｇがガウス・ボンネットの定理により曲面の位相的性質のみに依存し、オイラー標数が２であることが好ましい。
【００１８】
さらに、（３）式を使用する場合、前記液滴合体手段または前記液低合体工程は、自発曲率Ｈ０および表面エネルギーγが、下記の（４）式および（５）式で表されることが好ましい。（４）式および（５）式により、高精度かつ高速でシミュレーションを行うことができる。
【数４】
2014049054.tif 000006
【発明の効果】
【００１９】
本発明によれば、液滴の物性や初期分布から液滴の粒径や空間分布の時間変化を精度良く推定することができる液滴成長シミュレーション装置、液滴成長シミュレーション方法、液滴成長シミュレーションプログラムおよび、液滴成長シミュレーションプログラムを記録したコンピュータ読み取り可能な記録媒体を提供することができる。
【図面の簡単な説明】
【００２０】
【図１】本発明の実施の形態の液滴成長シミュレーション装置を示すブロック図である。
【図２】本発明の実施の形態の液滴成長シミュレーション方法のフローチャートである。
【図３】本発明の実施の形態の液滴成長シミュレーション装置および液滴成長シミュレーション方法による、油滴の濃度が３６０ｍＭおよび４８０ｍＭの２通りの場合についてのシミュレーション結果の、経過時刻と油滴の大きさの平均値との関係を示すグラフである。
【図４】本発明の実施の形態の液滴成長シミュレーション装置および液滴成長シミュレーション方法による、重力加速度が（ａ）１μＧおよび（ｂ）１Ｇの２通りの場合についてのシミュレーション結果の、経過時刻毎の油滴の分布状態を示す側面図である。
【発明を実施するための形態】
【００２１】
以下、図面に基づき、本発明の実施の形態について説明する。
図１乃至図４は、本発明の実施の形態の液滴成長シミュレーション装置、液滴成長シミュレーション方法、液滴成長シミュレーションプログラムおよび、液滴成長シミュレーションプログラムを記録したコンピュータ読み取り可能な記録媒体を示している。
【００２２】
本発明の実施の形態の液滴成長シミュレーション装置１０は、流体中の液滴の粒径や分布の時間変化（時間発展）を推定するための液滴成長シミュレーション装置１０であって、コンピュータから成っている。図１に示すように、本発明の実施の形態の液滴成長シミュレーション装置１０は、入力手段１１と液滴計算手段１２と出力手段１３とを有している。
【００２３】
なお、コンピュータとは、ハードウェアとＯＳ（オペレーティングシステム）とを含む概念であり、ＯＳの制御の下で動作するハードウェアを意味している。また、ＯＳが不要で、アプリケーションプログラム単独でハードウェアを動作させるような場合には、そのハードウェア自体がコンピュータに相当する。ハードウェアは、少なくとも、ＣＰＵ等のマイクロプロセッサと、記録媒体に記録されたコンピュータプログラムを読み取るための手段とを備えている。
【００２４】
入力手段１１は、キーボードやマウスなどから成り、初期状態での液滴の数、粒径、密度や動粘度などの物性、分布（座標）、液滴に作用する重力加速度および重力の方向、流体の物性、温度、シミュレーションにおける計算のステップ数、積分時間（経過時間）、セルサイズなどの初期情報を入力するよう構成されている。入力手段１１は、コンピュータに備え付けられたモニタに写し出された入力画面で初期情報を直接入力するよう構成されていてもよく、入力ファイルを用いて初期情報を一括入力するよう構成されていてもよい。
【００２５】
液滴計算手段１２は、コンピュータに内蔵されたＣＰＵやメモリ等から成り、入力手段１１で入力された初期情報に基づいて、所定の経過時間後の液滴の粒径と分布とを計算するよう構成されている。液滴計算手段１２は、液滴の分布を計算する運動計算手段１２ａと、液滴の粒径を求める液滴合体手段１２ｂとを有している。
【００２６】
運動計算手段１２ａは、（１）式の運動方程式に基づいて、液滴の分布を計算するよう構成されている。運動計算手段１２ａは、入力手段１１で入力された初期情報のうちの各液滴の位置情報や、所定の時刻での各液滴の位置情報から、（１）式に基づいて時刻ステップΔｔ経過後の各液滴の位置情報を算出するようになっている。具体的には、運動計算手段１２ａは、液滴に働く外力を重力のみとし、流体の抵抗が十分大きく流速を無視することができる（流速ゼロ）と仮定して、（１）式から（６）式を求め、時刻ステップΔｔ経過後の液滴の座標を（６）式に従って求めるよう構成されている。
【００２７】
【数５】
2014049054.tif 000007
【００２８】
また、運動計算手段１２ａは、液滴間に働く相互作用として、分子シミュレーションの分野などで用いられるポテンシャルを考慮可能に構成されている。具体的には、運動計算手段１２ａは、ＤＬＶＯ理論（Ｄｅｒｊａｇｕｉｎ，Ｌａｎｄａｕ，Ｖｅｒｗｅｙ，Ｏｖｅｒｂｅｅｋ  ｔｈｅｏｒｙ）に基づいて、相互作用の全エネルギーが、ファンデアワールス相互作用によるエネルギーと静電相互作用によるエネルギーとの和で与えられるものとして計算するよう構成されている。運動計算手段１２ａは、ファンデアワールス相互作用によるエネルギーＶＶＤＷ（Ｈ）を、（７）式に示す方程式に従って計算し、静電相互作用によるエネルギーＶｅｌを、（８）式に示す方程式に従って計算するようになっている。
【００２９】
【数６】
2014049054.tif 000008
【００３０】
液滴合体手段１２ｂは、運動計算手段１２ａで計算された液滴の分布から、メトロポリス法を利用して、（２）式の表面自由エネルギーＧから液滴が合体成長する確率を求め、その確率に基づいて、液滴の合体処理を行って液滴の粒径を求めるよう構成されている。このとき、液滴合体手段１２ｂは、（２）式の表面自由エネルギーＧを、（３）式に基づいて、単純な表面積の変化によるエネルギー変化のうち、減少する表面積分の項Ｇ１、増加する表面積分の項Ｇ２、曲がっていることによるエネルギー変化Ｇ３、ガウス曲率Ｇｇに分けて算出するようになっている。さらに、液滴合体手段１２ｂは、（３）式中のガウス曲率Ｇｇがガウス・ボンネットの定理により曲面の位相的性質のみに依存し、オイラー標数が２であると仮定して計算するようになっている。また、（３）式中の自発曲率Ｈ０が温度と比例関係にあり、表面エネルギーγが温度の自乗と比例関係にあると仮定して、それぞれを（４）式および（５）式で計算するようになっている。
【００３１】
液滴計算手段１２は、所定の経過時間になるまで、すなわち入力手段１１で入力された初期情報のうち、計算のステップ数の数だけ、運動計算手段１２ａと液滴合体手段１２ｂとを繰り返すよう構成されている。
出力手段１３は、モニタやプリンタなどから成り、液滴計算手段１２で計算された液滴の数と粒径と分布（座標）とを出力するよう構成されている。
【００３２】
本発明の実施の形態の液滴成長シミュレーション装置１０は、本発明の実施の形態の液滴成長シミュレーション方法により、流体中の液滴の粒径や分布の時間変化を推定することができる。図２に示すように、本発明の実施の形態の液滴成長シミュレーション方法は、まず、入力手段１１により、初期状態での液滴の数、粒径、密度や動粘度などの物性、分布（座標）、液滴に作用する重力加速度および重力の方向、流体の物性、温度、シミュレーションにおける計算のステップ数、積分時間（経過時間）、セルサイズなどの初期情報を入力する（ステップ２１）。なお、図２では、液滴として油滴についてシミュレーションする場合が記載されている。
【００３３】
初期情報を入力したならば、その初期情報から、シミュレーションに使用する各種の初期値を設定する（ステップ２２）。液滴計算手段１２により、ブラウン運動の計算を行うために乱数を発生させ（ステップ２３）、運動計算手段１２ａにより全ての液滴についてΔｔ経過後の位置を計算する（ステップ２４）。全ての液滴の位置を計算した後、液滴合体手段１２ｂにより液滴の合体処理を行い、液滴の粒径を計算する（ステップ２５）。所定の経過時間になるまで、時刻ステップΔｔ毎に運動計算手段１２ａと液滴合体手段１２ｂとを繰り返す（ステップ２６）。所定の経過時間後、出力手段１３により、液滴の座標（分布）、粒径、数などを出力する（ステップ２７）。
【００３４】
なお、本発明の実施の形態の液滴成長シミュレーションプログラムは、図１に示す液滴成長シミュレーション装置１０を構成するコンピュータを、図２に示すフローに従って機能させるプログラムから成っている。また、本発明の実施の形態の液滴成長シミュレーションプログラムは、例えば、ＣＤ（ＣＤ−ＲＯＭ、ＣＤ−Ｒ、ＣＤ−ＲＷなど）、ＤＶＤ（ＤＶＤ−ＲＯＭ、ＤＶＤ−ＲＡＭ、ＤＶＤ−Ｒ、ＤＶＤ−ＲＷ、ＤＶＤ＋Ｒ、ＤＶＤ＋ＲＷなど）等のコンピュータ読み取り可能な記録媒体に記録された形態で提供されてもよい。この場合、コンピュータは、その記録媒体から液滴成長シミュレーションプログラムを読み取ってコンピュータの内部記憶装置または外部記憶装置に転送し、格納して用いることができる。また、本発明の実施の形態の液滴成長シミュレーションプログラムを、例えば、磁気ディスク、光ディスク、光磁気ディスク等の記憶装置（記録媒体）に記録しておき、その記憶装置から通信回線を介してコンピュータに提供するようになっていてもよい。
【００３５】
本発明の実施の形態の液滴成長シミュレーションプログラムとしてのアプリケーションプログラムは、上述のようなコンピュータに実現させるプログラムコードを含んでいる。また、その機能の一部は、アプリケーションプログラムではなくＯＳによって実現されてもよい。なお、本発明の実施の形態の液滴成長シミュレーションプログラムを記録したコンピュータ読み取り可能な記録媒体としては、上述したフレキシブルディスク、ＣＤ、ＤＶＤ、磁気ディスク、光ディスク、光磁気ディスクのほか、コンピュータの内部記憶装置（ＲＡＭやＲＯＭなどのメモリ）、外部記憶装置等や、バーコードなどの符号が印刷された印刷物等の、コンピュータ読み取り可能な種々の媒体を利用することもできる。
【００３６】
本発明の実施の形態の液滴成長シミュレーション装置１０、液滴成長シミュレーション方法、液滴成長シミュレーションプログラムおよび、液滴成長シミュレーションプログラムを記録したコンピュータ読み取り可能な記録媒体は、ブラウン運動を考慮した（１）式のＳｔｏｋｅｓの抵抗法則を用いて液滴の運動を算出して液滴の分布を計算するだけでなく、（２）式のＨｏｏｋｅａｎの方程式を用いて液滴が合体成長する確率を求めることにより、液滴の物性や初期分布から、液滴の移動および合体による粒成長の挙動を考慮して、液滴の粒径や空間分布の時間変化を精度良く推定することができる。
【００３７】
また、（１）式の運動方程式での計算過程で、（７）式のファンデアワールス相互作用および（８）式の静電相互作用以外の相互作用を考慮しないため、高速でシミュレーションを行うことができる。また、ファンデアワールス相互作用のみ、または静電相互作用のみを考慮した場合に比べて、高精度でシミュレーションを行うことができる。比較的影響が大きいファンデアワールス相互作用および静電相互作用のみを考慮するため、高速で高精度なシミュレーションを行うことができる。また、（３）式を使用することにより、精度良く液滴の合体成長の確率を計算することができる。（４）式〜（６）式により、高精度かつ高速でシミュレーションを行うことができる。
【実施例１】
【００３８】
本発明の実施の形態の液滴成長シミュレーション装置１０および液滴成長シミュレーション方法を用いて、鉛直方向に細長い容器内における、水中での油滴（ＴＭＢ；１，３，５−トリメチルベンゼン）の粒径や分布の時間変化についてシミュレーションを行った。初期情報を、以下のように設定した。
系の温度                ４０［℃］
ＴＭＢの濃度            ３６０［ｍＭ］，４８０［ｍＭ］，６００［ｍＭ］
ＴＭＢ粒子の初期粒径    １２０［μｍ］
セルの初期サイズ        １２０×１２０×６００［μｍ］
重力加速度              １［Ｇ］，１［μＧ］
ＴＭＢ粒子の合体確率    １
流体の粘度              ６．５５×１０−４［Ｐａ・ｓ］
流体の密度              １．００×１０３［ｇ／Ｌ］
ＴＭＢの比重            ０．８６３７
重力の方向              鉛直下方：（ｘ，ｙ，ｚ）＝（０，０，−１）
【００３９】
また、重力加速度が１Ｇのとき、実際は容器の高さが有限であるため、油滴が容器上面に溜まると思われるが、今回の計算ではそのように扱わず、容器上面に達した油滴はそのまま上昇していく（計算上では消滅する）ものとした。
【００４０】
また、（２）〜（５）式の液滴合体手段１２ｂで使用するパラメータを、以下のように設定した。
曲げ係数κ：１．０×１０−２０［Ｊ］
鞍状係数  ：８．０×１０−２１［Ｊ］
δ＝−１．０×１０７［ｍ−１Ｋ−１］
Ｔ０＝３１．８［℃］
κ’α２＝３．６×１０−７［Ｊｍ−１Ｋ−１］
【００４１】
また、（７）式および（８）式の相互作用の計算で使用するパラメータを、以下のように設定した。
電解質濃度：          ｃ＝１３．５［ｍＭ］
Ｓｔｅｒｎ電位：      ΨＳ＝−３２．０［ｍＶ］
イオン価：            Ｚ＝１
Ｈａｍａｋｅｒ定数：  Ａ＝０．４１０×１０−２０［Ｊ］
溶液の比誘電率：      ８０．４
【００４２】
まず、重力加速度を１Ｇとし、油滴の濃度が３６０ｍＭおよび４８０ｍＭの２通りの場合についてシミュレーションを行った。このときの経過時刻と油滴の大きさの平均値との関係を、図３に示す。図３に示すように、油滴の濃度が４８０ｍＭの方が、３６０ｍＭの場合に比べて、経過時刻に従って油滴が大きく成長していく様子が確認できる。これは、油滴の濃度が高いほど、油滴が合体して成長する頻度が高くなるためであると考えられる。
【００４３】
次に、油滴の濃度を６００ｍＭとし、重力加速度が１μＧおよび１Ｇの２通りの場合についてシミュレーションを行った。このときのそれぞれの油滴の初期配置（０ｈ）、４時間経過後の油滴の分布状態（４ｈ）および８時間経過後の油滴の分布状態（８ｈ）を、図４（ａ）および（ｂ）に示す。図４に示すように、重力加速度が１Ｇの方が、１μＧの場合に比べて、油滴が早く上昇していることが確認できる。これは、重力加速度が大きいほど油滴に働く浮力が大きくなるためであると考えられる。
【００４４】
図３および図４に示すように、本発明の実施の形態の液滴成長シミュレーション装置１０および液滴成長シミュレーション方法によれば、実際の液滴の挙動に調和するシミュレーション結果が得られることから、液滴の物性や初期分布から液滴の粒径や空間分布の時間変化を精度良く推定することができるといえる。
【産業上の利用可能性】
【００４５】
本発明に係る液滴成長シミュレーション装置、液滴成長シミュレーション方法、液滴成長シミュレーションプログラムおよび、液滴成長シミュレーションプログラムを記録したコンピュータ読み取り可能な記録媒体は、宇宙ステーションなどの微小重力環境などで発生する特殊な液滴の粒成長や、地上での通常の重力環境下での油滴を含む各種液滴に関する粒成長シミュレーションに使用することができる。また、各種材料の溶液を用いた製造プロセスの設計に利用することもできる。
【符号の説明】
【００４６】
１０  液滴成長シミュレーション装置
１１  入力手段
１２  液滴計算手段
１２ａ  運動計算手段
１２ｂ  液滴合体手段
１３  出力手段
【図１】
2014049054.tif 000019
【図２】
2014049054.tif 000020
【図３】
2014049054.tif 000021
【図４】
2014049054.tif 000022