(57)【要約】      （修正有）【課題】マルチセルシステムダイナミクスモデルの自動化された格子形成用の方法及び装置を提供する。【解決手段】隣接する一列の多角形によってシステム表面を表すステップ、システム全体を取り囲む単一のモデルセルを創出するステップと、すべての表面多角形を単一のモデルセルに割り当てるステップ、１つ以上の細分面を使用して、割り当てられたすべての交差多角形を細分することを含む、表面多角形列によって取り囲まれた単一のモデルセルの部分を少なくとも細分するステップ、細分ステップから得られる細分された個々のセルに表面多角形を再び割り当てるステップ、共通のシステム表面から、２つ以上の表面多角形を有するそれらの細分されたセルを識別するステップ、それらのセルの有効な共通境界領域及び法線ベクトルを計算するステップ、「分割」セル及び「スリバ」セルを取り扱うステップを含む。【選択図】図１
2007179548.tif 000002
【特許請求の範囲】
【請求項１】
１つ以上の表面を有する流体ダイナミクスシステムのマルチセルコンピュータモデル用の固定幾何学的格子を形成するための方法であって、
一列の隣接多角形によってシステム表面を表すステップと、
システム全体を取り囲む単一のモデルセルを創出するステップと、
すべての表面多角形を単一のモデルセルに割り当てるステップと、
１つ以上の細分面を使用して、割り当てられたすべての交差多角形を細分することを含む、表面多角形列によって取り囲まれた単一のモデルセルの部分を少なくとも細分するステップと、
細分ステップから得られる細分された個々のセルに表面多角形のすべてを再び割り当てるステップと、
共通のシステム表面から、２つ以上の表面多角形を有するそれらの細分されたセルを識別するステップと、
それらの識別されたセルの表面多角形に関し有効な境界領域及び法線ベクトルを計算するステップと、
を含む方法。
【請求項２】
計算するステップが、それぞれの識別されたセルについて、
個々の表面多角形の領域に等しい大きさ有する共通の表面多角形のそれぞれから、法線ベクトルを計算するステップと、
単一の法線ベクトルの方向の個々の共通の表面多角形の投影面積に等しい大きさを有する単一の法線ベクトルに、法線ベクトルを結合するステップと、
を含む、請求項１に記載の方法。
【請求項３】
１つ以上の細分されたセルのそれぞれが、流体システム表面によって２つ以上の別個の部分に分割され、さらに、方法が、分割セルのそれぞれのセル部分と個々の隣接境界セル及び非境界セルとを接続する情報を決定するステップを含み、
分割セル部分のそれぞれに同一の格子位置情報が付与され、また
分割セル部分のそれぞれが、その個々の部分に関する輸送方程式を解く際に無関係に取り扱われる、請求項１に記載の方法。
【請求項４】
予め選択された量部分未満の容積を有する境界セルを識別するステップと、
それぞれの識別された部分境界セルと隣接セルとを対にするステップと、
対にしたそれぞれの隣接セルの幾何学的情報を調整して、個々の部分境界セルの容積情報、面情報、及び境界情報を少なくとも含めるステップとを含み、
形成された格子の使用中に、それぞれのセル対の調整された隣接セルのみへの、および該隣接セルのみからの熱物理学的な値フラックスを規定する輸送方程式が解かれる、請求項１に記載の方法。
【請求項５】
１つ以上の表面を有する流体ダイナミクスシステムのマルチセルコンピュータモデル用の固定幾何学的格子を自動的に形成するための装置であって、
デジタルコンピュータと、
コンピュータに記憶された計算流体ダイナミクスプログラムと、を含み、
計算流体ダイナミクスプログラムが、自動的に、
一列の隣接多角形によってシステム表面を表すデータセットを受け取り、
システム全体を取り囲む単一のモデルセルを創出し、
すべての表面多角形を単一のモデルセルに割り当て、
１つ以上の細分面を使用して、割り当てられたすべての交差多角形を細分することを含めて、表面多角形列によって取り囲まれた単一のモデルセルの部分を少なくとも細分し、
細分ステップから得られる細分された個々のセルに表面多角形のすべてを再び割り当て、
共通のシステム表面から、２つ以上の表面多角形を有するそれらの細分されたセルを識別し、また
それらの識別されたセルの表面多角形に関し有効な境界領域及び法線ベクトルを計算するためのソフトウェアを有するルーチンを含む、装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明の開示は、一般に、流体動的システムのような動的システムをコンピュータモデル化するための方法及び装置に関する。より詳しくは、本発明の開示は、マルチセルシステムダイナミクスモデルの自動化された格子形成用の方法、及び方法を実行するための装置に関する。
【背景技術】
【０００２】
コンピュータ、特に高速デジタルコンピュータを使用して、流体動的システムを含む動的システムをモデル化することは、実際のシステムを物理的に構成して、試験することを必要とすることなく、安定した熱物理学的及び過渡的な状態の両方についてシステム性能を予測する周知で費用効率的な方法である。コンピュータモデル化の利点は、システム構造及び構成における変化の性能上の効果を容易に算定でき、これによって、商業的プロトタイプを構成する前にシステム設計の最適化が行われることである。
【０００３】
公知のモデル化プログラムは、一般に、モデル化すべき構造が複数の別個の容積ユニット（セル）に分割される「マルチセル」方法を使用する。典型的に、コンピュータは、隣接セルへの又はそれからの及び／又はセル境界を通した熱物理学的な値ユニットの輸送を規定する保存方程式を解くことによって、それぞれのセル内のシステムの部分の熱物理学的な値、例えば質量、運動量、及びエネルギ値、ならびに密度、圧力、速度、及び温度のような関連の流体システム設計パラメータを計算するために使用される。デカルト座標を使用する固定幾何学的システムモデルでは、またシステム境界がない場合、それぞれのセルは、立方体形状のセルの６つの面に隣接して配置された６つのセル隣接部を有するであろう。計算流体ダイナミクス（「ＣＦＤ」）モデル化プログラムの例は、コンバージェントシンキング（Ｃｏｎｖｅｒｇｅｎｔ  Ｔｈｉｎｋｉｎｇ）有限会社（ウィスコンシン州マディソン（Ｍａｄｉｓｏｎ，Ｗｉｓｃｏｎｓｉｎ））から入手可能なＭｏＳＥＳプログラムである。しかし、改良は既存のモデル化プログラムで可能であり、かつ望ましい。
【０００４】
従来のＣＦＤプログラムでは、ユーザは、シミュレートすべき対象物を表す３次元格子を最初に供給しなければならない。これらのＣＦＤプログラムでは、格子又は「メッシュ」生成は、特に、それぞれの過渡的な時間ステップの後に格子を再形成する必要があり得る運動境界を有する複雑な幾何学的形状にとって、モデル化プロセスの最もユーザ集中的な部分である。
【０００５】
境界にフィットされた２種類の格子がある。より従来型の境界にフィットされた格子は、幾何学的形状の形に従うように境界の近くのセルをモーフィングし、例えば、境界の近くの六面のセルは必ずしも完全な立方体ではないであろう。他の方法は、共通に「カットセル」方法と呼ばれる。従来の「カットセル」方法では、境界セルは１回以上カットされて、モデル化される流体システムの固体表面とのより優れた「フィット」を提供する。固体表面は、固体表面にフィットされた一列の三角形によって表すことができ、次に、これらの三角形は、固体表面に隣接する個々のモデルセルと関連付けられる。従来のカットセル方法では、多数の表面三角形によってカットされるセルを処理する２つの方法がある。１つの方法では、それぞれの三角形の表面効果は、それぞれのカットセルについて別個に取り扱われる。このプロセスは、より長いＣＦＤシミュレーション実行時間をもたらす可能性がある。第２の方法では、多数の三角形のカットは単一の平坦なカットによって近似される。この方法では、情報が表面近似で失われるので、元の幾何学的形状は正確に表されない。
【０００６】
カットセルのデカルト方法は、２種類の問題のあるセル、すなわち「スリバ」セル及び「分割」セルを生じる可能性がある。スリバセルは、非境界セルと比較して非常に小さな容積を有する境界セルである。小さな容積を有するセルは、解を安定維持するために、過度に小さな時間ステップ又は過度の繰り返し数、したがって、より多くの計算を必要とする可能性がある。
【０００７】
分割セルはまた、論理的なブロック構造化格子にとって問題であり、その形態はデカルトである。デカルト格子では、それぞれのセルは立方体として形成され、その６つの面のそれぞれの方向に１つのみの隣接部を有すると考えられる。それぞれのセルは論理（ｘ、ｙ、ｚ）座標を有し、その隣接部のそれぞれは正確に１つの論理ステップだけ異なる。デカルト格子の単純さの多くは、それらの６つの隣接部へのセルのこの構造化接続の直接の結果である。このシステムが適切に動作するには、論理座標に１つのみのセルが存在しなければならない。モデル化すべき流体システムの表面が、デカルトセルを２つ以上の部分にカットする場合、このことは、同一の位置に２つの分離セルがあるので、従来のＣＦＤコードの問題を提起する。これらの２つの別個のセル部分は、同一の論理座標を有し、それゆえ、それらの個々の隣接部について、論理座標のみを観察することによって適切なセル部分を配置することを不可能にする。
【０００８】
【特許文献１】米国特許出願第１１／３１８，６３４号明細書
【特許文献２】米国特許出願第１１／３１８，６３３号明細書
【特許文献３】米国特許出願第１１／３１８，６３２号明細書
【発明の開示】
【発明が解決しようとする課題】
【０００９】
ある場合には、ユーザの介入によっても、より伝統的なカットセル方法に関わるこれらの問題の克服が不可能であることを認識すべきである。開示した方法及び装置は、上述の問題の１つ以上を緩和することに向けられる。
【課題を解決するための手段】
【００１０】
一形態では、本明細書に具体化されかつ広く記載されるように、１つ以上の表面を有する流体ダイナミクスシステムのマルチセルコンピュータモデル用の固定幾何学格子を形成するための方法は、隣接する一列の多角形によってシステム表面を表すステップを含む。本方法はまた、システム全体を取り囲む単一のモデルセルを創出するステップと、すべての表面多角形を単一のモデルセルに割り当てるステップとを含む。さらに、システムは、１つ以上の細分面を使用して、割り当てられたすべての交差多角形を細分することを含む、表面多角形列によって取り囲まれた単一のモデルセルの部分を少なくとも細分するステップを含む。なおさらに、本方法は、細分ステップから得られる細分された個々のセルに表面多角形のすべてを再び割り当てるステップと、共通のシステム表面から、２つ以上の表面多角形を有するそれらの細分されたセルを識別するステップと、識別されたセルの共通の表面多角形に関し有効な共通の境界領域及び法線ベクトルを計算するステップとを含む。
【００１１】
他の形態において、本明細書に具体化されかつ広く記載されるように、１つ以上の表面を有する流体ダイナミクスシステムのマルチセルコンピュータモデル用の固定幾何学的格子を自動的に形成するための装置は、デジタルコンピュータと、コンピュータに記憶された計算流体ダイナミクスプログラムとを含む。ＣＦＤプログラムは、一列の隣接多角形によってシステム表面を表すデータセットを受け取るための、システム全体を取り囲む単一のモデルセルを創出するための、またすべての表面多角形を単一のモデルセルに割り当てるためのソフトウェアを含む。プログラムはまた、１つ以上の細分面を使用して、割り当てられたすべての交差多角形を細分することを含む、表面多角形列によって取り囲まれた単一のモデルセルの部分を少なくとも細分するためのソフトウェアを含む。さらに、プログラムは、細分ステップから得られる細分された個々のセルに表面多角形のすべてを再び割り当て、共通のシステム表面から、２つ以上の表面多角形を有するそれらの細分されたセルを識別し、またそれらの識別されたセルについて有効な境界領域及び法線ベクトルを計算するためのソフトウェアを含む。
【発明を実施するための最良の形態】
【００１２】
一形態では、１つ以上の表面を有する流体動的システムのマルチセルコンピュータモデル用の固定幾何学的格子を形成するための方法は、隣接する一列の多角形によってシステム表面を表すステップを含む。本明細書に具体化されるように、図１は、内燃機関（図示せず）用のヘッダパイプ１０を示しており、パイプ１０の表面１２は、パイプ１０のモデルを形成する際に使用するために「三角形に分けられる」。すなわち、流体（ガス、液体等）用の境界として使用される表面１２のそれぞれの部分は、一列の個々の三角形１４によって覆われる。しかし、本発明の方法は、システム表面を表すための三角形の類の多角形の使用に限定されない。
【００１３】
列１４の個々の三角形は、パイプ１０の比較的直線状の円筒面部分１２ａと１２ｂにおける細長い三角形１４ａと１４ｂのように、またパイプ１０の湾曲面部分１２ｃにおける三角形１４ｃのようなよりコンパクトな三角形のように、異なる寸法であることができる。本発明の開示した方法では、表面全体を閉じなければならず、それぞれの多角形は、図１の三角形の類の多角形のように、他の１つの多角形のみと辺を共有しなければならない。
【００１４】
これらの三角形のそれぞれは、図１Ａに拡大して図示した三角形１４ｂの頂点１６ａ、１６ｂ、１６ｃのようなそれぞれの三角形の３つの頂点のデカルト座標のような座標を規定するデータセットによって、３次元空間で一意に識別することができる。同様に、それぞれの表面が流体ダイナミクスシステムで異なって実行する多数の表面を有するシステムでは、識別子（例えば、数）を三角形のデータセットの部分として確立することが好ましいかもしれない。一般に、「一意な」表面のこのような識別子は、モデル化の段階で境界条件を適用するときに使用することができる。当業者は、パラメトリック・テクノロジー社（Ｐａｒａｍｅｔｒｉｃ  Ｔｅｃｈｎｏｌｏｇｙ  Ｃｏｒｐ．）（マサチューセッツ州ニーダム（Ｎｅｅｄｈａｍ，Ｍａｓｓａｃｈｕｓｅｔｔｓ））から入手可能なプロ／エンジニア（Ｐｒｏ／Ｅｎｇｉｎｅｅｒ）のような公知のＣＡＤプログラムを使用して、このような三角形に分けられた表面及び識別データセットを用意できるであろう。
【００１５】
さらに、この方法は、モデル化すべきフローシステム全体を取り囲む単一のモデルセルを創出するステップを含む。本明細書に具体化されるように、図２は、パイプ１０の内側境界を画定する代表的な表面三角形１４ｂを例えば含む三角形に分けられた表面列１４を完全に取り囲むように寸法決めされた単一のデカルトセル２０を示している。
【００１６】
さらに、記載した方法は、すべての表面多角形を単一のモデルセルに割り当てるステップを含む。本明細書に具体化されるように、また図３Ａを参照すると、列１４からの代表的な表面三角形が、モデル格子内に「マッピングされて」いるのが示されており、このモデル格子は、第１の例では、形成すべき格子の前駆体として図２の単一セル２０を含む。実際に、マッピングされた列１４は、格子内の閉じた流体システム、すなわち示した実施形態のパイプ１０の輪郭を描いている。したがって、オクツリーデータセットであり得る格子を表すデータセットは、以下に説明するようにモデル格子を提供するためにさらなる処理において使用するために、パイプ１０の表面に対応する列１４の個々の三角形の上の情報を含む。
【００１７】
なおさらに、開示した方法は、１つ以上の細分面を使用して、割り当てられたすべての交差多角形を細分することを含む表面多角形列を取り囲む単一のモデルセルの部分を少なくとも細分するステップを含む。本明細書に具体化されるように、図３Ａと図３Ｂは、三角形に分けられた表面列１４によって画定されたヘッダパイプモデル容積内に個々の部分を有する８つの細分されたセル２８、３０、３２、３４、３６、３８、４０、４２を創出するために、３つの相互直交面２２、２４、２６によって分割される単一セル２０を示している。
【００１８】
当業者に公知のように、細分面２２、２４、２６はまた、図３Ａの代表的な三角形１４ｂのような１つ以上の面によって交差される任意の表面三角形を分割し、三角形１４ｂは、表面多角形５０と５２を形成するために面２４によって交差されているのが示されている。これらの表面多角形は、多角形５２のような三角形、又は四辺形の表面多角形５０のような４つ以上の辺を有する幾何学的図形であり得る。
【００１９】
さらに、記載した方法は、すべての表面多角形、すなわち細分された多角形及び分割されていない任意の残りの多角形を、細分するステップから得られる細分されたセルの個々のセルに再び割り当てるステップを含む。本明細書に具体化されるように、図３Ａと図３Ｂの代表的な三角形１４ａのような分割されていない表面三角形、及び表面多角形５０と５２は、細分されたセルに「マッピング」される。具体的に、図３Ｂの表面三角形１４ａが、下方右の細分されたセル２８にマッピングされる。表面多角形５０はまた、細分されたセル２８にマッピングされ、一方、表面多角形５２は、上方右の細分されたセル３６にマッピングされる。例えば、マッピングは、オクツリーフォーマットのような適切なデータ記録フォーマットで、格子に対するそれぞれの表面三角形及び表面多角形の頂点の一意なデカルト（ｘ、ｙ、ｚ）位置を考慮することによって達成することができる。
【００２０】
細分するステップ及び再び割り当てるステップは、細分されたセルが所望の格子サイズに等しいか、又はそれ未満になるまで繰り返されることが好ましいかもしれない。例えば、図４Ａは、２次元で、図１のヘッダパイプ１０用のモデルについて得られる格子５４を示している。格子５４は、表面多角形の細分された、しかしなお閉じた隣接列１４’内に閉じ込められる。
【００２１】
特に、温度、圧力、速度等において予想される大きな勾配を有するシステム位置に、格子５４の平均の非境界セル５６に対しより微細な格子の寸法を有するセル（図示せず）を包埋するような格子５４の別の修正を行うことが可能である。当業者は、このような修正を行う方法を知っているであろう。
【００２２】
なおさらに、この方法は、２つ以上の表面多角形を有するそれらの細分されたセルを共通のシステム表面から識別するステップを含む。本明細書に具体化されるように、図４Ａ〜図４Ｃのセル５８は、１回以上細分された表面三角形１４ｃから生じた複数の多角形を有するこのような境界セルである。明快さ及び理解しやすくするため、２つのみの表面多角形６０と６２が図４Ｂに示されている。開示した実施形態において、２つの表面多角形は、共通のシステム表面、すなわち、曲がりの近傍の、しかし、パイプ１０を貫通するポート（番号なし）の表面１２ｄから離間したパイプ１０の表面１２ｃの部分から生じると考えることができる。
【００２３】
記載した本発明の方法を使用する当業者は、モデル化されるシステムに対する予想される流体力学的及び／又は熱力学的効果に基づき、このような共通の表面を識別できるであろう。例えば、パイプ１０内の鋭角隅部を形成する表面、又は相対的に運動する表面は、「共通の表面」であると予想されないであろう。同様に、パイプ１０の表面１２の外部加熱されるか又は冷却される局所的な部分は、特定の用途及び実施される試験に応じて、表面１２の隣接する未加熱／未冷却の部分と「共通」の表面であると考えてもよく、考えなくてもよい。
【００２４】
なおさらに、開示した方法は、識別されたセルの共通の表面多角形に関し有効な境界領域及び境界法線ベクトルを計算するステップを含む。当業者が理解するように、それぞれの別個の多角形を所定のセル内の別個の境界として取り扱うことは、モデルの使用中に、そのセルの輸送方程式を解くときにそれぞれの多角形に適切な境界条件を適用するために、別個の計算時間（またデジタルコンピュータによって実施される場合には別個のメモリ）を必要とするであろう。共通の表面多角形用の単一の有効表面積及び法線ベクトルをそれぞれのセルに設けることによって、例えば、次の手順によって、幾何学的形状を変更することなく、計算時間及びメモリの両方の節減を達成することができる。
【００２５】
対流輸送の有限体積解法では、境界を横切るフラックスが、次の方法で、辺領域フラックスと共に加算される。
【００２６】
【数１】
2007179548.tif 000003
【００２７】
同一の境界のすべての表面多角形について、φｂ及びｕｂは同一である。次に、残りの領域及び法線ベクトルに次式を結合することができる。
【００２８】
【数２】
2007179548.tif 000004
【００２９】
【数３】
2007179548.tif 000005
ベクトルはまた、単位ベクトル
【００３０】
【数４】
2007179548.tif 000006
として扱うことができ、大きさ、Ａｓｕｍは、次式によって定義されるような平均された有効境界の面積に等しい。
【００３１】
【数５】
2007179548.tif 000007
【００３２】
この結合方法は、単一セルに関する境界の幾何学的形状情報のすべてが、結合された境界の投影面積に等しい大きさを有する単一の法線ベクトルに変形されることを可能にする。これにより、メモリ及び計算時間も節約されるが、この理由は、セルの変数の変化が次式のように表すことができるからであり、
【００３３】
【数６】
2007179548.tif 000008
ここで、すべての「同様の」又は共通の境界は、単一の境界として取り扱われる。同様の関係が拡散フラックスに関し存在し、同一の有効境界領域及び法線ベクトルをもたらす。
【００３４】
本明細書に具体化されるように、計算するステップは、それぞれの識別されたセルについて、共通の表面多角形のそれぞれから法線ベクトルを計算することを含み、それぞれのベクトルの大きさは、個々の共通の表面多角形の面積である。図４Ｂで最善に理解されるように、ベクトル６４は、表面多角形６０に対し法線に（直角に）計算され、その多角形の面積に等しい大きさ（長さ）を有する。同様に、ベクトル６６は、表面多角形６２に法線に計算され、多角形６２の面積に等しい大きさを有する。ベクトル方法は著しく有用であるが、この理由は、一般に、共通の表面多角形が、例えばパイプ１０の湾曲面１２を表す多角形について、互いに又は格子セルの立方面と同一平面でなくてもよいからである。
【００３５】
さらに、計算するステップは、それぞれの識別されたセルについて、共通の表面多角形のそれぞれから、単一の法線ベクトルの方向の表面多角形の結合された投影面積に等しい大きさを有する単一の法線ベクトルに、すべての法線ベクトルを結合するステップを含む。本明細書に具体化されるように、ベクトル和
【００３６】
【数７】
2007179548.tif 000009
は上述のように計算される。ベクトル和は、図４Ｂの多角形６０と６２の幾何学的特性を表す。当業者は、
【００３７】
【数８】
2007179548.tif 000010
のｘ、ｙとｚ成分を使用して、多角形６０と６２に有効な境界領域を通してフラックスを配分でき、これらの成分は、セル５８の面７２（辺）、７４（頂部）及び／又は７６（後部）を通して隣接セル（図示せず）に又はそこから輸送された熱物理学的な値の量（例えば、質量、運動量、エネルギ）に加えるか又は差し引くことができることを認識するであろう。
【００３８】
開示した本発明の方法はまた、従来のいくつかの計算による流体動的モデル化プロセスにおける問題及び／又は過度の労力を生じるであろうモデルセルタイプを有利に処理することが可能である。具体的に、流体システムの（セル幅に対して）薄い固体表面が、表面によって分割される２つ以上の部分にセルを「分割する」場合、従来のカットセルＣＦＤ方法は、それらが接続されたかのように、分離された分割セル部分を取り扱う。このような従来の方法は、計算性能の非常に大きなエラーをもたらすことがある。例えば、表面の片側のセル部分が、このような従来のＣＦＤ方法で隣接セルからエネルギ（例えば、熱）フラックスを受け取った場合、当該のフラックスは、セルの輸送方程式を解くときに、結合されたセル部分が全体として扱われる場合、他のセル部分によって瞬時に「感じられる」であろう。
【００３９】
このような不調和な結果を回避するために、開示した本発明の方法は、分割セルのそれぞれのセル部分と個々の隣接境界セル及び非境界セルとを接続する情報を決定するステップを含んでもよい。本明細書に具体化されるように、図５は（説明及び視覚化を容易にするため２次元で）、薄い固体表面８４によって交差される細分セル８２からなるモデル格子８０の部分を示している。例えば、このような表面は、ダクトに配置された可動のスロットルプレートの部分であり得る（両方とも図示せず）。セルＳは、表面８４によって２つの部分Ｓ１とＳ２に分割される。セル部分Ｓ１は、隣接部として境界セル部分Ａ１とＣ１、及び非境界セルＢを有する。セル部分Ｓ２は、隣接部として境界セル部分Ａ２とＣ２、及び非境界セルＤを有する。本開示において、「隣接する」セルは、共通のインターフェイス領域を共有するセルであり、境界セルは、１つ以上の表面三角形又は他の多角形を有するセルである。
【００４０】
開示した本発明の方法で、フラックスをセル部分Ｓ１に算入しまたそれから算出するために、Ｓ１及びＡ２、Ｃ２及びＤの間の輸送方程式でなく、Ｓ１及びＡ１、Ｃ１及びＢの間の輸送方程式が解かれるであろう。同様に、セル部分Ｓ２について決定された接続情報は、部分Ａ１、Ｃ１、及びセルＢでなく、境界セル部分Ａ２、Ｃ２及び非境界セルＤで解かれる輸送方程式を提供するであろう。すなわち、分割セルのそれぞれの部分は無関係に、それ自体の熱物理学的な特性を有すると取り扱われる。
【００４１】
しかし、開示した本発明の方法によれば、両方のセル部分は同一の格子位置情報を保持するであろう。すなわち、図５の実施形態の格子８０の実装中に、Ｓ１とＳ２の両方は、予め分割されたセルＳと同一の格子位置を有するであろう。同様に、Ａ１とＡ２は、予め分割されたセルＡと同一の格子位置を有し、Ｃ１とＣ２は、予め分割されたセルＣと同一の格子位置を有するであろう。当業者は、複雑な薄い表面は、図５に示した２つの部分よりも多く、また３次元格子の部分よりも多くセルを分割でき、分割セル部分は、セル部分Ｓ１とＳ２について示されるよりも多くの隣接部を有することが可能であることを理解するであろう。
【００４２】
従来のＣＦＤプロセスに問題を生じることがある他のモデルセルタイプは、いわゆる「スリバ」セルである。スリバセルは、格子の典型的な非境界セルよりも著しく小さな容積を有する境界セルである。このような著しく小さなセルは、当業者が認識するように、モデルを使用して過渡の又は準定常の熱物理学的な計算中に安定性を保証するために、非常に小さな時間ステップ又は多数の反復を必要とすることがある。本方法のデジタルコンピュータによる実施における小さな時間ステップ又は多数の反復は、ＣＦＤ方法にとって非常に長い実行時間を意味する可能性がある。
【００４３】
本発明の開示によれば、固定幾何学的格子を形成する本方法は、非境界セルの容積の予め選択された部分のような予め選択された量未満の容積を有する境界セルを識別するステップを含む。本明細書に具体化されるように、図６は、ヘッダパイプ１０をモデル化するための完成した格子５４の図４Ａの２次元表示の詳細図を示している。細分されたセルＴ、Ｖ、Ｗ、Ｙ、Ｚは、細分された表面三角形列１４’によって境界付けられないような境界セルではなく、すべては同一の体積を有し、すなわち、最終の格子サイズセル５６に対応する。しかし、境界セルＲ、Ｕ、Ｘのそれぞれは、セル５６の容積部分を有し、セルＵは最小部分を有する。以下の説明のために、セルＵは、予め選択された最小の所望のセルの部分容積未満の容積を有するが、両方のセルＸ及びセルＲは、予め選択された部分容積よりも大きな部分容積を有することが想定される。
【００４４】
さらに本開示によれば、本方法は、識別されたそれぞれの部分境界セルと隣接セルとを対にするステップと、対にしたそれぞれの隣接セルの幾何学的情報を調整して、識別された個々の部分境界セルの容積及び境界領域情報を少なくとも含めるステップとを含む。本明細書に具体化されるように、図６の部分境界セルＵは非境界セルＶと対にされ、このセルＶと共通面９０を共有する。しかし、対にされたセルは境界セルであることができる。フラックスをセルＶに算入しまたそれから算出するために、セルＶの容積は、セルＵの容積を含むように調整されるであろう。同様に、Ｖの面情報は、境界セルＸと共通のセルＵの面９２の情報を含むように調整されるであろう。セルＵとセルＶとの間の共通面９０に関する情報は、本方法の実行中に輸送方程式がセルＵとセルＶとの間で解かれないので、消去してもよい。しかし、隣接する非境界セルＶ及び識別された部分境界セルＵに関する別個の格子位置情報は、保持される。
【００４５】
本方法の実行中、フラックスをセル対（Ｕ、Ｖ）の内の調整されたセルＶのみに又はそれから決定するために、輸送方程式が解かれ、次に、セルＶの計算された特定のシステム設計値（例えば、密度、温度、圧力等）がセルＵに帰せられる。当業者は、図４Ａと図６の実際の３次元格子に関し、セルＵのセル面９２に追加してセル面が隣接セルＶに加えられ、すなわち、セルＵは、その隣接セルと共通に、図４Ａと図６の２次元面内にまたそれから出る方向に面することを理解するであろう。また、２つよりも多くのセルを「対に」してもよい。
【産業上の利用可能性】
【００４６】
他の形態によれば、１つ以上の表面を有する流体ダイナミクスシステムのマルチセルコンピュータモデル用の固定幾何学的格子を自動的に形成するための装置は、デジタルコンピュータと、固定幾何学的格子を有するマルチセルモデルを実行するためにコンピュータにロードされた計算流体ダイナミクスプログラムとを有する。本明細書に具体化されるように、図７は、前述の方法を実施するように修正されたＭｏＳＥＳプログラムのような計算による流体ダイナミクスプログラム１１２でプログラミングされたＡＭＤ「オプテロン（Ｏｐｔｅｒｏｎ）」コンピュータのようなデジタルコンピュータ１１０を含む装置１００を概略的に示している。プログラム１１２は、図４Ａの格子５４によって表されるように、図１に示したインテークパイプ１０の３次元マルチセルモデルを実行するように構成することができる。
【００４７】
本明細書にさらに具体化されるように、プログラム１１２は、自動化された格子形成のための上に開示した方法を実施するためのソフトウェアを有するプログラムルーチン１２４を含むように修正される。ルーチン１２４のソフトウェアは別個のプリプロセッサプログラムに含めることができるであろうが、本明細書に開示した自動化された格子形成方法は、格子配向、解像度、又は表面の幾何学的形状を変更するための手動介入を必要とせず、ＭｏＳＥＳ  ＣＦＤプログラムと容易に一体化して、ＣＦＤプログラムの実行中に、特に固定幾何学的格子を有する運動面のモデル化中に、方法への容易なアクセスを可能にすることができる。このような場合、運動境界に当てはまるように、実行時中に格子を再形成する必要があり得るので、商業的に実行可能な方法にとって、格子形成の速度が重要であり得る。本明細書に開示した方法は、慣例的な方法よりも相当高速であることが示されている。
【００４８】
具体的に、ルーチン１２４は、図１に示した三角形列１４のような一列の隣接多角形によってヘッダパイプの固体表面を表すデータセットを受け取るためのソフトウェアを含む。さらに、ルーチン１２４は、図２のセル２０のようなシステム全体を取り囲む単一のモデルセルを創出するためのソフトウェアを含む。なおさらに、ルーチン１２４は、図２に示した代表的な三角形１４ａ、１４ｂ、１４ｃのようなすべての表面多角形を単一のモデルセルに割り当てるためのソフトウェアを含む。
【００４９】
なおさらに、ルーチン１２４は、図３Ａの面２２、２４、２６のような１つ以上の細分面を使用して割り当てられたすべての交差多角形を細分することを含む、表面多角形列によって取り囲まれた単一のモデルセルの部分を少なくとも細分するためのソフトウェアを含む。図３Ａはまた、図３Ａの代表的なセル１４ｃを交差する面から形成された表面多角形５０と５２を示し、図３Ｂは、前に説明したように、単一セル２０から細分された８つのセル２８、３０、３２、３４、３６、３８、４０を示している。
【００５０】
再び、三角形列１４によって表されたシステム内の細分されたセルの部分が対象とされ、細分されたセルの立方構造全体は、見やすくするために図３Ａと図３Ｂに示されている。
【００５１】
なおさらに、ルーチン１２４は、もしあれば、残りの表面三角形のようなすべての表面多角形、及び図３Ａと図３Ｂに示した細分プロセスによって創出された表面多角形を細分された個々のセルに再び割り当てるためのソフトウェアを含む。ソフトウェアはまた、図４Ａの格子５４のように、所望の格子サイズが達成されるまで細分及び再割当の繰り返しを行うことが可能である。
【００５２】
なおさらに、ルーチン１２４は、共通のシステム表面から、２つ以上の表面多角形を有するそれらの細分されたセルを識別して、このようなセルの有効な共通境界領域及び法線ベクトルを計算するためのソフトウェアを含む。前に説明したように、図４Ｂは、図１のヘッダパイプ１０の共通表面、すなわち表面三角形１４ｄによって表される表面部分１２ｄから遠い表面１２ｃの部分から、表面多角形６０と６２を有する識別されたセル５８を示している。特定の表面に適切な境界条件（例えば、分割された多角形６０と６２を通したセル５８からの熱フラックスを規定するための熱伝達係数）が設けられた場合、前に説明したように計算される有効境界領域及び法線ベクトルは、表面多角形６０と６２を別個に扱うであろう従来のＣＦＤコードに対し、セル５８の輸送方程式を解くときに計算時間を低減する。当業者は、本発明の開示を前提として、共通の有効境界領域及び法線ベクトルを提供するように、ソフトウェアを創出するか又は適合させることができるであろう。
【００５３】
さらに、ルーチン１２４は、図５と図６に関し前に説明した方法をそれぞれ使用して、「スリバ」型境界セル及び／又は「分割」型境界セルを取り扱うためのソフトウェアを含んでもよい。再び、当業者は、本発明の開示を前提として、ルーチン１２４に適切なソフトウェアを提供できるであろう。
【００５４】
本発明の方法及び装置は、本出願と同時に出願された（特許文献１）に開示された埋め込みセルを有するマルチセルコンピュータモデル用のマルチグリッド計算を実施するための方法及び装置（Ｍｅｔｈｏｄ  ａｎｄ  Ａｐｐａｒａｔｕｓ  ｆｏｒ  Ｉｍｐｌｅｍｅｎｔｉｎｇ  Ｍｕｌｔｉ−Ｇｒｉｄ  Ｃｏｍｐｕｔａｔｉｏｎ  ｆｏｒ  Ｍｕｌｔｉ−Ｃｅｌｌ  Ｃｏｍｐｕｔｅｒ  Ｍｏｄｅｌｓ  ｗｉｔｈ  Ｅｍｂｅｄｄｅｄ  Ｃｅｌｌｓ）と関連して使用することが好ましいかもしれない。
【００５５】
同様に、本発明の方法及び装置は、本出願と同時に出願された（特許文献２）に開示された流体動的システムのマルチセルコンピュータモデルの運動境界を取り扱うための方法及び装置（Ｍｅｔｈｏｄ  ａｎｄ  Ａｐｐａｒａｔｕｓ  ｆｏｒ  Ｔｒｅａｔｉｎｇ  Ｍｏｖｉｎｇ  Ｂｏｕｎｄａｒｉｅｓ  ｉｎ  Ｍｕｌｔｉ−Ｃｅｌｌ  Ｃｏｍｐｕｔｅｒ  Ｍｏｄｅｌｓ  ｏｆ  Ｆｌｕｉｄ  Ｄｙｎａｍｉｃ  Ｓｙｓｔｅｍｓ）と関連して使用することが好ましいかもしれない。
【００５６】
なお同様に、本発明の方法及び装置は、本出願と同時に出願された（特許文献３）に開示された動的システムのマルチセルコンピュータモデルの輸送方程式を解くための方法及び装置（Ｍｅｔｈｏｄ  ａｎｄ  Ａｐｐａｒａｔｕｓ  ｆｏｒ  Ｓｏｌｖｉｎｇ  Ｔｒａｎｓｐｏｒｔ  Ｅｑｕａｔｉｏｎｓ  ｉｎ  Ｍｕｌｔｉ−Ｃｅｌｌ  Ｃｏｍｐｕｔｅｒ  Ｍｏｄｅｌｓ  ｏｆ  Ｄｙｎａｍｉｃ  Ｓｙｓｔｅｍｓ）と関連して使用することが好ましいかもしれない。
【００５７】
本発明の他の実施形態は、本明細書に開示した本発明の説明と実施を考慮すれば当業者には明白であろう。説明および実施例は模範的なものに過ぎないと考えられ、本発明の真の範囲及び精神は、次の特許請求の範囲によって示されることが意図される。
【図面の簡単な説明】
【００５８】
【図１】自動化された格子形成用の開示した模範的な方法に使用するための、内燃機関のヘッダパイプの三角形に分けられた内面の斜視図である。
【図１Ａ】図１の一部の詳細図である。
【図２】自動化された格子形成用の開示した模範的な方法に従って形成された格子の単一セル前駆体の斜視図である。
【図３Ａ】自動化された格子形成用の開示した模範的な方法による、３つの直交面によって分割される図２の単一セル前駆体の斜視図である。
【図３Ｂ】自動化された格子形成用の開示した模範的な方法による、図３Ａの単一セルの分割面の動作から得られる細分されたセルの斜視図である。
【図４Ａ】自動化された格子形成用の開示した模範的な方法による、図１に示したような内面を有するヘッダパイプをモデル化するために形成された格子の２次元の切り取り平面図である。
【図４Ｂ】自動化された格子形成用の開示した模範的な方法による、表面に多角形が割り当てられた図４Ａの形成された格子の境界セルの詳細斜視図である。
【図５】自動化された格子形成用の開示した模範的な方法による、プレート面により分割境界線セルが形成される可動のスロットルプレートを有するダクトをモデル化するために形成された格子の詳細平面図である。
【図６】自動化された格子形成用の開示した模範的な方法による、図４Ａの格子のスリバ境界セルの切り取り平面図である。
【図７】自動化された格子形成用の開示した模範的な方法を実施するための開示した模範的な装置の斜視図である。
【符号の説明】
【００５９】
１０  ヘッダパイプ
１２  パイプ表面
１２ａ  パイプ表面
１２ｂ  パイプ表面
１２ｃ  パイプ表面
１２ｄ  ポート表面
１４  三角形列
１４’  細分された三角形列
１４ａ  表面１２ａの個々の三角形
１４ｂ  表面１２ｂの個々の三角形
１４ｃ  表面１２ｃの個々の三角形
１６ａ  代表的な三角形１４ｂの頂点
１６ｂ  代表的な三角形１４ｂの頂点
１６ｃ  代表的な三角形１４ｂの頂点
２０  単一のモデルセル
２２  分割面
２４  分割面
２６  分割面
２８  細分された８つのセルの１つ
３０  細分された８つのセルの１つ
３２  細分された８つのセルの１つ
３４  細分された８つのセルの１つ
３６  細分された８つのセルの１つ
３８  細分された８つのセルの１つ
４０  細分された８つのセルの１つ
４２  細分された８つのセルの１つ
５０  表面多角形
５２  表面多角形
５４  形成された格子
５６  非境界セル
５８  境界セル
６０  表面多角形
６２  表面多角形
６４  表面多角形の法線ベクトル
６６  表面多角形の法線ベクトル
７２  セル５８の面
７４  セル５８の面
７６  セル５８の面
８０  格子部分
８２  格子８０の細分されたセル
８４  薄い表面
９０  共通のセル面
９２  セル面
Ａ  格子８０の境界セル
Ａ１  境界セルＡの部分
Ａ２  境界セルＡの部分
Ｂ  格子８０の非境界セル
Ｃ  格子８０の境界セル
Ｃ１  境界セルＣの部分
Ｃ２  境界セルＣの部分
Ｄ  格子８０の非境界セル
Ｒ  格子５４の境界セル
Ｓ  格子８０の境界セル
Ｓ１  セルＳの部分
Ｓ２  セルＳの部分
Ｔ  格子５４の非境界セル
Ｕ  格子５４の境界セル
Ｖ  格子５４の非境界セル
Ｗ  格子５４の非境界セル
Ｘ  格子５４の境界セル
Ｙ  格子５４の非境界セル
Ｚ  格子５４の非境界セル
【図１】
2007179548.tif 000011
【図１Ａ】
2007179548.tif 000012
【図２】
2007179548.tif 000013
【図３Ａ】
2007179548.tif 000014
【図３Ｂ】
2007179548.tif 000015
【図４Ａ】
2007179548.tif 000016
【図４Ｂ】
2007179548.tif 000017
【図５】
2007179548.tif 000018
【図６】
2007179548.tif 000019
【図７】
2007179548.tif 000020