(57)【要約】【課題】キャビティにおける流体の流動解析において、計算負荷をより低減させる。【解決手段】流体が流動するキャビティを複数の微小要素に分割する分割部と、前記流体に関する物理量の保存則を表す方程式を、前記キャビティの厚さの分布に基づいて、調整する調整部と、前記調整部により調整された前記方程式を用いて、前記微小要素の各々について、流速を算出する解析部と、を備える、情報処理装置を提供する。【選択図】図４
2017162269.tif 000002
【特許請求の範囲】
【請求項１】
流体が流動するキャビティを複数の微小要素に分割する分割部と、
前記流体に関する物理量の保存則を表す方程式を、前記キャビティの厚さの分布に基づいて、調整する調整部と、
前記調整部により調整された前記方程式を用いて、前記微小要素の各々について、流速を算出する解析部と、
を備える、情報処理装置。
【請求項２】
前記調整部は、前記方程式の空間微分項における被演算子に前記厚さの空間微分及び前記厚さの逆数を乗じて得られる項の前記方程式への付加に相当する調整を行う、請求項１に記載の情報処理装置。
【請求項３】
前記調整部は、前記方程式の時間微分項における被演算子に前記厚さの時間微分及び前記厚さの逆数を乗じて得られる項の前記方程式への付加に相当する調整を行う、請求項１に記載の情報処理装置。
【請求項４】
2017162269.tif 000038
【数１】
2017162269.tif 000039
【請求項５】
2017162269.tif 000040
【数２】
2017162269.tif 000041
【請求項６】
前記方程式は、前記流体の質量保存則を表す連続の式、及び前記流体の運動量保存則を表すナビエストークスの運動方程式を含む、請求項１に記載の情報処理装置。
【請求項７】
前記連続の式、及び前記ナビエストークスの運動方程式は、前記流体の湧き出しに関する項を含む、請求項６に記載の情報処理装置。
【請求項８】
前記ナビエストークスの運動方程式は、前記流体へ掛かる重力に関する項を含む、請求項６に記載の情報処理装置。
【請求項９】
前記ナビエストークスの運動方程式は、前記流体の流動抵抗に関する項を含む、請求項６に記載の情報処理装置。
【請求項１０】
前記ナビエストークスの運動方程式は、前記流体の表面張力に関する項を含む、請求項６に記載の情報処理装置。
【請求項１１】
前記調整部は、前記方程式における前記厚さと相関を有する項に前記厚さの逆数を乗じることによって、当該項を正規化する調整を行う、請求項１に記載の情報処理装置。
【請求項１２】
前記解析部は、前記微小要素の各々について、圧力を算出する、請求項１に記載の情報処理装置。
【請求項１３】
請求項１に記載の情報処理装置を備える、電子機器。
【請求項１４】
流体が流動するキャビティを複数の微小要素に分割することと、
前記流体に関する物理量の保存則を表す方程式を、前記キャビティの厚さの分布に基づいて、情報処理装置によって、調整することと、
調整された前記方程式を用いて、前記微小要素の各々について、流速を算出することと、
を含む情報処理方法。
【請求項１５】
コンピュータを、
流体が流動するキャビティを複数の微小要素に分割する分割部と、
前記流体に関する物理量の保存則を表す方程式を、前記キャビティの厚さの分布に基づいて、調整する調整部と、
前記調整部により調整された前記方程式を用いて、前記微小要素の各々について、流速を算出する解析部と、
として機能させるためのプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本開示は、情報処理装置、電子機器、情報処理方法、及びプログラムに関する。
【背景技術】
【０００２】
近年、射出成形における溶融した樹脂の流動に代表されるキャビティにおける流体の流動についての流動解析が広く行われている。また、そのような流動解析に関する技術が提案されている。
【０００３】
例えば、特許文献１では、２平面で囲まれた空間の熱流体解析を一段と容易に行うために、平行２平面間の流体に作用する流動抵抗を流体に作用する外力として定義すると共に上記２平面を流体に抵抗を与えないすべり壁として定義し、定義された条件に基づいて上記２平面間の流体及び熱の流れを解析する技術が提案されている。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特開２０００−３２９７８０号公報
【発明の概要】
【発明が解決しようとする課題】
【０００５】
ところで、流体の流動解析に関する分野において、計算負荷をより低減させることが望ましいと考えられる。例えば、キャビティの各位置における厚さが均一でない場合等に、解析精度を確保するために、解析モデルにおいて当該キャビティを比較的細かい微小要素に分割する必要が生じる場合がある。それにより、微小要素の数が増加することによって、計算負担が増大し得る。
【０００６】
そこで、本開示では、キャビティにおける流体の流動解析において、計算負荷をより低減させることが可能な、新規かつ改良された情報処理装置、電子機器、情報処理方法、及びプログラムを提案する。
【課題を解決するための手段】
【０００７】
本開示によれば、流体が流動するキャビティを複数の微小要素に分割する分割部と、前記流体に関する物理量の保存則を表す方程式を、前記キャビティの厚さの分布に基づいて、調整する調整部と、前記調整部により調整された前記方程式を用いて、前記微小要素の各々について、流速を算出する解析部と、を備える、情報処理装置が提供される。
【０００８】
また、本開示によれば、当該情報処理装置を備える、電子機器が提供される。
【０００９】
また、本開示によれば、流体が流動するキャビティを複数の微小要素に分割することと、前記流体に関する物理量の保存則を表す方程式を、前記キャビティの厚さの分布に基づいて、情報処理装置によって、調整することと、調整された前記方程式を用いて、前記微小要素の各々について、流速を算出することと、を含む情報処理方法が提供される。
【００１０】
また、本開示によれば、コンピュータを、流体が流動するキャビティを複数の微小要素に分割する分割部と、前記流体に関する物理量の保存則を表す方程式を、前記キャビティの厚さの分布に基づいて、調整する調整部と、前記調整部により調整された前記方程式を用いて、前記微小要素の各々について、流速を算出する解析部と、として機能させるためのプログラムが提供される。
【発明の効果】
【００１１】
以上説明したように本開示によれば、キャビティにおける流体の流動解析において、計算負荷をより低減させることが可能である。
【００１２】
なお、上記の効果は必ずしも限定的なものではなく、上記の効果とともに、又は上記の効果に代えて、本明細書に示されたいずれかの効果又は本明細書から把握され得る他の効果が奏されてもよい。
【図面の簡単な説明】
【００１３】
【図１】キャビティにおける流体の流動解析における要素分割の一例について説明するための説明図である。
【図２】解析対象のキャビティの一例を示す概念図である。
【図３】キャビティにおける流体の流動解析における要素分割の一例について説明するための説明図である。
【図４】本開示の実施形態に係る情報処理装置の機能構成の一例を示す説明図である。
【図５】キャビティの各位置における実際の厚さ及び質量流束密度について説明するための説明図である。
【図６】同実施形態に係る情報処理装置による流動解析における解析モデルの一例を示す概念図である。
【図７】キャビティにおける流体に掛かる表面張力について説明するための説明図である。
【図８】キャビティの所定の位置における実際の厚さと解析モデルにおける厚さとの差異について説明するための説明図である。
【図９】キャビティの所定の位置における厚さの経時変化を示す説明図である。
【図１０】同実施形態に係る情報処理装置が行う流動解析処理の流れの一例を示すフローチャートである。
【図１１】同実施形態に係る情報処理装置が行う方程式調整処理の流れの一例を示すフローチャートである。
【図１２】本開示に係る情報処理装置のハードウェア構成の一例を示す説明図である。
【発明を実施するための形態】
【００１４】
以下に添付図面を参照しながら、本開示の好適な実施の形態について詳細に説明する。なお、本明細書及び図面において、実質的に同一の機能構成を有する構成要素については、同一の符号を付することにより重複説明を省略する。
【００１５】
なお、説明は以下の順序で行うものとする。
１．導入
２．情報処理装置
２−１．機能構成
２−２．方程式の調整
３．動作
４．ハードウェア構成
５．まとめ
【００１６】
＜１．導入＞
キャビティにおける流体の流動解析において、計算負荷を低減させるために、微小要素の分割数を低減することが考えられる。例えば、キャビティの厚さ方向への流動の影響が比較的小さくなる程度の厚さを有する薄肉部について、二次元的な解析モデルを適用することにより、微小要素への分割である要素分割における分割数を低減することが提案されている。
【００１７】
具体的には、特許文献１では、２平面によって挟まれた上記薄肉部に相当する領域において、流動方向に直交する厚さ方向について速度分布を均一なものとして取り扱うことによって、厚さ方向についての要素分割を省略することを提案している。それにより、図１に示すように、厚さ方向について複数の微小要素への分割が行われた解析モデルＭ１０と比較して、厚さ方向についての要素分割が省略された二次元的な解析モデルＭ２０では、要素分割における分割数を低減することができる。なお、図１では、各解析モデルにおける流体に相当する部分Ｂ１、キャビティの分割により得られた各微小要素Ｂ３、及び内部空間としてのキャビティを形成する金型Ｂ５が示されている。
【００１８】
しかしながら、キャビティの薄肉部の各位置における厚さが均一でない場合には、上述した二次元的な解析モデルを適用することが困難な場合がある。例えば、図２に概念的に示すように、キャビティ１０の薄肉部１２０において、互いに厚さが異なる第１の領域１２２と第２の領域１２４が隣接して設けられ、第１の領域１２２と第２の領域１２４との境界に段差部Ｃ１０が形成される場合がある。このような場合には、段差部Ｃ１０の近傍において、解析精度を確保するために、比較的細かい微小要素への要素分割を行う必要が生じ得る。ゆえに、二次元的な解析モデルを適用することが困難な場合があるので、図３に示すように、厚さ方向について複数の微小要素への分割が行われた解析モデルＭ３０を用いて流動解析を行う必要があった。
【００１９】
それにより、流動解析における要素分割における分割数が増大することによって、計算負担が増大する場合があった。そこで、本明細書では、キャビティにおける流体の流動解析において、計算負荷をより低減させることが可能な仕組みを提案する。
【００２０】
＜２．情報処理装置＞
以下では、本開示の実施形態に係る情報処理装置４０について説明する。まず、情報処理装置４０の概要について説明した後に、図４を参照して、情報処理装置４０が有する各機能について説明する。
【００２１】
情報処理装置４０は、キャビティにおける流体の流動解析を行い、解析結果を出力する。情報処理装置４０の流動解析における解析対象の流体は、例えば、溶融樹脂であってもよい。なお、解析対象の流体は、特に限定されず、非圧縮性流体であってもよく、圧縮性流体であってもよい。また、解析対象の流体は、ニュートン流体であってもよく、非ニュートン流体であってもよい。情報処理装置４０は、例えば、オペレータから入力された入力情報を用いて、流動解析を行い、解析結果を表示させることによって、オペレータへ当該解析結果を提示する。
【００２２】
情報処理装置４０は、例えば、オペレータからの入力を受け付ける入力装置から入力された情報を受信することによって、入力情報を取得し得る。その場合には、当該入力装置と情報処理装置４０とは有線又は無線のネットワークを介して通信可能に構成される。また、情報処理装置４０は、当該入力装置と一体として構成されてもよい。なお、情報処理装置４０は、オペレータにより入力された入力情報以外の情報を用いて流動解析を行ってもよい。例えば、情報処理装置４０は、予め情報処理装置４０の記憶素子に記憶された情報を用いて、流動解析を行ってもよい。
【００２３】
また、情報処理装置４０は、例えば、画像を表示する表示装置へ解析結果を送信することによって、当該表示装置に解析結果を表示させ得る。その場合には、当該表示装置と情報処理装置４０とは有線又は無線のネットワークを介して通信可能に構成される。また、情報処理装置４０は、当該表示装置と一体として構成されてもよい。なお、情報処理装置４０は、解析結果の表示以外の方法によって、オペレータへ当該解析結果を提示してもよい。例えば、情報処理装置４０は、解析結果を印刷装置に印刷させることによって、オペレータへ当該解析結果を提示してもよい。
【００２４】
［２−１．機能構成］
続いて、図４を参照して、本実施形態に係る情報処理装置４０の機能構成について説明する。図４は、本実施形態に係る情報処理装置４０の機能構成の一例を示す説明図である。図４に示したように、情報処理装置４０は、取得部４０２と、抽出部４０４と、調整部４０６と、分割部４０８と、解析部４１０と、出力部４１２と、を備える。なお、情報処理装置４０は、記憶素子を有してもよく、各種処理において用いられる情報を当該記憶素子に記憶させてもよい。例えば、情報処理装置４０は、取得部４０２により取得された情報、抽出部４０４による抽出の結果、調整部４０６により調整される前の方程式、調整部４０６により調整された後の方程式、分割部４０８により構築された解析モデル、解析部４１０による解析の結果等の情報を記憶素子に記憶させてもよい。当該記憶素子に記憶された情報は、各種処理において適宜参照され得る。
【００２５】
（取得部）
取得部４０２は、流動解析に用いられる情報を取得する機能を有する。例えば、流動解析に用いられる情報は、オペレータによって入力される。取得部４０２は、オペレータからの入力を受け付ける入力装置から送信される情報を受信することによって、流動解析に用いられる情報を取得し得る。また、情報処理装置４０は、当該入力装置と一体として構成されてもよい。その場合、取得部４０２の当該機能は、当該入力装置によって実現される。また、流動解析に用いられる情報は、予め情報処理装置４０の記憶素子に記憶されてもよい。その場合、取得部４０２は、当該記憶素子から流動解析に用いられる情報を取得する。
【００２６】
流動解析に用いられる情報は、例えば、キャビティの形状、流体の粘性、密度、熱伝導率等の物性値、及び流体の流入時並びに流出時における流体の速度や温度等の各種条件値に関する情報を含む。取得部４０２は、取得した情報を抽出部４０４及び分割部４０８へ出力する。
【００２７】
（抽出部）
抽出部４０４は、キャビティの形状を示す情報を取得部４０２から取得し、当該情報に基づいて、キャビティの厚さの分布を抽出し、抽出結果を調整部４０６へ出力する。
【００２８】
（調整部）
調整部４０６は、流体に関する物理量の保存則を表す方程式を、キャビティの厚さの分布に基づいて、調整する。具体的には、当該方程式は、流体の質量保存則を表す式及び流体の運動量保存則を表す式を含んでもよい。より具体的には、当該方程式は、流体の質量保存則を表す連続の式及び流体の運動量保存則を表すナビエストークスの運動方程式を含んでもよい。
【００２９】
調整部４０６は、調整後の方程式を示す情報を解析部４１０へ出力する。調整後の方程式は、解析部４１０による流動解析において、用いられる。調整部４０６は、例えば、キャビティの薄肉部についての流動解析における解析モデルとして、二次元的な解析モデルを適用可能とするように、方程式を調整する。本実施形態によれば、調整部４０６が行う方程式の調整によって、計算負荷の低減が実現される。なお、調整部４０６による方程式の調整の詳細については、後述する。
【００３０】
（分割部）
分割部４０８は、流体が流動するキャビティを複数の微小要素に分割する。分割部４０８は、キャビティの形状を示す情報に基づいて、キャビティの分割を行うことによって、流動解析における解析モデルを構築する。具体的には、分割部４０８は、キャビティの薄肉部について、厚さ方向についての要素分割が省略された二次元的な解析モデルを構築する。分割部４０８は、例えば、特許文献１に提案された技術を用いて、二次元的な解析モデルの構築を行ってもよい。分割部４０８は、構築した解析モデルを解析部４１０へ出力する。
【００３１】
（解析部）
解析部４１０は、調整部４０６により調整された方程式及び分割部４０８によって構築された解析モデルに基づいて、キャビティにおける流体の流動解析を行い、解析結果を出力部４１２へ出力する。具体的には、解析部４１０は、調整部４０６により調整された方程式を用いて、当該解析モデルにおける微小要素の各々について、流速を算出する。より具体的には、解析部４１０は、当該方程式を微小要素の各々について、解くことによって、キャビティにおける流速分布を算出する。また、解析部４１０は、当該解析モデルにおける微小要素の各々について、さらに圧力を算出してもよい。具体的には、解析部４１０は、当該方程式を微小要素の各々について、解くことによって、キャビティにおける圧力分布を算出する。
【００３２】
また、解析部４１０は、キャビティにおける流体の各時刻についての流動挙動を解析してもよい。この場合、解析部４１０は、例えば、解析対象の時刻のうち最先の時刻についての流動挙動を、取得部４０２により取得された情報に基づいて、解析した後、当該解析の結果を用いて直後の時刻についての流動挙動の解析を行う。また、調整部４０６は、解析部４１０に用いられる方程式を、解析対象の時刻の各々について、互いに異ならせるように調整してもよい。その場合、解析部４１０は、解析対象の時刻の各々について、互いに異なる方程式を用いてもよい。
【００３３】
（出力部）
出力部４１２は、解析部４１０により得られた流動解析の結果を出力する機能を有する。出力部４１２は、流動解析の結果を出力することによって、オペレータへ解析結果を提示する。例えば、出力部４１２は、画像を表示する表示装置へ解析結果を送信することによって、当該表示装置に解析結果を表示させてもよい。それにより、解析結果がオペレータへ提示される。また、情報処理装置４０は、当該表示装置と一体として構成されてもよい。その場合、出力部４１２の当該機能は、当該表示装置によって実現される。
【００３４】
出力部４１２は、媒体への印刷を行う印刷装置へ解析結果を送信することによって、当該印刷装置に解析結果を印刷させてもよい。それにより、解析結果がオペレータへ提示される。また、情報処理装置４０は、当該印刷装置と一体として構成されてもよい。その場合、出力部４１２の当該機能は、当該印刷装置によって実現される。
【００３５】
なお、出力部４１２により流動解析の結果が他の情報処理装置へ送信され、当該他の情報処理装置によって当該結果が二次的な処理に利用されてもよい。換言すると、出力部４１２は、必ずしもオペレータへ解析結果を提示しなくてもよい。
【００３６】
［２−２．方程式の調整］
続いて、図５〜図９を参照して、調整部４０６による方程式の調整の詳細について、説明する。以下、調整部４０６により調整される流体に関する物理量の保存則を表す方程式として、流体の質量保存則を表す連続の式、及び流体の運動量保存則を表すナビエストークスの運動方程式が用いられる例について説明する。また、以下の説明において、流体の流動方向をｘ軸方向とし、厚さ方向をｙ軸方向とし、ｘ軸及びｙ軸に直交する方向をｚ軸方向（以下、奥行き方向とも呼ぶ。）とする。
【００３７】
連続の式及びｘ軸方向、ｙ軸方向、並びにｚ軸方向についてのナビエストークスの運動方程式は、例えば、それぞれ下記式（３）〜式（６）により表される。
【００３８】
【数１】
2017162269.tif 000003
【００３９】
なお、式（３）〜式（６）において、ρは流体の密度を表し、ｕはｘ軸方向の流速成分を表し、ｖはｙ軸方向の流速成分を表し、ｗはｚ軸方向の流速成分を表す。また式（４）〜式（６）の右辺は、流体中の流体塊の表面に働く力である表面力に相当し、σｘ、σｙ、及びσｚは表面に垂直な引張応力を表し、τｘｙ、τｘｚ、及びτｙｚは表面に沿うせん断応力を表す。
【００４０】
ここで、説明を容易にするために、キャビティにおける流体の流動挙動を流動方向であるｘ軸方向についての一方向の流れとし、ｘ座標で記述されると考えると、流動方向について、連続の式及びナビエストークスの運動方程式は、下記式（７）及び式（８）で表される。
【００４１】
【数２】
2017162269.tif 000004
【００４２】
以下では、式（７）及び式（８）により表される連続の式及びナビエストークスの運動方程式を用いて流動解析を行う例について主に説明する。
【００４３】
2017162269.tif 000005
【００４４】
【数３】
2017162269.tif 000006
【００４５】
2017162269.tif 000007
【００４６】
図５は、キャビティの各位置における実際の厚さ及び質量流束密度について説明するための説明図である。図５では、実際のキャビティにおける流体Ｒ１及び内部空間としてのキャビティを形成する実際の金型Ｒ５が示されている。図５に示したように、キャビティの段差部の近傍において、段差部の上流側であるｘ座標「−Δｘ」の位置と、段差部の下流側であるｘ座標「Δｘ」の位置とでは、キャビティの厚さが異なる。ここで、質量保存則によれば、所定の領域内の流体の質量が時間的に保存されると仮定すると、所定の領域に単位時間当たりに流入する質量である流入質量と、所定の領域から単位時間当たりに流出する質量である流出質量とは等しい。
【００４７】
図５に示すように、ｘ座標「−Δｘ」の位置において、流速と密度との積により表される質量流束密度及びキャビティの厚さがそれぞれρｕ（−Δｘ）及びｈ−Δｈである場合、ｘ座標「−Δｘ」から「Δｘ」までの間の領域についての流入質量は、奥行き方向の単位長さあたりρｕ（−Δｘ）×（ｈ−Δｈ）である。また、ｘ座標「Δｘ」の位置において質量流束密度及びキャビティの厚さがρｕ（Δｘ）及びｈである場合、ｘ座標「−Δｘ」から「Δｘ」までの間の領域についての流出質量は、奥行き方向の単位長さあたりρｕ（Δｘ）×ｈである。よって、下記式（９）が成立する。
【００４８】
【数４】
2017162269.tif 000008
【００４９】
また、式（９）を整理すると、下記式（１０）が導出される。
【００５０】
【数５】
2017162269.tif 000009
【００５１】
式（１０）は、ｘ座標「Δｘ」の位置における質量流束密度ρｕ（Δｘ）がｘ座標「−Δｘ」の位置における質量流束密度ρｕ（−Δｘ）の（ｈ−Δｈ）／ｈ倍であることを示す。
【００５２】
図６に概念的に示すように、本実施形態に係る分割部４０８により構築される二次元的な解析モデルＭ４０では、厚さ方向についての要素分割が省略される。このように、解析モデルＭ４０では、厚さ方向の次元が省略されるので、互いに隣接する各微小要素間について、例えば、式（７）及び式（８）により表される連続の式及びナビエストークスの運動方程式が用いられる。ここで、ｘ座標「−Δｘ」の位置とｘ座標「Δｘ」の位置との間について、調整前の式（７）により表される連続の式を、時間微分項を０として、離散化して適用すると、下記式（１１）が導出される。
【００５３】
【数６】
2017162269.tif 000010
【００５４】
また、式（１１）を整理すると、下記式（１２）が導出される。
【００５５】
【数７】
2017162269.tif 000011
【００５６】
このように、調整前の式（７）により表される連続の式を用いて導出された式（１２）では、式（１０）と異なり、段差部の上流側と下流側との質量流束密度の関係性に対してキャビティの厚さの各位置における差異が与える影響が表されていない。ゆえに、厚さ方向の次元を省略した二次元的な解析モデルについて、調整前の式（７）により表される連続の式を直接的に用いて流動解析を行った場合、実際のキャビティの厚さ分布が流動挙動へ与える影響を反映した計算を実行することが困難となる場合がある。
【００５７】
2017162269.tif 000012
【００５８】
【数８】
2017162269.tif 000013
【００５９】
また、式（１３）を整理すると、式（１０）が導出される。ゆえに、調整後の連続の式を用いて導出された式（１３）では、式（１０）と同様に、段差部の上流側と下流側との質量流束密度の関係性に対してキャビティの厚さの各位置における差異が与える影響が表されている。
【００６０】
このように、厚さ方向の次元を省略した二次元的な解析モデルについて、調整部４０６による調整後の連続の式を用いて流動解析を行った場合、段差部の上流側と下流側との質量流束密度の関係性に対してキャビティの厚さの各位置における差異が与える影響を反映した計算を実行することができる。
【００６１】
また、運動量保存則によれば、所定の領域に単位時間当たりに流入する運動量である流入運動量と、所定の領域から単位時間当たりに流出する運動量である流出運動量との差は、所定の領域に働く力と等しい。以下では、所定の領域に働く力として、表面力を適用した場合について、説明する。
【００６２】
ｘ座標「−Δｘ」の位置において質量流束密度と流速との積により表される運動量流束密度がρｕ２（−Δｘ）である場合、ｘ座標「−Δｘ」から「Δｘ」までの間の領域についてのｘ軸方向の流入運動量は、奥行き方向の単位長さあたりρｕ２（−Δｘ）×（ｈ−Δｈ）である。また、ｘ座標「Δｘ」の位置において運動量流束密度がρｕ２（Δｘ）である場合、ｘ座標「−Δｘ」から「Δｘ」までの間の領域についてのｘ軸方向の流出運動量は、奥行き方向の単位長さあたりρｕ２（Δｘ）×ｈである。
【００６３】
ｘ座標「−Δｘ」の位置におけるｘ方向の引張応力がσｘ（−Δｘ）である場合、ｘ座標「−Δｘ」から「Δｘ」までの間の領域の上流側の面に働くｘ軸方向の表面力は、奥行き方向の単位長さあたりσｘ（−Δｘ）×（ｈ−Δｈ）である。また、ｘ座標「Δｘ」の位置におけるｘ方向の引張応力がσｘ（Δｘ）である場合、ｘ座標「−Δｘ」から「Δｘ」までの間の領域の上流側の面に働くｘ軸方向の表面力は、奥行き方向の単位長さあたり−σｘ（Δｘ）×ｈである。よって、下記式（１４）が成立する。
【００６４】
【数９】
2017162269.tif 000014
【００６５】
ここで、ｘ座標「−Δｘ」の位置とｘ座標「Δｘ」の位置との間について、調整前の式（８）により表されるナビエストークスの運動方程式を、時間微分項を０として、離散化して適用すると、下記式（１５）が導出される。
【００６６】
【数１０】
2017162269.tif 000015
【００６７】
このように、調整前の式（８）により表されるナビエストークスの運動方程式を用いて導出された式（１５）では、式（１４）と異なり、段差部の上流側と下流側との運動量流束密度及び表面力の関係性に対してキャビティの厚さの各位置における差異が与える影響が表されていない。ゆえに、厚さ方向の次元を省略した二次元的な解析モデルについて、調整前の式（８）により表されるナビエストークスの運動方程式を直接的に用いて流動解析を行った場合、実際のキャビティの厚さ分布が流動挙動へ与える影響を反映した計算を実行することが困難となる場合がある。
【００６８】
2017162269.tif 000016
【００６９】
【数１１】
2017162269.tif 000017
【００７０】
また、式（１６）を整理すると、式（１４）が導出される。ゆえに、調整後のナビエストークスの運動方程式を用いて導出された式（１６）では、式（１４）と同様に、段差部の上流側と下流側との運動量流束密度及び表面力の関係性に対してキャビティの厚さの各位置における差異が与える影響が表されている。
【００７１】
このように、厚さ方向の次元を省略した二次元的な解析モデルについて、調整部４０６による調整後のナビエストークスの運動方程式を用いて流動解析を行った場合、段差部の上流側と下流側との運動量流束密度及び表面力に対してキャビティの厚さの各位置における差異が与える影響を反映した計算を実行することができる。
【００７２】
以上説明したように、本実施形態に係る調整部４０６による方程式の調整によれば、キャビティの薄肉部の各位置における厚さが均一でない場合であっても、厚さ方向の次元を省略した二次元的な解析モデルを用いて、解析精度を確保しつつ、流動解析を行うことができる。ゆえに、キャビティにおける流体の流動解析において、計算負荷をより低減させることが可能である。
【００７３】
（湧き出しに関する項の決定）
調整部４０６は、流体の湧き出しに関する項を決定し、流動解析において用いられる連続の式、及びナビエストークスの運動方程式へ当該項を追加する調整を行ってもよい。例えば、調整部４０６は、各微小要素に対応する位置における湧き出しによる流体の質量変化量に相当する項Ｓ及び、湧き出しにより流体へ掛かるｘ軸方向の力成分に相当する項Ｓｘを決定し、各方程式へ追加する。具体的には、当該項Ｓ及び項Ｓｘの追加によって、式（７）及び式（８）により表される連続の式及びナビエストークスの運動方程式は、下記式（１７）及び式（１８）によって表される式へ調整される。
【００７４】
【数１２】
2017162269.tif 000018
【００７５】
このように、情報処理装置４０が行う流動解析において用いられる連続の式、及びナビエストークスの運動方程式は、流体の湧き出しに関する項を含んでもよい。
【００７６】
（重力に関する項の決定）
調整部４０６は、流体へ掛かる重力に関する項を決定し、流動解析において用いられるナビエストークスの運動方程式へ当該項を追加する調整を行ってもよい。例えば、調整部４０６は、流体へ掛かる重力のｘ方向成分に相当する項を決定し、ナビエストークスの運動方程式へ追加する。具体的には、当該項の追加によって、式（８）により表されるナビエストークスの運動方程式は、下記式（１９）によって表される式へ調整される。なお、式（１９）において、Ｋｘは重力加速度のｘ方向成分を表す。
【００７７】
【数１３】
2017162269.tif 000019
【００７８】
このように、情報処理装置４０が行う流動解析において用いられるナビエストークスの運動方程式は、流体へ掛かる重力に関する項を含んでもよい。なお、情報処理装置４０が行う流動解析において用いられるナビエストークスの運動方程式は、流体へ掛かる重力以外の質量力を含んでもよい。質量力は、質量や電荷を支える物質そのものに作用する力であり、例えば、電気力やローレンツ力を含む。
【００７９】
（流動抵抗に関する項の決定）
調整部４０６は、流体の流動抵抗に関する項を決定し、流動解析において用いられるナビエストークスの運動方程式へ当該項を追加する調整を行ってもよい。例えば、調整部４０６は、各微小要素に対応する位置における流体へ掛かる流動抵抗のｘ方向成分に相当する項を決定し、ナビエストークスの運動方程式へ追加する。具体的には、当該項の追加によって、式（８）により表されるナビエストークスの運動方程式は、下記式（２０）によって表される式へ調整される。なお、式（２０）において、キャビティの厚さｈの２乗にｘ軸方向の流速ｕを乗じて１２で除し、−１を乗じた項が、流動抵抗のｘ方向成分に相当する項に該当する。
【００８０】
【数１４】
2017162269.tif 000020
【００８１】
このように、情報処理装置４０が行う流動解析において用いられるナビエストークスの運動方程式は、流体の流動抵抗に関する項を含んでもよい。なお、当該項は、流体の流動抵抗に相当する項であればよく、式（２０）により例示した項に特に限定されない。例えば、調整部４０６は、流動抵抗を算出するための数式として経験的に得られた数式に基づいて、各微小要素に対応する位置における流体へ掛かる流動抵抗のｘ方向成分に相当する項を決定してもよい。
【００８２】
（表面張力の決定）
調整部４０６は、流体の表面張力に関する項を決定し、流動解析において用いられるナビエストークスの運動方程式へ当該項を追加する調整を行ってもよい。例えば、調整部４０６は、各微小要素に対応する位置における流体へ掛かる表面張力のｘ方向成分に相当する項を決定し、ナビエストークスの運動方程式へ追加する。具体的には、当該項の追加によって、式（８）により表されるナビエストークスの運動方程式は、下記式（２１）によって表される式へ調整される。なお、式（２１）において、Ｐｓは流体へ掛かる表面張力のｘ方向成分を示し、関数ｆ（ｘ）は微小要素内に金型表面と流体の自由表面とが接する箇所が存在するか否かを判定するための関数を示す。
【００８３】
【数１５】
2017162269.tif 000021
【００８４】
図７は、キャビティにおける流体に掛かる表面張力について説明するための説明図である。図７に示したように、空間としてのキャビティを形成する金型に対して濡れ広がった流体の端部において、流体の自由表面は、金型表面に対して接触角θを有する曲面形状を有し得る。この場合、流体を挟む２つの金型表面のそれぞれと流体の自由表面とが接する箇所において、接触角θに応じた方向に表面張力が生じる。ゆえに、調整部４０６は、下記式（２２）で表されるＰｓを、流体へ掛かる表面張力のｘ方向成分として算出してもよい。なお、式（２２）において、σは流体へ掛かる単位長さあたりの表面張力の絶対値を示す。
【００８５】
【数１６】
2017162269.tif 000022
【００８６】
また、図７に示したように、関数ｆ（ｘ）は、例えば、金型表面と流体の自由表面とが接する位置では１を返し、金型表面と流体の自由表面とが接しない位置では０を返す関数である。ゆえに、式（２１）によれば、微小要素内に金型表面と流体の自由表面とが接する箇所が存在する場合には、Ｐｓに１が乗じられ、微小要素内に金型表面と流体の自由表面とが接する箇所が存在しない場合には、Ｐｓに０が乗じられる。
【００８７】
このように、情報処理装置４０が行う流動解析において用いられるナビエストークスの運動方程式は、流体の表面張力に関する項を含んでもよい。
【００８８】
2017162269.tif 000023
【００８９】
【数１７】
2017162269.tif 000024
【００９０】
上述したように、分割部４０８により構築される二次元的な解析モデルＭ４０では、厚さ方向の次元が省略されるので、各微小要素Ｂ３の厚さは同一の値に設定される。例えば、図８に示したように、実際のキャビティの厚さｈに対して、解析モデルＭ４０における微小要素Ｂ３の厚さはｈ０に設定される。ここで、実際のキャビティの各位置における厚さは均一でない場合がある。一方、解析モデルＭ４０における微小要素Ｂ３の厚さは、キャビティの位置によらず、ｈ０に設定される。
【００９１】
ゆえに、流動解析において用いられる方程式におけるキャビティの厚さと相関を有する項をキャビティの厚さに基づいて正規化することによって、各方程式を二次元的な解析モデルに対して、より精度良く適合することができる。
【００９２】
なお、調整部４０６は、式（２０）において例示された流体の流動抵抗に関する項又は式（２１）において例示された流体の表面張力に関する項についても、同様に正規化を行ってもよい。
【００９３】
2017162269.tif 000025
【００９４】
【数１８】
2017162269.tif 000026
【００９５】
2017162269.tif 000027
【００９６】
実際のキャビティの厚さは、時間に応じて変化する場合がある。例えば、図９に示したように、時刻「−Δｔ」においてｈ＋Δｈであったキャビティの所定の位置における厚さが、時間経過後の時刻「Δｔ」においてｈへ減少する場合がある。時刻「−Δｔ」における密度がρ（−Δｔ）である場合、時刻「−Δｔ」における当該所定の位置におけるｘｙ平面の単位面積当たりの流体の質量は、ρ（−Δｔ）×（ｈ＋Δｈ）である。また、時刻「Δｔ」における密度がρ（Δｔ）である場合、時刻「Δｔ」における当該所定の位置におけるｘｙ平面の単位面積当たりの流体の質量は、ρ（Δｔ）×ｈである。ここで、上述した空間微分演算子の置換についての説明において述べたように所定の領域内の流体の質量が時間的に保存されると仮定すると、下記式（２５）が成立する。
【００９７】
【数１９】
2017162269.tif 000028
【００９８】
上述したように、分割部４０８により構築される二次元的な解析モデルＭ４０では、上述したように、例えば、式（７）及び式（８）により表される連続の式及びナビエストークスの運動方程式が用いられる。ここで、図９に示したキャビティの所定の位置において、時刻「−Δｔ」から時刻「Δｔ」の間について、調整前の式（７）により表される連続の式を、空間微分項を０として、離散化して適用すると、下記式（２６）が導出される。
【００９９】
【数２０】
2017162269.tif 000029
【０１００】
また、式（２６）を整理すると、下記式（２７）が導出される。
【０１０１】
【数２１】
2017162269.tif 000030
【０１０２】
このように、調整前の式（７）により表される連続の式を用いて導出された式（２７）では、式（２５）と異なり、流体の密度の時間変化に対してキャビティの厚さの時間変化が与える影響が表されていない。ゆえに、厚さ方向の次元を省略した二次元的な解析モデルについて、調整前の式（７）により表される連続の式を直接的に用いて流動解析を行った場合、実際のキャビティの厚さの時間変化が流動挙動へ与える影響を反映した計算を実行することが困難となる場合がある。
【０１０３】
2017162269.tif 000031
【０１０４】
【数２２】
2017162269.tif 000032
【０１０５】
また、式（２８）を整理すると、式（２５）が導出される。ゆえに、調整後の連続の式を用いて導出された式（２８）では、式（２５）と同様に、流体の密度の時間変化に対してキャビティの厚さの時間変化が与える影響が表されている。
【０１０６】
このように、厚さ方向の次元を省略した二次元的な解析モデルについて、調整部４０６による調整後の連続の式を用いて流動解析を行った場合、実際のキャビティの厚さの時間変化が流動挙動へ与える影響を反映した計算を実行することができる。
【０１０７】
調整部４０６は、このような調整部４０６による時間微分演算子の置換を、同様に、式（８）により表されるナビエストークスの運動方程式について行ってもよい。それにより、厚さ方向の次元を省略した二次元的な解析モデルを用いた流動解析において、調整後のナビエストークスの運動方程式を適用することによって、実際のキャビティの厚さの時間変化が流動挙動へ与える影響を反映した計算を実行することができる。
【０１０８】
本実施形態と異なり、厚さ方向に複数の微小要素の分割を行った解析モデルを用いた流動解析では、実際のキャビティの厚さが時間に応じて変化する場合、厚さの時間変化に伴って、微小要素の厚さを変化させる場合がある。そのような場合、各微小要素のアスペクト比が過剰に大きくなり、解析精度を確保することが困難となり得る。アスペクト比の増大に伴う解析精度の低下を抑制するために、より細かい微小要素に分割する場合がある。それにより、計算負担が増大し得る。
【０１０９】
一方、本実施形態では、実際のキャビティの厚さが時間に応じて変化する場合であっても、厚さ方向の次元を省略した二次元的な解析モデルを用いて、解析精度を確保しつつ、アスペクト比を一定に保ちながら流動解析を行うことができる。ゆえに、キャビティにおける流体の流動解析において、計算負荷をより低減させることが可能である。
【０１１０】
＜３．動作＞
続いて、図１０及び図１１を参照して、本実施形態に係る情報処理装置４０が行う各種処理の流れについて説明する。
【０１１１】
図１０は、本実施形態に係る情報処理装置４０が行う流動解析処理の流れの一例を示すフローチャートである。図１０に示したように、まず、取得部４０２は、流動解析に用いられる情報を取得し（ステップＳ５０２）、抽出部４０４及び分割部４０８へ出力する。次に、抽出部４０４は、キャビティの厚さの分布を抽出し（ステップＳ５０４）、抽出結果を調整部４０６へ出力する。そして、調整部４０６は、流体に関する物理量の保存則を表す方程式を、キャビティの厚さの分布に基づいて、調整し（ステップＳ６００）、調整後の方程式を示す情報を解析部４１０へ出力する。具体的には、当該方程式は、流体の質量保存則を表す連続の式及び流体の運動量保存則を表すナビエストークスの運動方程式である。
【０１１２】
次に、分割部４０８は、流体が流動するキャビティを複数の微小要素に分割する（ステップＳ５０６）。また、分割部４０８は、キャビティの分割を行うことによって構築された解析モデルを解析部４１０へ出力する。そして、解析部４１０は、調整部４０６により調整された方程式を用いて、当該解析モデルにおける微小要素の各々について、流速及び圧力を算出し（ステップＳ５０８）、解析結果を出力部４１２へ出力する。次に、出力部４１２は、解析部４１０により得られた流動解析の結果を出力し（ステップＳ５１０）、図１０に示した処理は終了する。
【０１１３】
続いて、図１１を参照して、調整部４０６が行う方程式調整処理（図１０に示したステップＳ６００）について、より詳細に説明する。
【０１１４】
2017162269.tif 000033
【０１１５】
＜４．ハードウェア構成＞
以上、本開示の実施形態を説明した。上述した情報処理装置４０の処理は、ソフトウェアと、以下に説明する情報処理装置４０のハードウェアとの協働により実現される。
【０１１６】
図１２は、本開示に係る情報処理装置４０のハードウェア構成の一例を示す説明図である。図１２に示したように、情報処理装置４０は、ＣＰＵ（Ｃｅｎｔｒａｌ  Ｐｒｏｃｅｓｓｉｎｇ  Ｕｎｉｔ）１４２と、ＲＯＭ（Ｒｅａｄ  Ｏｎｌｙ  Ｍｅｍｏｒｙ）１４４と、ＲＡＭ（Ｒａｎｄｏｍ  Ａｃｃｅｓｓ  Ｍｅｍｏｒｙ）１４６と、ブリッジ１４８と、バス１５０と、インタフェース１５２と、入力装置１５４と、出力装置１５６と、ストレージ装置１５８と、ドライブ１６０と、接続ポート１６２と、通信装置１６４とを備える。
【０１１７】
ＣＰＵ１４２は、演算処理装置および制御装置として機能し、各種プログラムと協働して情報処理装置４０内の各機能構成の動作を実現する。また、ＣＰＵ１４２は、マイクロプロセッサであってもよい。ＲＯＭ１４４は、ＣＰＵ１４２が使用するプログラムまたは演算パラメータ等を記憶する。ＲＡＭ１４６は、ＣＰＵ１４２の実行において使用するプログラムまたは実行において適宜変化するパラメータ等を一時記憶する。ＣＰＵ１４２、ＲＯＭ１４４およびＲＡＭ１４６は、ＣＰＵバスなどから構成される内部バスにより相互に接続されている。
【０１１８】
入力装置１５４は、例えば、マウス、キーボード、タッチパネル、ボタン、マイクロフォン、スイッチおよびレバーなどオペレータが情報を入力するための入力手段であり、オペレータによる入力に基づいて入力信号を生成し、ＣＰＵ１４２に出力する入力制御回路などから構成されている。情報処理装置４０のオペレータは、入力装置１５４を操作することにより、情報処理装置４０に対して各種のデータを入力したり処理動作を指示したりすることができる。
【０１１９】
出力装置１５６は、例えば、液晶ディスプレイ（ＬＣＤ）装置、ＯＬＥＤ（Ｏｒｇａｎｉｃ  Ｌｉｇｈｔ  Ｅｍｉｔｔｉｎｇ  Ｄｉｏｄｅ）装置、ランプなどの装置への出力を行う。さらに、出力装置１５６は、スピーカおよびヘッドフォンなどの音声出力を行ってもよい。
【０１２０】
ストレージ装置１５８は、データ格納用の装置である。ストレージ装置１５８は、記憶媒体、記憶媒体にデータを記録する記録装置、記憶媒体からデータを読み出す読出し装置および記憶媒体に記録されたデータを削除する削除装置等を含んでもよい。ストレージ装置１５８は、ＣＰＵ１４２が実行するプログラムや各種データを格納する。
【０１２１】
ドライブ１６０は、記憶媒体用リーダライタであり、情報処理装置４０に内蔵、あるいは外付けされる。ドライブ１６０は、装着されている磁気ディスク、光ディスク、光磁気ディスク、または半導体メモリ等のリムーバブル記憶媒体に記録されている情報を読み出して、ＲＡＭ１４４に出力する。また、ドライブ１６０は、リムーバブル記憶媒体に情報を書込むこともできる。
【０１２２】
接続ポート１６２は、例えば、情報処理装置４０の外部の情報処理装置または周辺機器と接続するためのバスである。また、接続ポート１６２は、ＵＳＢ（Ｕｎｉｖｅｒｓａｌ  Ｓｅｒｉａｌ  Ｂｕｓ）であってもよい。
【０１２３】
通信装置１６４は、例えば、ネットワークに接続するための通信デバイスで構成された通信インタフェースである。また、通信装置１６４は、赤外線通信対応装置であっても、無線ＬＡＮ（Ｌｏｃａｌ  Ａｒｅａ  Ｎｅｔｗｏｒｋ）対応通信装置であっても、ＬＴＥ（Ｌｏｎｇ  Ｔｅｒｍ  Ｅｖｏｌｕｔｉｏｎ）対応通信装置であっても、有線による通信を行うワイヤー通信装置であってもよい。
【０１２４】
なお、上述のような本実施形態に係る情報処理装置４０の各機能を実現するためのコンピュータプログラムを作製し、ＰＣ等に実装することが可能である。本実施形態に係る情報処理装置４０は、本開示に係るコンピュータに相当し得る。また、このようなコンピュータプログラムが格納された、コンピュータで読み取り可能な記録媒体も提供することができる。記録媒体は、例えば、磁気ディスク、光ディスク、光磁気ディスク、フラッシュメモリ等である。また、上記のコンピュータプログラムは、記録媒体を用いずに、例えばネットワークを介して配信されてもよい。また、本実施形態に係る情報処理装置４０の各機能は複数のコンピュータにより分割されてもよく、その場合、当該複数のコンピュータが有する各機能は、上記のコンピュータプログラムにより実現され得る。
【０１２５】
＜５．まとめ＞
以上説明したように、本開示の実施形態によれば、調整部４０６は、流体に関する物理量の保存則を表す方程式を、キャビティの厚さの分布に基づいて、調整する。例えば、調整部４０６は、キャビティの厚さの分布に基づく調整として、上述した空間微分演算子の置換及び時間微分演算子の置換の少なくとも一方を行う。それにより、厚さ方向の次元を省略した二次元的な解析モデルを用いて、解析精度を確保しつつ、流動解析を行うことができる。ゆえに、キャビティにおける流体の流動解析において、計算負荷をより低減させることが可能である。
【０１２６】
なお、上述した実施形態に係る情報処理装置は、様々な電子機器に利用することができる。具体的には、パーソナルコンピュータやスマートフォン等の端末の一部として利用されてもよい。あるいは、射出成形装置や生産工程における検査装置等の産業用機器の一部として利用されてもよい。本実施形態に係る電子機器によっても、上述した実施形態に係る情報処理装置と同様の効果を得ることが可能である。
【０１２７】
なお、本明細書において説明した各装置による一連の制御処理は、ソフトウェア、ハードウェア、及びソフトウェアとハードウェアとの組合せのいずれを用いて実現されてもよい。ソフトウェアを構成するプログラムは、例えば、各装置の内部又は外部に設けられる記憶媒体（非一時的な媒体：non-transitory  media）に予め格納される。そして、各プログラムは、例えば、実行時にＲＡＭに読み込まれ、ＣＰＵなどのプロセッサにより実行される。各プログラムを実行するプロセッサは、単数であっても複数であってもよい。
【０１２８】
また、本明細書においてフローチャートを用いて説明した処理は、必ずしもフローチャートに示された順序で実行されなくてもよい。いくつかの処理ステップは、並列的に実行されてもよい。また、追加的な処理ステップが採用されてもよく、一部の処理ステップが省略されてもよい。
【０１２９】
以上、添付図面を参照しながら本開示の好適な実施形態について詳細に説明したが、本開示の技術的範囲は係る例に限定されない。本開示の技術分野における通常の知識を有する者であれば、特許請求の範囲に記載された技術的思想の範疇内において、各種の変更例または修正例に想到し得ることは明らかであり、これらについても、当然に本開示の技術的範囲に属するものと了解される。
【０１３０】
また、本明細書に記載された効果は、あくまで説明的または例示的なものであって限定的ではない。つまり、本開示に係る技術は、上記の効果とともに、または上記の効果に代えて、本明細書の記載から当業者には明らかな他の効果を奏しうる。
【０１３１】
なお、以下のような構成も本開示の技術的範囲に属する。
（１）
流体が流動するキャビティを複数の微小要素に分割する分割部と、
前記流体に関する物理量の保存則を表す方程式を、前記キャビティの厚さの分布に基づいて、調整する調整部と、
前記調整部により調整された前記方程式を用いて、前記微小要素の各々について、流速を算出する解析部と、
を備える、情報処理装置。
（２）
前記調整部は、前記方程式の空間微分項における被演算子に前記厚さの空間微分及び前記厚さの逆数を乗じて得られる項の前記方程式への付加に相当する調整を行う、前記（１）に記載の情報処理装置。
（３）
前記調整部は、前記方程式の時間微分項における被演算子に前記厚さの時間微分及び前記厚さの逆数を乗じて得られる項の前記方程式への付加に相当する調整を行う、前記（１）又は（２）に記載の情報処理装置。
2017162269.tif 000034
【数２３】
2017162269.tif 000035
2017162269.tif 000036
【数２４】
2017162269.tif 000037
（６）
前記方程式は、前記流体の質量保存則を表す連続の式、及び前記流体の運動量保存則を表すナビエストークスの運動方程式を含む、前記（１）〜（５）のいずれか一項に記載の情報処理装置。
（７）
前記連続の式、及び前記ナビエストークスの運動方程式は、前記流体の湧き出しに関する項を含む、前記（６）に記載の情報処理装置。
（８）
前記ナビエストークスの運動方程式は、前記流体へ掛かる重力に関する項を含む、前記（６）又は（７）に記載の情報処理装置。
（９）
前記ナビエストークスの運動方程式は、前記流体の流動抵抗に関する項を含む、前記（６）〜（８）のいずれか一項に記載の情報処理装置。
（１０）
前記ナビエストークスの運動方程式は、前記流体の表面張力に関する項を含む、前記（６）〜（９）のいずれか一項に記載の情報処理装置。
（１１）
前記調整部は、前記方程式における前記厚さと相関を有する項に前記厚さの逆数を乗じることによって、当該項を正規化する調整を行う、前記（１）〜（１０）のいずれか一項に記載の情報処理装置。
（１２）
前記解析部は、前記微小要素の各々について、圧力を算出する、前記（１）〜（１１）のいずれか一項に記載の情報処理装置。
（１３）
前記（１）〜（１２）のいずれか一項に記載の情報処理装置を備える、電子機器。
（１４）
流体が流動するキャビティを複数の微小要素に分割することと、
前記流体に関する物理量の保存則を表す方程式を、前記キャビティの厚さの分布に基づいて、情報処理装置によって、調整することと、
調整された前記方程式を用いて、前記微小要素の各々について、流速を算出することと、
を含む情報処理方法。
（１５）
コンピュータを、
流体が流動するキャビティを複数の微小要素に分割する分割部と、
前記流体に関する物理量の保存則を表す方程式を、前記キャビティの厚さの分布に基づいて、調整する調整部と、
前記調整部により調整された前記方程式を用いて、前記微小要素の各々について、流速を算出する解析部と、
として機能させるためのプログラム。
【符号の説明】
【０１３２】
１０  キャビティ
４０  情報処理装置
１２０  薄肉部
１２２  第１の領域
１２４  第２の領域
１４２  （Ｃｅｎｔｒａｌ  Ｐｒｏｃｅｓｓｉｎｇ  Ｕｎｉｔ）
１４４  （Ｒｅａｄ  Ｏｎｌｙ  Ｍｅｍｏｒｙ）
１４６  （Ｒａｎｄｏｍ  Ａｃｃｅｓｓ  Ｍｅｍｏｒｙ）
１４８  ブリッジ
１５０  バス
１５２  インタフェース
１５４  入力装置
１５６  出力装置
１５８  ストレージ装置
１６０  ドライブ
１６２  接続ポート
１６４  通信装置
４０２  取得部
４０４  抽出部
４０６  調整部
４０８  分割部
４１０  解析部
【図１】
2017162269.tif 000042
【図２】
2017162269.tif 000043
【図３】
2017162269.tif 000044
【図４】
2017162269.tif 000045
【図５】
2017162269.tif 000046
【図６】
2017162269.tif 000047
【図７】
2017162269.tif 000048
【図８】
2017162269.tif 000049
【図９】
2017162269.tif 000050
【図１０】
2017162269.tif 000051
【図１１】
2017162269.tif 000052
【図１２】
2017162269.tif 000053