(57)【要約】【課題】  界面での占有率関数のにじみを抑制しつつ、より物理的に接触角を表現しうる表面張力の計算を行うための技術を提供すること。【解決手段】  分布関数ｆ、ｇについて、ランダウ・ギンツブルグ・ポテンシャルによる気液界面のエネルギー汎関数と、固液界面の面積を近似的にエネルギーとするエネルギー汎関数と、固気界面の面積を近似的にエネルギーとするエネルギー汎関数の３者の総和からなるエネルギー汎関数を定義し、定義されたエネルギー汎関数を分布関数ｆ、ｇに対して最小化することにより得られる微分方程式を差分化した式に従った演算処理を行うことで表面張力を近似的に計算する（Ｓ１０３）。【選択図】  図１
2006284456.tif 000002
【特許請求の範囲】
【請求項１】
固体、液体、気体の三相が共存する空間における流れを、仮想空間内に設けた格子点群で構成される計算領域内でもって計算する情報処理装置であって、
前記計算領域内における液体及び固体の分布関数を定義した場合に、隣接する格子点間を結ぶ線分で囲まれた小領域毎に、小領域の位置に応じた前記液体及び固体の分布関数のデータを保持する保持手段と、
前記保持手段が保持するデータを用いて、前記液体及び固体の分布関数について、気液界面のエネルギー汎関数と、固液界面のエネルギー汎関数と、固気界面のエネルギー汎関数の３者の総和からなるエネルギー汎関数を定義し、当該定義されたエネルギー汎関数を前記液体及び固体の分布関数に対して最小化することにより得られる微分方程式を差分化した式に従った演算処理を行うことで、前記小領域毎に、小領域の中心から外部に向かう力を、表面張力として近似的に計算する計算手段と
を備えることを特徴とする情報処理装置。
【請求項２】
前記気液界面のエネルギー汎関数は、ランダウ・ギンツブルグ・ポテンシャルによるエネルギー汎関数であることを特徴とする請求項１に記載の情報処理装置。
【請求項３】
前記固液界面のエネルギー汎関数は、固液界面の面積を近似的にエネルギーとするエネルギー汎関数であることを特徴とする請求項１に記載の情報処理装置。
【請求項４】
前記気液界面のエネルギー汎関数は、気液界面の面積を近似的にエネルギーとするエネルギー汎関数であることを特徴とする請求項１に記載の情報処理装置。
【請求項５】
固体、液体、気体の三相が共存する空間における流れを、仮想空間内に設けた格子点群で構成される計算領域内でもって計算する情報処理装置が行う情報処理方法であって、
前記計算領域内における液体及び固体の分布関数を定義した場合に、隣接する格子点間を結ぶ線分で囲まれた小領域毎に、小領域の位置に応じた前記液体及び固体の分布関数のデータを保持する保持工程と、
前記データを用いて、前記液体及び固体の分布関数について、気液界面のエネルギー汎関数と、固液界面のエネルギー汎関数と、固気界面のエネルギー汎関数の３者の総和からなるエネルギー汎関数を定義し、当該定義されたエネルギー汎関数を前記液体及び固体の分布関数に対して最小化することにより得られる微分方程式を差分化した式に従った演算処理を行うことで、前記小領域毎に、小領域の中心から外部に向かう力を、表面張力として近似的に計算する計算工程と
を備えることを特徴とする情報処理方法。
【請求項６】
コンピュータに請求項５に記載の情報処理方法を実行させることを特徴とするプログラム。
【請求項７】
請求項６に記載のプログラムを格納したことを特徴とする、コンピュータ読み取り可能な記憶媒体。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、固体、液体、気体の三相が共存する流れを数値的に模擬する方法為の技術に関するものであり、特に、各相間の界面における表面張力を近似的に計算する為の技術に関するものである。
【背景技術】
【０００２】
異なる相が混在する流れを数値的に解く手法のひとつとして、占有率関数を用いる手法が知られている。占有率関数とは、各格子点における気体、液体、固体のそれぞれの占有率をあらわしたものである。この占有率関数を用いた手法は、ＶＯＦ法（非特許文献１を参照）の流れをくむものであり、液滴の分裂・合体といった問題や、界面が大変形するような問題に対しても適用可能であることから、現在に至っても種々の改良を施されながら用いられている。
【０００３】
占有率関数を用いた手法のもとで表面張力を取り扱う手法として、よく知られたものにＣＳＦモデル（非特許文献２を参照）がある。ＣＳＦモデルの考え方は次のようなものである。簡単のため２つの物質、たとえば水と空気からなる系を考え、物質１（たとえば水）の占有率関数をｆとする。この系全体のエネルギー汎関数Ｅ１２を、
【０００４】
【数１】
2006284456.tif 000003
【０００５】
のように、物質１と物質２の界面の面積に比例するとして定義する。なおｄｖｏｌは微小体積要素であって、以後の数式においても同様である。ここで、比例係数σ１２が物質１と物質２の界面に働く表面エネルギーの係数と考えれば、結局、エネルギー汎関数Ｅ１２は物質１と物質２の界面の界面エネルギーを近似的に表現しているといえる。よって、ＣＳＦモデルにおける表面張力は、界面の表面積を最小にしようとするいわゆる最小原理から、エネルギー汎関数Ｅ１２を最小化するように働く力として定義される。
【０００６】
一方、占有率関数を用いた手法におけるよく知られた問題点として、占有率関数の等高線の間隔が計算ステップの増加とともに広がってしまい界面の位置が不明瞭になる、といういわゆる占有率関数のにじみの問題がある。この問題の解決策として近年注目されているのが、ランダウ・ギンツブルグ・ポテンシャルを用いる手法（非特許文献３）である。これは先のＣＳＦモデルと同様、２つの物質からなる系において、物質１の占有率関数をｆとし、エネルギー汎関数Ｅ１２を
【０００７】
【数２】
2006284456.tif 000004
【０００８】
と定義する手法である。この手法における表面張力は、先のＣＳＦモデルの場合と同様にエネルギー汎関数Ｅ１２を最小化するように働く力として定義される。同時に、ポテンシャルφをエネルギー汎関数Ｅ１２のｆに対する変分として
【０００９】
【数３】
2006284456.tif 000005
【００１０】
のように定義し、このポテンシャルφを用いて、移流方程式を
【００１１】
【数４】
2006284456.tif 000006
【００１２】
のように修正する。ここでα、βは占有率関数のにじみ具合および表面張力の大きさを制御するためのパラメータ、κは適当な時定数であって、これらを適当に制御することにより、占有率関数のにじみを抑制することができる。また、ξとしては２つの極小をもつような関数、たとえば
【００１３】
【数５】
2006284456.tif 000007
【００１４】
が用いられる。この計算手法を用いて、２流体の系の場合に限っては信頼に足る計算結果が得られている。
【００１５】
しかしながら、上記説明したようなランダウ・ギンツブルグ・ポテンシャルを用いる手法をもってしても、依然として次のような問題点が残されている。
【００１６】
第一に、上記非特許文献２にも示されているように、三相界面を含む系においてダイナミクスが絡んできた場合等には、計算方法に任意性が入り込む問題がある。上記説明したような方法は、接触角を幾何学的に捕らえるという観点から構成されているが、非定常状態においては接触角に物理的な意味はないからである。ここで接触角とは、図６に示したような三相界面（８０２は液滴、８０３は固体壁、８０４は空気）における角度８０１として定義される角度のことであって、三相界面の特徴をとらえるためのよい指標となる量である。
【００１７】
なお、ダイナミクスが絡んできた場合における接触角の取り扱いについては、三相界面近傍の微小領域においてレベルセット関数を用いて接触角を監視し、あらかじめパラメータとして与えてある臨界接触角との関係から境界条件を場合分けするという方法（特許文献１）も提案されている。
【非特許文献１】C. W. Hirt and B. D. Nichols, Volume of fluid (VOF) method for the dynamics of free boundaries, J. Comp. Phys., 39, 201(1981)
【非特許文献２】J. U. Brackbill, D. B. Kothe, and C. Zemach, A continuum method for modeling surface tension, J. Comp. Phys, 100, 335(1992)
【非特許文献３】D. Jacqmin, Calculation of Two-Phase Navier-Stokes Flows Using Phase-Field Modeling, J. Comp. Phys., 155, 96(1999)
【特許文献１】特開２００３−２８５４２８号公報
【発明の開示】
【発明が解決しようとする課題】
【００１８】
しかしながら特許文献１に開示されている手法では、レベルセット関数を用いるため、特に３次元化した場合に計算コストが高くついてしまうほか、三相界面近傍の微小領域の取り扱いや、人工的に導入される表面張力等、理論的に不明確な点がいくつか見られ、その信頼性には疑問が残る。
【００１９】
また、重力等の力を受けている系の問題がある。図７に示したように、空気９０６中におけるある斜度９０１をなす傾いた固体壁９０５上の液滴９０２は、重力の影響により接触角９０３と接触角９０４が異なる値を持つと考えるのが自然である。しかしながら、上記手法ではこの違いを取り込むことができない。
【００２０】
本発明は以上の問題に鑑みてなされたものであり、界面での占有率関数のにじみを抑制しつつ、より物理的に接触角を表現しうる表面張力の計算を行うための技術を提供することを目的とする。
【課題を解決するための手段】
【００２１】
本発明の目的を達成するために、例えば本発明の情報処理装置は以下の構成を備える。
【００２２】
即ち、固体、液体、気体の三相が共存する空間における流れを、仮想空間内に設けた格子点群で構成される計算領域内でもって計算する情報処理装置であって、
前記計算領域内における液体及び固体の分布関数を定義した場合に、隣接する格子点間を結ぶ線分で囲まれた小領域毎に、小領域の位置に応じた前記液体及び固体の分布関数のデータを保持する保持手段と、
前記保持手段が保持するデータを用いて、前記液体及び固体の分布関数について、気液界面のエネルギー汎関数と、固液界面のエネルギー汎関数と、固気界面のエネルギー汎関数の３者の総和からなるエネルギー汎関数を定義し、当該定義されたエネルギー汎関数を前記液体及び固体の分布関数に対して最小化することにより得られる微分方程式を差分化した式に従った演算処理を行うことで、前記小領域毎に、小領域の中心から外部に向かう力を、表面張力として近似的に計算する計算手段と
を備えることを特徴とする。
【００２３】
本発明の目的を達成するために、例えば本発明の情報処理方法は以下の構成を備える。
【００２４】
即ち、固体、液体、気体の三相が共存する空間における流れを、仮想空間内に設けた格子点群で構成される計算領域内でもって計算する情報処理装置が行う情報処理方法であって、
前記計算領域内における液体及び固体の分布関数を定義した場合に、隣接する格子点間を結ぶ線分で囲まれた小領域毎に、小領域の位置に応じた前記液体及び固体の分布関数のデータを保持する保持工程と、
前記データを用いて、前記液体及び固体の分布関数について、気液界面のエネルギー汎関数と、固液界面のエネルギー汎関数と、固気界面のエネルギー汎関数の３者の総和からなるエネルギー汎関数を定義し、当該定義されたエネルギー汎関数を前記液体及び固体の分布関数に対して最小化することにより得られる微分方程式を差分化した式に従った演算処理を行うことで、前記小領域毎に、小領域の中心から外部に向かう力を、表面張力として近似的に計算する計算工程と
を備えることを特徴とする。
【発明の効果】
【００２５】
本発明の構成により、界面での占有率関数のにじみを抑制しつつ、より物理的に接触角を表現しうる表面張力の計算を行うことができる。
【発明を実施するための最良の形態】
【００２６】
以下添付図面を参照して、本発明を好適な実施形態に従って詳細に説明する。
【００２７】
先ず、本発明に係る表面張力計算方法について説明する。本発明に係る表面張力計算方法は、気体、液体、固体の三相が共存する二次元以上の系に対して定義される、以下に示すエネルギー汎関数Ｅを基礎としている。
【００２８】
【数６】
2006284456.tif 000008
【００２９】
ここで、ｆ、ｇは［０，１］の値を持つ関数で、系全体で定義されている。ｆｇは液体の占有率を表し、ｆｇが１のとき、その場所には液体が存在しているとする。同様に、（１−ｆ）ｇは気体の占有率関数を表し、（１−ｆ）ｇが１のとき、その場所には気体が存在しているとする。またｇは固体の占有率関数を表し、ｇが０のとき、その場所には固体が存在しているとする。なお、σ２は固体と液体の間の表面張力係数、σ３は気体と固体の間の表面張力係数である。
【００３０】
上記エネルギー汎関数の第二項の積分は固液の界面の表面積を、同様に第三項の積分は固気の界面の表面積を、それぞれ近似的に表しており、表面積が小さい方がエネルギーが低いことを意味している。また第一項の積分は、ランダウ・ギンツブルグ・ポテンシャルの項であって、
【００３１】
【数７】
2006284456.tif 000009
【００３２】
のように定義されている。上式を見るとわかるように、エネルギー汎関数の第一項は、気液の界面の表面積を近似的に表す項である。また、上式中のα、βは占有率関数のにじみと表面張力の強さに関係するパラメータである。
【００３３】
続いてこのエネルギー汎関数に対し、ネーターの定理を適用する。すると、
【００３４】
【数８】
2006284456.tif 000010
【００３５】
を得る。なお、ネーターの定理に関してはたとえば「量子場を学ぶための場の解析力学入門」（高橋康著，講談社サイエンティフィック１９８２年発行）に詳しい説明がある。δＥ／δｘｉを系に対してｉ方向に与えた摂動と考えれば、上式右辺は摂動に対する反作用による力を表していると考えることができる。本発明におけるエネルギー汎関数は界面のエネルギーを近似的に表しており、この界面のエネルギーは近似的に界面の面積に比例している。したがってこの反作用による力は、表面積を最小にする最小原理から導かれた力、すなわち表面張力であるということができる。
【００３６】
［第１の実施形態］
図２は、本実施形態に係る情報処理装置として機能するコンピュータの基本構成を示すブロック図である。なお、このコンピュータには、例えば一般のＰＣ（パーソナルコンピュータ）やＷＳ（ワークステーション）等が適用可能であるし、その他にも、スーパーコンピュータなどの計算専用装置を適用することのもちろん可能である。
【００３７】
２０１はＣＰＵで、メモリ２０２にロードされているプログラムやデータを用いて本コンピュータ全体の制御を行うと共に、コンピュータが行う後述の計算処理を実行する。
【００３８】
２０２はメモリで、ＨＤＤ２０３からロードされたプログラムやデータを一時的に記録するためのエリアを備えると共に、ＣＰＵ２０１が各処理を実行する際に使用するワークエリアを備える。
【００３９】
２０３はＨＤＤ（ハードディスクドライブ）で、ここにＯＳ（オペレーティングシステム）や、ＣＰＵ２０１にコンピュータが行う後述の計算処理を実行させるためのプログラムやデータが保存されており、これらの一部もしくは全部はＣＰＵ２０１の制御に従ってメモリ２０２にロードされる。
【００４０】
２０４は表示装置で、ＣＲＴや液晶画面などにより構成されており、ＣＰＵ２０１による実行結果を文字や画像でもって表示することができる。
【００４１】
２０５は入力装置で、キーボードやマウス等により構成されており、本コンピュータの操作者が操作することで、ＣＰＵ２０１に対して各種の指示を入力することができる。
【００４２】
２０６は上述の各部を繋ぐバスである。
【００４３】
次に、このような構成を有するコンピュータによって実行される計算処理について説明する。ここでの計算処理とは、固体、液体、気体の三相が共存する空間における流れを、仮想空間内に設けた格子点群で構成される計算領域内でもって計算する処理のことである
図１は、このような計算処理のフローチャートである。なお、同図のフローチャートに従った処理をＣＰＵ２０１に実行させるためのプログラムやデータはＨＤＤ２０３に保存されており、これをＣＰＵ２０１の制御に基づいてメモリ２０２にロードし、ＣＰＵ２０１の処理対象とすることで、本コンピュータは以下説明する各処理を実行することになる。
【００４４】
先ず、同図のフローチャートに示す如く、この計算処理は、初期化処理（ステップＳ１０２）、与えられた分布関数ｆ、ｇに対する表面張力計算処理（ステップＳ１０３）、表面張力を含む流体計算処理（ステップＳ１０４）、時間更新処理（ステップＳ１０５）、終了判定処理（ステップＳ１０６）により構成されている。以下では、このステップＳ１０３における処理、即ち、表面張力計算処理について詳細に説明する。
【００４５】
以下の説明では、計算領域として２次元平面を考える。従って以下の説明では、２次元の系における、上記三相が共存する空間の流れを計算することになる。よってこのために、メモリ２０２には、２次元平面上に配置された個々の格子点（スタガード格子）の位置データ等、個々の格子点に係るデータが保持されている。格子点に係るデータについては以下の各処理の説明に伴って説明する。
【００４６】
図３は、２次元平面上における格子群の一部を示す図である。同図において３０３〜３０６はそれぞれ格子で、４つの格子で囲まれた領域３０１を以下「胞」と呼称する。また、この胞３０１は４つの辺により構成されており、３０２はその辺の１つを示す。
【００４７】
胞３０１には圧力や占有率関数ｆ，ｇ等のスカラー量が定義される。胞の位置はｘ方向にｍ番目、ｙ方向にｎ番目というラベル付けを行うので、胞の位置は（ｍ，ｎ）のように表される。上述の通り、占有率関数ｆ、ｇは計算領域全体で定義されるものであるので、個々の胞に対する占有率関数ｆ、ｇのスカラー量のデータは、その胞の位置に応じたものとなる。
【００４８】
また、辺３０２には、速度や運動量などのベクトル量がそれぞれ定義される。この辺の位置は（ｍ＋１／２，ｎ）のように表される。
【００４９】
また格子点３０３〜３０６には、必要に応じて種々の物理量が補間されて配置される。格子点３０３を例に取ると、その位置は（ｍ＋１／２，ｎ＋１／２）と表される。格子間隔はｘ方向、ｙ方向共に一様であるとする。
【００５０】
このように、胞の位置を用いて格子点、辺の位置が決まると共に、胞、辺には上述の物理量が定義される。そしてそれぞれの格子点にはそれぞれの物理量が補間されて定義される。なお、胞、辺、そして格子点のそれぞれについて定義された上述の物理量はデータとしてＨＤＤ２０３に保存されており、必要に応じてメモリ２０２にロードされるものである。
【００５１】
上記式（８）によって定義される表面張力は微小体積要素あたりに働く力、すなわち体積力である。この様子を図４に示す。図４は体積力を説明する図である。同図において４０１は微小体積要素で、４０２，４０３はそれぞれｘ方向、ｙ方向への表面張力Ｆｘ、Ｆｙを示す。本実施形態では胞を微小体積要素として用いるので、４０１は胞である。また、表面張力Ｆｘ、Ｆｙは上記式（８）を用いて以下の式でもって表すことができる。
【００５２】
【数９】
2006284456.tif 000011
【００５３】
この表面張力Ｆｘ、Ｆｙは、スタガード格子の辺上に定義されるものであり、より具体的には、辺で囲む胞の中心から辺に直交する方向（即ち胞の外部に向けて）に表面張力を定義する。
【００５４】
本実施形態ではこの表面張力Ｆｘ、Ｆｙを求めるのであるが、Ｆｘ、Ｆｙが対象性を有することから以下ではＦｘを求める方法についてのみ説明するが、Ｆｙを求める方法についても同様である。
【００５５】
先ず、ｘ方向の表面張力Ｆｘを具体的に書き下すと、以下の式（１０）のようになる。
【００５６】
【数１０】
2006284456.tif 000012
【００５７】
ここで、上式の各項の微分演算子が作用する部分（カッコ内）を順にそれぞれτ１１、τ１２、η１、ξ１１、ξ１２、χ１とする。すなわち式（１０）を
【００５８】
【数１１】
2006284456.tif 000013
【００５９】
のように書き換える。ここでτ１１、η１、ξ１１、χ１は胞に、τ１２、ξ１２は格子に定義される値であるとする。これらの値がそれぞれ微分演算を受け、ｘ方向の表面張力として図４におけるＦｘの値となる。このことを考慮し、図３に従い、式（１１）の微分演算子∂１、∂２の差分式を以下のように定義する。
【００６０】
【数１２】
2006284456.tif 000014
【００６１】
ここでκは任意の物理量であり、以後も同様である。同様にしてξ１１、χ１の中にある、胞に配置される勾配の絶対値の差分式を
【００６２】
【数１３】
2006284456.tif 000015
【００６３】
と定義する。また、ξ１２の中にある、格子に配置される勾配の絶対値の差分式は、
【００６４】
【数１４】
2006284456.tif 000016
【００６５】
と定義される。さらに、τ１１、η１、ξ１１に含まれる微分演算子∂１の差分式を、
【００６６】
【数１５】
2006284456.tif 000017
【００６７】
と定義する。同様にτ１２，ξ１２に含まれるふたつの微分演算子∂１、∂２の差分式をそれぞれ、
【００６８】
【数１６】
2006284456.tif 000018
【００６９】
と定義する。また、τ１２、ξ１２における占有率関数ｇは、単純な平均操作
【００７０】
【数１７】
2006284456.tif 000019
【００７１】
によって定義される。
【００７２】
このようにして得られた差分式に従って計算処理を行うことで、表面張力Ｆｘの値を計算することができる。また上述の通りＦｙについても、Ｆｘとの対称性を考慮して上記と同様の差分化を行うことにより計算ができる。
【００７３】
上述した手法を用い、実際に計算を行った例を図５に示した。図５は重力を受けている系における斜面上の液滴を計算した結果を示す図であって、図７と対応するものである。５０１は液滴、５０２は固体壁、５０３は空気であり、重力の働く方向を矢印５０４で示した。このように、上述の表面張力計算方法を用いることにより、界面のにじみを低減しつつ、二つの異なる接触角が計算できた。
【００７４】
以上のようにして、式（６）をもとに導出された式（８）によって定義される表面張力を二次元スタガード格子上で計算でき、ステップＳ１０３における処理がなされる。
【００７５】
なお、以上の説明は計算領域を２次元平面としたので、添え字の組を（ｍ，ｎ）＝｛１，２｝としたが、計算領域を３次元空間へと拡張するためには、添え字の組を（ｍ，ｎ）＝｛１，２，３｝とするだけでよい。また、差分化についても、以上の説明から容易に類推することが可能であるため、三次元の計算への適用も容易に行うことができる。
【００７６】
［第２の実施形態］
なお、式（８）の第一項を−ｆ▽φに置き換えたものを表面張力の式として用いるようにしても良い。ここでφは式（３）で定義されるポテンシャルである。
【００７７】
［第３の実施形態］
なお、式（８）の第二項を０としたものを表面張力の式として用いるようにしても良い。このように、項の省略により、計算時間の短縮が期待される。
【００７８】
また、以上の説明の趣旨に添ったものであれば、第２の実施形態、第３の実施形態に示したような表面張力の式の変形を行っても同様の結果が得られることは自明である。また、表面張力の表式の差分化においても、上記実施形態に限定されるものではなく、差分法のみならず有限体積法や有限要素法に対しても以上の実施形態を適用することが可能である。
【００７９】
以上のことから、界面での占有率関数のにじみを抑制しつつ、より物理的に接触角を表現しうる表面張力の計算が可能になる。特に上記各実施形態においては、エネルギーの観点から表面張力を評価しているため、非平衡状態においてや、重力の働きを考慮した場合においても正しく表面張力の効果を計算することが可能になる。
【００８０】
［その他の実施形態］
また、本発明の目的は、前述した実施形態の機能を実現するソフトウェアのプログラムコードを記録した記録媒体（または記憶媒体）を、システムあるいは装置に供給し、そのシステムあるいは装置のコンピュータ（またはＣＰＵやＭＰＵ）が記録媒体に格納されたプログラムコードを読み出し実行することによっても、達成されることは言うまでもない。この場合、記録媒体から読み出されたプログラムコード自体が前述した実施形態の機能を実現することになり、そのプログラムコードを記録した記録媒体は本発明を構成することになる。
【００８１】
また、コンピュータが読み出したプログラムコードを実行することにより、前述した実施形態の機能が実現されるだけでなく、そのプログラムコードの指示に基づき、コンピュータ上で稼働しているオペレーティングシステム（ＯＳ）などが実際の処理の一部または全部を行い、その処理によって前述した実施形態の機能が実現される場合も含まれることは言うまでもない。
【００８２】
さらに、記録媒体から読み出されたプログラムコードが、コンピュータに挿入された機能拡張カードやコンピュータに接続された機能拡張ユニットに備わるメモリに書込まれた後、そのプログラムコードの指示に基づき、その機能拡張カードや機能拡張ユニットに備わるＣＰＵなどが実際の処理の一部または全部を行い、その処理によって前述した実施形態の機能が実現される場合も含まれることは言うまでもない。
【００８３】
本発明を上記記録媒体に適用する場合、その記録媒体には、先に説明したフローチャートに対応するプログラムコードが格納されることになる。
【図面の簡単な説明】
【００８４】
【図１】固体、液体、気体の三相が共存する空間における流れを、仮想空間内に設けた格子点群で構成される計算領域内でもって計算する処理のフローチャートである。
【図２】本発明の第１の実施形態に係る情報処理装置として機能するコンピュータの基本構成を示すブロック図である。
【図３】２次元平面上における格子群の一部を示す図である。
【図４】体積力を説明する図である。
【図５】重力を受けている系における斜面上の液滴を計算した結果を示す図である。
【図６】接触角を示す図である。
【図７】空気９０６中におけるある斜度９０１をなす傾いた固体壁９０５上の液滴９０２を示す図である。
【図１】
2006284456.tif 000020
【図２】
2006284456.tif 000021
【図３】
2006284456.tif 000022
【図４】
2006284456.tif 000023
【図５】
2006284456.tif 000024
【図６】
2006284456.tif 000025
【図７】
2006284456.tif 000026