JP2014174990A - 流体領域に対して移動メッシュを用いる空間−時間ｃｅ／ｓｅソルバーを利用したｆｓｉの数値的シミュレーション - Google Patents

流体領域に対して移動メッシュを用いる空間−時間ｃｅ／ｓｅソルバーを利用したｆｓｉの数値的シミュレーション Download PDF

Info

Publication number
JP2014174990A
JP2014174990A JP2014036188A JP2014036188A JP2014174990A JP 2014174990 A JP2014174990 A JP 2014174990A JP 2014036188 A JP2014036188 A JP 2014036188A JP 2014036188 A JP2014036188 A JP 2014036188A JP 2014174990 A JP2014174990 A JP 2014174990A
Authority
JP
Japan
Prior art keywords
fluid
mesh
time
space
fsi
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
JP2014036188A
Other languages
English (en)
Inventor
Zhang Zeng-Chan
ジャン ゼン−チャン
cook Grant
クック グラント
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Livermore Software Technology LLC
Original Assignee
Livermore Software Technology LLC
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Livermore Software Technology LLC filed Critical Livermore Software Technology LLC
Publication of JP2014174990A publication Critical patent/JP2014174990A/ja
Pending legal-status Critical Current

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00—Details relating to CAD techniques
    • G06F2111/10—Numerical modelling

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • Algebra (AREA)
  • Computing Systems (AREA)
  • Fluid Mechanics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Physics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

【課題】構造力学を数値的にシミュレートする方法と連係する流体メッシュを用いる空間−時間ＣＥ／ＳＥ法を用いて、ＦＳＩを数値的にシミュレートするシステムおよび方法を提供する。
【解決手段】ＦＳＩ境界面が、コンピュータシステムにおいて受け取った流体領域および構造体定義に基づいて決定される。ＦＳＩ境界面に作用する流体力が初期化される。シミュレートされた構造体挙動が、ＦＳＩ境界面における受け取った流体力に応じて、ＦＥＡを用いて取得される。構造体挙動には、空間−時間流体メッシュのＦＳＩ境界の更新に用いられる構造体の外部の境界上のノード位置が含まれる。流体メッシュの内部ノードは、ユーザに選択されたメッシュ調整法を用いて調整される。シミュレートされた流体挙動は、調整された流体メッシュで、ＣＥ／ＳＥソルバーを用いて流体ソリューションを更新することにより、取得される。流体力がＦＥＡモデルに再び印加され、次のソリューションサイクルに対するシミュレートされた構造体挙動を取得する。
【選択図】図５Ａ

Description

本発明は、コンピュータ支援エンジニアリング（ＣＡＥ）解析に関し、特に構造力学（つまり構造体挙動）を数値的にシミュレートする方法と連係する移動空間−時間流体メッシュを用いる空間−時間保存要素・ソリューション要素（conservation element and solution element（ＣＥ／ＳＥ））法を利用したＦＳＩ（流体−構造体相互作用）の数値的シミュレーションに関する。
コンピュータ支援エンジニアリング解析は、対象の数値的にシミュレートされた応答／結果（例えば構造体挙動、流体運動など）を取得するよう構成される。そして、シミュレートされた応答／結果は、製品（例えば自動車、飛行機など）を改良するよう設計判断するために、あるいは視覚化するのが困難または不可能な特定の物理的現象を調査するために、技術者および／または科学者によって用いられる。
コンピューティング技術の進歩によって、個々の数値的シミュレーションにおいて構造体挙動または流体運動のいずれかを取得する代わりに、流体および構造体モデリングの組み合わせシステムが、流体−構造体相互作用（ＦＳＩ）（例えば飛行中の飛行機、海上の船など）の数値的シミュレーションにおいて用いられている。
数値的にＦＳＩをシミュレートするための従来技術アプローチは、空間と時間とを別々に処理する方法で行われていた。しかしながら、高速流体の場合、流体シミュレーションでの不正確さが問題となる。空間−時間ＣＥ／ＳＥ（保存要素／ソリューション要素）法と呼ばれる別のアプローチが、流体シミュレーションに用いられる。しかし、空間−時間ＣＥ／ＳＥ法における従来技術アプローチは、オイラーグリッド／メッシュすなわち固定グリッド／メッシュ（すなわち、メッシュが数値的シミュレーション全体で一定のままである）に基づいて、空間−時間領域において流体（すなわち、空気）を、オイラーグリッドを移動している他のグリッドモデル（例えば有限要素解析モデル）によって表される構造体（すなわち、航空機）とともに、表している。しかしながら、固定のオイラーグリッドにおけるＦＳＩインタフェースの結果、ある程度の精度が失われる。
したがって、構造力学を数値的にシミュレートする方法と連係する移動空間−時間流体メッシュを用いる空間−時間ＣＥ／ＳＥ法を用いて、ＦＳＩを数値的にシミュレートする改良技術が望まれよう。
構造力学を数値的にシミュレートする方法と連係する流体メッシュを用いる空間−時間ＣＥ／ＳＥ法を用いて、流体−構造体相互作用（ＦＳＩ）を数値的にシミュレートするシステムおよび方法が提供する。本発明の一の例示的な実施形態では、流体領域定義と、構造体定義（例えば飛行機、自動車など）と、がコンピュータシステムに受け取られる。流体領域は空間−時間流体メッシュによって表される。また、構造体は有限要素解析（ＦＥＡ）モデルによって表される。流体領域にはさらに、流体変数（例えば、密度、速度、圧力、粘度など）が含まれる。ＦＳＩ境界面は、受け取った定義から決定される。次に、ソルバーの状態変数が初期化される。その後、ＦＳＩ境界面に作用する流体力が、ＦＳＩの時間進行数値的シミュレーションの開始において初期化される。
構造体の数値的にシミュレートされた構造体挙動は、ＦＳＩ境界面における受け取った流体力に応じてＦＥＡモデルを用いるＦＥＡによって取得される。限定するものではないが、構造体挙動には、構造体の外部の境界上のノード位置が含まれる。構造体の外部境界上のノード位置は、空間−時間ＣＥ／ＳＥ流体メッシュのＦＳＩ境界面境界を更新するのに用いることができる。それに応じて、流体メッシュの内部ノードは、ユーザに選択されたメッシュ調整法を用いて、境界条件としてＦＳＩ境界面における運動を使って、調整される。数値的にシミュレートされた流体挙動（例えばＦＳＩ境界面の流体力）が、調整された空間−時間流体メッシュで、ＣＥ／ＳＥソルバーを用いて流体ソリューションを更新することにより、取得される。流体力がＦＥＡモデルに再び印加され、ソリューション時間が進められた次のソリューションサイクルに対するシミュレートされた構造体挙動を取得する。予め定義された終了条件に達するまで、ＦＳＩの数値的シミュレーションが継続される。
本発明の目的、特徴および効果は、添付した図面を参照し、以下の本発明の実施の形態の詳細な説明を考察することによって明らかとなろう。
本発明のこれらおよび他の特徴、面および利点は、以下の説明、添付の特許請求の範囲および添付した図面を考慮してより理解されよう。図面は次の通りである。
種々の例示的な流体領域定義および構造体定義を示す図である。 種々の例示的な流体領域定義および構造体定義を示す図である。 種々の例示的な流体領域定義および構造体定義を示す図である。 例示的な構造体のＦＥＡモデルを示す図である。 本発明の一の実施形態にかかる、一つの空間的次元に対する空間−時間ＣＥ／ＳＥソルバーの例示的な機構を示す図である。 本発明の一の実施形態にかかる、一つの空間的次元に対する空間−時間ＣＥ／ＳＥソルバーの例示的な機構を示す図である。 本発明の一の実施形態にかかる、一つの空間的次元に対する空間−時間ＣＥ／ＳＥソルバーの例示的な機構を示す図である。 本発明の一の実施形態にかかる、一つの空間的次元に対する空間−時間ＣＥ／ＳＥソルバーの例示的な機構を示す図である。 本発明の一の実施形態にかかる、一つの空間的次元に対する空間−時間ＣＥ／ＳＥソルバーの例示的な機構を示す図である。 本発明の一の実施形態にかかる、二つの空間的次元に対する空間−時間ＣＥ／ＳＥソルバーの例示的な機構を示す図である。 本発明の一の実施形態にかかる、二つの空間的次元に対する空間−時間ＣＥ／ＳＥソルバーの例示的な機構を示す図である。 本発明の実施形態にかかる、空間−時間ＣＥ／ＳＥ法において用いることができる、固定オイラーメッシュと例示的な移動流体メッシュとの間の比較を示す図である。 本発明の実施形態にかかる、移動流体メッシュで、空間−時間ＣＥ／ＳＥソルバーを用いて、流体−構造体相互作用を数値的にシミュレートする例示的なプロセスを図示するフローチャートである。 本発明の実施形態にかかる、移動流体メッシュで、空間−時間ＣＥ／ＳＥソルバーを用いて、流体−構造体相互作用を数値的にシミュレートする例示的なプロセスを図示するフローチャートである。 本発明の一の実施形態にかかる、空間−時間流体メッシュ調整の例示的な一連のシーケンスを示す概略図である。 本発明の一の実施形態にかかる、空間−時間流体メッシュ調整の例示的な一連のシーケンスを示す概略図である。 本発明の一の実施形態にかかる、空間−時間流体メッシュ調整の例示的な一連のシーケンスを示す概略図である。 本発明の実施形態を実現可能である例示的なコンピュータの主要な部品を示すブロック図である。
構造力学を数値的にシミュレートする方法と連係する流体メッシュを用いる空間−時間ＣＥ／ＳＥ法を用いて、流体−構造体相互作用（ＦＳＩ）を数値的にシミュレートするシステムおよび方法を開示する。ユーザによって指定された構造体と流体領域との間のＦＳＩの時間進行シミュレーションを実行する。図１Ａ〜図１Ｃは、種々の例示的な流体領域１２０ａ〜１２０ｃと構造体１１０ａ〜１１０ｃとを示す図である。構造体１１０ａは全体的が流体領域１２０ａ内に位置することができ、あるいは、構造体１１０ｂは部分的に流体領域１２０ｂ内に位置することもでき、またあるいは、構造体１１０ｃは流体領域１２０ｃの直ぐ隣りに位置することもできる。
図１Ａに示す通り、ＦＳＩ境界面は、構造１１０ａの外部表面１３０ａ全体（楕円状の点線で示す）である。図１Ｂに示すＦＳＩ境界面は、流体領域１２０ｂと重なる、構造体１１０ｂの外部表面１３０ｂの一部分（円弧状の点線で示す）である。図１Ｃにおいては、ＦＳＩ境界面は、構造体１１０ｃと流体領域１２０ｃとの間の接触位置１３０ｃ（楕円状ドットで示す）である。
図１Ｄを参照して、構造体（例えば飛行機、自動車など）を表す例示的な有限要素解析（ＦＥＡ）モデル１００を説明する。荷重条件下の構造体挙動は、コンピュータシステム（例えば図７のコンピュータ７００）においてＦＥＡモデル用いる有限要素解析を用いて、数値的にシミュレートすることができる。本発明の一の実施形態において、構造体挙動は、構造体に隣接するまたは構造体を囲んでいる流体領域に対する空間−時間流体メッシュで、空間−時間保存要素／ソリューション要素（ＣＥ／ＳＥ）ソルバーを用いて取得される流体−構造体相互作用（ＦＳＩ）境界における流体荷重や力に応じて、ＦＥＡモデルを用いて数値的にシミュレートされる。他の物理学的機構では、構造体が移動するかつ／または形状が変化する（例えば熱膨張、化学反応など）場合もある。
空間−時間ＣＥ／ＳＥソルバーはＦＳＩ相互作用（例えば構造体変形）に応じて空間−時間流体メッシュによって表される流体領域の流体挙動を取得するように構成される。図２Ａ〜図２Ｅは、一つの空間的次元に対する空間−時間ＣＥ／ＳＥ法を例示する概略図である。図２Ａに、ＣＥ／ＳＥソルバーにおいて用いられるよう構成されたメッシュ２００を示す。メッシュ２００は、時刻軸（ｔ）２０１および空間軸（ｘ）２０２の二つの軸を用いて、流体領域の空間−時間領域（メッシュによって覆われた四角形領域）を表す。ＣＥ／ＳＥ法は、以下の偏微分方程式（ＰＤＥ）を考慮して記述できる。
Figure 2014174990
 ここで、ａは定数であり、uは流体領域の保存量(例えば密度、運動量、エネルギーなど)である。
各メッシュ点(j,n)２０４（黒丸で示す）は、ソリューション要素ＳＥ(j,n)２１４の中心に配置される。指数ｊおよびｎはそれぞれ空間軸２０２および時刻軸２０１に対する指数である。ＳＥ(j,n)２１４は、図２Ｂに示す点線の曲線によって境界付けられた空間−時間領域の内部であると定義される。ＳＥ(j,n)には、水平線セグメント、垂直線セグメント、およびそれらの直近部分が含まれる。
［００２３］
ＳＥ(j,n)内の任意の点(x,t)に対して、u(x,t)(保存量)およびh(x,t)(流束ベクトル)はそれぞれ以下の式により近似される。
Figure 2014174990
Figure 2014174990
 ここで、aは定数であり、un j、(ux)n j、(ut)n jはＳＥ(j,n)における定数であり、(xj、tn)はメッシュ点(j，n)２０４の座標である。数式３は定義h = (au,u)に類似した数値的な量である。
u = u*(x,t;j,n)はＳＥ(j,n)内において数式１を満たすものである。したがって、(ut)n j = -a(ux)n jを得る。これにより、数式２を以下の式に
Figure 2014174990
 ＳＥ(j,n)内の各点(x,t)に対して簡単に書ける。つまり、un jおよび(ux)n jだけがメッシュ点(j,n)２０４と関連する独立した進行変数である。
流体領域を、保存要素（ＣＥ）と呼ばれる、互いに重なり合わない四角形領域（図２Ａを参照）に分割する。図２Ｃおよび図２Ｄにそれぞれ示す通り、ＣＥ１(j,n)２２１およびＣＥ２(j,n)２２２が内部メッシュ点(j,n)２０４と関連付けられる。これらの二つのＣＥは、基本保存要素（ＢＣＥ）と呼ばれる。一方、図２Ｅに示すＣＥ(j,n)２２４は、合成保存要素（ＣＣＥ）と呼ばれ、ＣＥ１(j,n)２２１とＣＥ２(j,n)２２２との結合状態である。
ＣＥ1(j,n)２２１の境界を形成する線セグメントのなかで、ABおよびADはＳＥ(j，n)２１４に属し、CBおよびCDはＳＥ(j-1/2，n-1/2)に属している。同様に、ＣＥ２(j,n)２２２の境界は、ＳＥ(j，n)２１４とＳＥ(j+1/2，n-1/2)とに属する。その結果、各メッシュ点(j,n)２０４においては二つの保存条件を課する、つまり、
Figure 2014174990
 である。数式３および数式４を用いると、（ｉ）として、
Figure 2014174990
が得られる。1 - ν2 ≠ 0とすると、（ｉｉ）として、
Figure 2014174990
 となる。ここで、ν ≡ aΔt/Δxであり、(u+ x)n j ≡ (Δx/4(ux)n j)である。ソリューション法は数式６および数式７によって形成される。
数式５から、各ＢＣＥの境界から流出するh*の全流束はゼロである。二つの隣接するＢＣＥを分ける境界にわたる面積分は一のＳＥからの情報を用いて評価されるので、明らかに、局所的な保存関係（数式５）から全流束保存関係が得られる。つまり、ＢＣＥの任意の組み合わせの結合状態である任意の各空間−時間領域の境界から流出するh*の全流束もまたゼロになる。具体的には、ＣＥ(j,n)２２４は、ＣＥ１(j,n)２２１とＣＥ２(j,n)２２２との結合状態であり、
Figure 2014174990
 が、数式５から得られる。実際、数式８が数式６と等価であることを示すことができる。
ＣＥ／ＳＥ法の上記の説明は、簡単なＰＤＥに基づくものである。しかしながら、一つ、二つまたは三つの空間的次元における保存則のシステムを含むことができる一般的なＣＥ／ＳＥ開発の本質を表している。
図３Ａ〜図３Ｂは、二つの空間的次元に対する例示的な空間−時間流体メッシュを示す概略図である。図３Ａに示す通り、ｘ−ｙ面は、互いに重なり合わない凸四角形へと分割され、隣接する任意の二つの四角形は共通の辺を共有する。さらに、（ｉ）四角形の頂点（バーテックス）および重心は、黒丸および白丸でそれぞれ印付けしている。（ｉｉ）Qは典型的な四角形B1B2B3B4の重心である。（ｉｉｉ）A1、A2、A3およびA4はそれぞれ、四辺形B1B2B3B4に隣接する四角形の重心である。(iv)Q*(エックス"X"で印付ける)は多角形A1B1A2B2A3B3A4B4の重心である。点Q*は、一般に点Qと一致せず、重心Qに関連付けられたソリューション点と呼ばれる。点A1 *、A2 *、A3 *およびA4 *("x"で印付ける)は、点A1、A2、A3およびA4に対するそれぞれのソリューション点である。
任意のn > 0について、ｎ番目の時間レベル(t = nΔt (n = 0, 1/2, 1, 3/2, ・・・)における例示的なＣＥ/ＳＥメッシュを図3Bに示す。点Q、Q'およびQ"はそれぞれn番目の時間レベル、(n - 1/2)番目の時間レベルおよび(n + 1/2)番目の時間レベルの点を表示する(図3Aを参照)。Q'およびQ"は、点Qの共通の空間的投影である。図3Bに図示したものや、図示されていないもの（例示するまでもなく明瞭であるために）などの、他の空間−時間メッシュ点も同様に定義される。具体的には、点Q*、A1 *、A2 *、A3 *およびA4 *は、定義から、n番目の時間レベルにあり、それぞれが、点Q、A1、A2、A3およびA4に関連付けられた空間−時間ソリューション点であり、点Q'*、A1'*、A2'*、A3'*およびA4'*は、定義から、(n - 1/2)番目の時間レベルにあり、それぞれが、点Q'、A1'、A2'、A3'およびA4'に関連付けられた空間−時間ソリューション点である。
上記の定義を用いると、ＳＥ(Q*)で示す点Q*のソリューション要素は、五つの面Q'Q"B1"B1'、Q'Q"B2"B2'、Q'Q"B3"B3、Q'Q"B4"B4'ならびにA1B1A2B2A3B3A4B4、およびそれらの直ぐ隣接する部分の結合状態として、定義される。さらに、ＣＥk(Q)で示す（k = 1, 2, 3, 4)、点Qの四つの基本保存要素(ＢＣＥ)はそれぞれ、空間−時間柱A1B1QB4A1'B1'Q'B4'、A2B2QB1A2'B2'Q'B1'、A3B3QB2A3'B3'Q'B2'、およびA4B4QB3A4'B4'Q'B3'となるよう定義される。さらに、ＣＥ(Q)で示す、点Qの合成保存要素は、空間−時間柱A1B1A2B2A3B3A4B4A1'B1'A2'B2'A3'B3'A4'B4'、つまり、上記の四つのＢＣＥの結合状態となるよう定義される。
固定オイラーメッシュ４１０と例示的な移動空間−時間メッシュ４００とを比較する図を、図４に示す。ｘ−ｙ面におけるABCD内の各ソリューション要素ＳＥ(Q)に対して、Δtだけ隔たっている二つの時間のレベルtn-1とtnとの間には、保存要素ＣＥ(Q)（つまり、空間−時間多面体ABCDA'B'C'D'）がある。Δtは、時間進行シミュレーションの二つのソリューションサイクル間の時間インクリメント（増分）である。ＣＥ(Q)は、四つのＢＣＥとしてA'S'Q'R'ASQR、B'P'Q'S'BPQS、C'W'Q'P'CWQPおよびD'R'Q'W'DRQWを含む。一の実施形態において、中点規則がＣＥ/ＳＥ法における各ＣＥ(Q)に関する数式８の積分計算において用いられる。この場合、各ＣＥ(Q)に関する、上面、底面および側面(例えば側面A'S'AS)を含むすべての面の面積および外向きの単位法線ベクトルが必要となる。
固定オイラーメッシュ４１０において、すべての幾何学的なデータは一度だけ（初期化の際に）計算に必要とされる。移動メッシュ４００においては、幾何学的データは、時間進行シミュレーションの際に一定ではなく、したがって、すべてのソリューションサイクルにおいて更新されたデータが必要となる。さらに、移動メッシュ４００では、側面はすべて、二次元における空間−時間面や、三次元における空間−時間多面体と見なされる。さらに、固定オイラーメッシュ４１０では、法線ベクトルの時間方向成分はゼロである。移動メッシュ４００では、法線ベクトルの時間方向成分はゼロではないので、数式８の値を求める際に追加項が一つ追加される。
次に、図５Ａおよび図５Ｂに、構造力学を数値的にシミュレートする方法と連係する移動流体メッシュを用いる空間−時間保存要素／ソリューション要素（ＣＥ／ＳＥ）ソルバーを用いて、流体−構造体相互作用（ＦＳＩ）を数値的にシミュレートする例示的なプロセス５００を合わせて例示するフローチャートを示す。プロセス５００は、好ましくはソフトウェアで実行される。
プロセス５００は、ステップ５０２において、流体領域定義および構造体定義を、関連するアプリケーションモジュール（例えばＦＥＡソフトウェア、空間−時間ＣＥ／ＳＥソルバー・ソフトウェアなど）がインストールされたコンピュータシステム（例えば図７のコンピュータ７００）に受け取ることによって、スタートする。流体領域は、ＣＥ／ＳＥ法に基づくソルバーに対して構成された空間−時間流体メッシュによって表される。構造体は、有限要素解析（ＦＥＡ）モデル（例えば図１ＤのＦＥＡモデル１００）によって表される。限定するものではないが、流体領域にはさらに、流体密度、圧力、速度、粘度などが含まれる。ユーザは、空間−時間流体メッシュおよびＦＥＡモデルを流体領域定義および構造体定義として定義することができる。例えば、ユーザは、構造体または流体領域を表す体積（ボリューム）メッシュを指定することができる。あるいは、体積メッシュを、受け取った定義に基づいて、コンピュータシステムにインストールされたアプリケーションモジュールによって生成することができる。例えば、ユーザは、流体領域または構造体の外部表面を定義でき、定義された外部表面がコンピュータシステムにおいて受け取られる。その後、対応するＣＥ／ＳＥ流体メッシュまたはＦＥＡモデル（体積モデル）が、受け取った面定義に基づいて生成される。
ステップ５０４において、受け取った定義を用いて、流体領域と構造体との間の流体−構造体相互作用（ＦＳＩ）境界が決定される。本発明の一の実施形態では、共通のノードまたはエッジ、すなわち位置合わせされたノードまたはエッジを、空間−時間流体メッシュとＦＥＡモデルの間には必要としない。必要なのは、ＦＳＩ境界面を有する流体領域および構造体が、ほぼ同じ面（例えば図１Ａ〜図１Ｃに示すＦＳＩ境界面１３０ａ〜１３０ｃ）に位置することだけである。言い換えれば、ＦＳＩ境界面は構造体の外部境界の一部または全部と一致する。次に、ステップ５０６において、ソルバーの状態変数をすべて初期化した後、ＦＳＩの時間進行シミュレーションのパラメータが、例えばＦＥＡモデルのＦＳＩ境界面に作用する初期の流体力が、初期化される、
ステップ５０８において、シミュレートされた構造体挙動が、ＦＳＩ境界面における受け取った流体力に応じてＦＥＡモデルを用いるＦＥＡを実行するによって、取得される。限定するものではないが、シミュレートされた構造体挙動には、構造体の外部境界（例えばＦＳＩ境界面）のノード位置が含まれる。ＦＥＡは、明示的または暗黙的有限要素解析とできる。一の例示的なＦＥＡソフトウェアパッケージには、リヴァーモア・ソフトウェア・テクノロジー社によって提供されている製品ＬＳ−ＤＹＮＡＲ（登録商標）がある。
このようにして、プロセス５００では、ステップ５１０において、空間−時間流体メッシュが、ＦＥＡからＦＳＩ境界面において新しく取得されたノード位置（つまり構造体挙動）を用いて、更新される。次に、ステップ５１２において、流体メッシュの内部メッシュノードは、ユーザに選択されたメッシュ調整法を用いて（限定するものではないが、ボール−バーテックス法、逆距離加重法（インバース・ディスタンス・ウェイティング法）、放射基底関数法など）、ＦＳＩ境界面における更新されたノード位置に応じて調整される。
新しい流体メッシュが更新された後、ステップ５１４において、シミュレートされた流体挙動は、新しく調整された流体メッシュにおいてＣＥ／ＳＥソルバーを用いて流体ソリューションを実行することにより、取得される。シミュレートされた流体挙動には、ＦＳＩ境界面に作用する流体力が含まれる。具体的には、図５Ｂに示す通り、プロセス５００は、ステップ５１４ａにおいて、直前のソリューション時間および現在のソリューション時間の両方における空間−時間流体メッシュに基づいて、流体領域の幾何学的パラメータの計算を実行する。そして、ステップ５１４ｂにおいて、流体領域変数（例えば流体密度、圧力、速度、粘度など）および対応する空間的導関数が計算される。直前のソリューション時間と現在のソリューション時間とは、時間増分Δtだけ隔たっている。
次に、ステップ５１６において、時間進行シミュレーションの現在のソリューション時間が次のソリューションサイクルにインクリメントされる（例えば、時間増分Δｔを加えることにより現在のソリューション時間をインクリメントする）。プロセス５００は、ＦＳＩの数値的シミュレーションを終了するか否かを判断する判断ステップ５１８に移行する。ＦＳＩの数値的シミュレーションを終了しないと判断した場合、プロセス５００は、他のソリューションサイクルにおけるＦＳＩを取得するよう、ステップ５０８〜５１６に戻って、これらのステップを繰り返す。ＦＳＩの数値的シミュレーションを終了すると判断した場合、プロセス５００は終了する。限定するものではないが、終了条件には、予め定義された総シミュレーション時間への到達が含まれる。
図６Ａ〜図６Ｃは、本発明の一の実施形態にかかる、シミュレートされた構造体挙動（例えば構造体変形および新しいノード位置）に応じた空間−時間流体メッシュ調整の例示的な一連のシーケンスを示す概略図である。図６Ａにおいて、構造体６０２（楕円状の点線で示する）を表すＦＥＡモデルが、空間−時間流体メッシュ６１２ａ（図の簡単化のために二次元メッシュとして示す）に隣接している。ＦＥＡモデルおよび空間−時間流体メッシュ６１２ａは互いに重なり合う。
図６Ｂにおいて、変形した構造体６０４（楕円状の実線）は、ＦＳＩ境界面において受ける流体力に応じた構造体６０２のシミュレーション結果である。空間−時間流体メッシュ６１２ｂは、シミュレートされた構造体挙動から得られた新しいノード位置／速度を反映するよう、更新される。
最後に、図６Ｃにおいて、空間−時間流体６１２メッシュｃの内部ノードが、ＦＳＩ境界面における新しいノード位置に応じて、ユーザに選択されたメッシュ調整法を用いて、調整される。
一の側面において、本発明は、ここに説明した機能を実行可能な一つ以上のコンピュータシステムに対してなされたものである。コンピュータシステム７００の一例を、図７に示す。コンピュータシステム７００は、プロセッサ７０４など一つ以上のプロセッサを有する。プロセッサ７０４は、コンピュータシステム内部通信バス７０２に接続されている。種々のソフトウェアの実施形態を、この例示的なコンピュータシステムの点から説明する。この説明を読むと、いかにして、他のコンピュータシステムおよび／またはコンピューターアーキテクチャーを用いて、本発明を実行するかが、関連する技術分野に習熟している者には明らかになるであろう。
コンピュータシステム７００は、また、メインメモリ７０８好ましくはランダムアクセスメモリ（ＲＡＭ））を有しており、そして二次メモリ７１０を有することもできる。二次メモリ７１０は、例えば、一つ以上のハードディスクドライブ７１２、および／またはフレキシブルディスクドライブ、磁気テープドライブ、光ディスクドライブなどを表わす一つ以上のリムーバブルストレージドライブ７１４を有することができる。リムーバブルストレージドライブ７１４は、よく知られている方法で、リムーバブルストレージユニット７１８を読み取りおよび／またはリムーバブルストレージユニット７１８に書き込む。リムーバブルストレージユニット７１８は、リムーバブルストレージドライブ７１４によって読み取り・書き込みされるフレキシブルディスク、磁気テープ、光ディスクなどを表わす。以下にわかるように、リムーバブルストレージユニット７１８は、コンピューターソフトウェアおよび／またはデータを内部に記憶しているコンピュータで使用可能な記憶媒体を有している。
代替的な実施形態において、二次メモリ７１０は、コンピュータプログラムあるいは他の命令をコンピュータシステム７００にロードすることを可能にする他の同様な手段を有することもできる。そのような手段は、例えば、リムーバブルストレージユニット７２２とインタフェース７２０とを有することができる。そのようなものの例には、プログラムカートリッジおよびカートリッジのインタフェース（ビデオゲーム機に見られるようなものなど）と、リムーバブルメモリチップ（消去可能なプログラマブルＲＯＭ（ＥＰＲＯＭ）、ユニバーサルシリアルバス（ＵＳＢ）フラッシュメモリ、あるいはＰＲＯＭなど）および関連するソケットと、ソフトウェアおよびデータをリムーバブルストレージユニット７２２からコンピュータシステム７００に転送することを可能にする他のリムーバブルストレージユニット７２２およびインタフェース７２０と、が含まれうる。一般に、コンピュータシステム７００は、プロセススケジューリング、メモリ管理、ネットワーキングおよびＩ／Ｏサービスなどのタスクを行なうオペレーティングシステム（ＯＳ）ソフトウェアによって、制御され連係される。
通信用インタフェース７２４も、また、バス７０２に接続することができる。通信用インタフェース７２４は、ソフトウェアおよびデータをコンピュータシステム７００と外部装置との間で転送することを可能にする。通信用インタフェース７２４の例には、モデム、ネットワークインタフェイス（イーサネット（登録商標）・カードなど）、コミュニケーションポート、ＰＣＭＣＩＡ（Ｐｅｒｓｏｎａｌ Ｃｏｍｐｕｔｅｒ Ｍｅｍｏｒｙ Ｃａｒｄ Ｉｎｔｅｒｎａｔｉｏｎａｌ Ａｓｓｏｃｉａｔｉｏｎ）スロットおよびカードなど、が含まれうる。コンピュータ７００は、専用のセットの規則（つまりプロトコル）に基づいて、データネットワーク上の他の演算装置と通信する。一般的なプロトコルのうちの一つは、インターネットにおいて一般に用いられているＴＣＰ／ＩＰ（伝送コントロール・プロトコル／インターネット・プロトコル）である。一般に、通信インタフェース７２４は、データファイルをデータネットワーク上で伝達される小さいパケットへのアセンブリングを管理し、あるいは受信したパケット元のデータファイルへと再アセンブルする。さらに、通信インタフェース７２４は、正しい宛先に届くようそれぞれのパケットのアドレス部分に対処し、あるいはコンピュータ７００が宛先となっているパケットを他に向かわせることなく受信する。この書類において、「コンピュータプログラム媒体」および「コンピュータで使用可能な媒体」という用語は、リムーバブルストレージドライブ７１４および／またはハードディスクドライブ７１２に組み込まれたハードディスクなどの媒体を概ね意味して用いられている。これらのコンピュータプログラム製品は、コンピュータシステム７００にソフトウェアを提供する手段である。本発明は、このようなコンピュータプログラム製品に対してなされたものである。
コンピュータシステム７００は、また、コンピュータシステム７００をアクセスモニタ、キーボード、マウス、プリンタ、スキャナ、プロッタなどに提供する入出力（Ｉ／Ｏ）インタフェース７３０を有することができる。
コンピュータプログラム（コンピュータ制御ロジックともいう）は、メインメモリ７０８および／または二次メモリ７１０にアプリケーションモジュール７０６として記憶される。コンピュータプログラムを、通信用インタフェース７２４を介して受け取ることもできる。このようなコンピュータプログラムが実行された時、コンピュータプログラムによって、コンピュータシステム７００がここに説明した本発明の特徴を実行することが可能になる。詳細には、コンピュータプログラムが実行された時、コンピュータプログラムによって、プロセッサ７０４が本発明の特徴を実行することが可能になる。したがって、このようなコンピュータプログラムは、コンピュータシステム７００のコントローラを表わしている。
ソフトウェアを用いて発明が実行される実施形態において、ソフトウェアをコンピュータプログラム製品に記憶でき、リムーバブルストレージドライブ７１４、ハードドライブ７１２あるいは通信用インタフェース７２４を用いてコンピュータシステム７００へとロードすることができる。アプリケーションモジュール７０６は、プロセッサ７０４によって実行された時、アプリケーションモジュールによって、プロセッサ７０４がここに説明した本発明の機能を実行する。
所望のタスクを達成するために、Ｉ／Ｏインタフェース７３０を介したユーザ入力によってあるいはよることなしに、一つ以上のプロセッサ７０４によって実行することができる一つ以上のアプリケーションモジュール７０６を、メインメモリ７０８に、ロードすることもできる。動作においては、少なくとも一つのプロセッサ７０４がアプリケーションモジュール７０６のうちの一つが実行されると、結果が演算されて二次メモリ７１０（つまりハードディスクドライブ７１２）に記憶される。有限要素解析および／または空間−時間ＣＥ／ＳＥソルバーの状況は、テキストあるいはグラフィック表現で、Ｉ／Ｏインタフェース７３０を介してユーザに報告される。
本発明を具体的な実施形態を参照しながら説明したが、これらの実施形態は単なる例示であって、本発明を限定するものではない。開示した例示的な実施形態に対する種々の変更あるいは変形を、当業者は思いつくであろう。例えば、空間−時間流体メッシュを二次元（空間的次元が一つ）や三次元（空間的次元が二つ）として示したが、空間−時間流体メッシュを、図に示すことが容易ではない四次元（空間的次元が三つ）とすることもできる。つまり、発明の範囲は、ここで開示した具体的で例示的な実施形態に限定されず、当業者が容易に想到するあらゆる変更が、本願の精神および認識範囲そして添付の特許請求の範囲の権利範囲に含まれる。
１００ 有限要素解析モデル
１１０ａ〜１１０ｃ 構造体
１２０ａ〜１２０ｃ 流体領域
１３０ａ〜１３０ｃ外部表面
２００ メッシュ
２０４ メッシュ点
２１４ ソリューション要素
２２２ 基本保存要素
２２４ 合成保存要素
４００ 移動空間−時間メッシュ
４１０ 固定オイラーメッシュ
６０２ 構造体
６０４ 変形した構造体
６１２ａ〜６１２ｃ 空間−時間流体メッシュ
７０２ バス
７０４ プロセッサ
７０６ モジュール
７０８ メインメモリ（ＲＡＭ）
７１０ 二次メモリ
７１２ ハードディスクドライブ
７１４ リムーバブルストレージドライブ
７１８ リムーバブルストレージユニット
７２０ インタフェース
７２２ リムーバブルストレージユニット
７２４ 通信インタフェース
７３０ Ｉ／Ｏインタフェース

Claims (20)

  1. 移動流体を用いる空間−時間保存要素／ソリューション要素（ＣＥ／ＳＥ）ソルバーを用いて、流体構造体相互作用（ＦＳＩ）を数値的にシミュレートする方法であって、
    （ａ）有限要素解析（ＦＥＡ）モジュールおよびＣＥ／ＳＥソルバーアプリケーションモジュールがインストールされたコンピュータシステムに、流体領域定義と構造体定義とを受け取るステップであって、前記流体領域定義はＣＥ／ＳＥソルバーのために構成された空間−時間流体メッシュを備えており、構造体定義はＦＥＡモデルを備えている、ステップと、
    （ｂ）前記空間−時間流体メッシュと前記ＦＥＡモデルとを用いてＦＳＩ境界面を決定するステップと、
    （ｃ）前記ＦＳＩ境界面上で前記ＦＥＡモデルに作用する流体力を初期化するステップと、
    （ｄ）前記受け取った流体力に応じて前記ＦＥＡモデルを用いるＦＥＡを実行することにより、数値的にシミュレートされた構造体挙動を取得するステップであって、前記数値的にシミュレートされた構造体挙動が構造体の外部境界上のＦＥＡモデルのノード位置を含んでいるステップと、
    （ｅ）ＦＳＩ境界面における流体メッシュを更新して、前記構造体の外部境界のノード位置と整合させるステップと、
    （ｆ）ユーザに選択されたメッシュ調整法を用いて、前記ＦＳＩ境界面におけるノード位置に応じて前記流体メッシュの内部ノード位置を調整するステップと、
    （ｇ）前記調整された流体メッシュで、前記ＣＥ／ＳＥソルバーを用いて流体ソリューションを実行することにより、数値的にシミュレートされた流体挙動を取得するステップであって、該数値的にシミュレートされた流体挙動にはＦＳＩ境界面における流体力が含まれているステップと、
    （ｈ）現在のソリューション時間をインクリメントするステップと、
    （ｉ）終了条件が満たされるまでステップ（ｄ）〜（ｈ）を繰り返すステップと、
    を備える方法。
  2. 請求項１に記載の方法であって、流体領域が三つの空間的次元を含んでいる場合、前記空間−時間流体メッシュは四次元である、方法。
  3. 請求項１に記載の方法であって、前記ＦＳＩ境界面は構造体の外部境界の一部または全部と一致する、方法。
  4. 請求項１に記載の方法であって、前記空間−時間流体メッシュおよび前記ＦＥＡモデルは共通ノードまたは共通辺を共有する必要はない、方法。
  5. 請求項１に記載の方法であって、前記流体ソリューションの実行には、直前のソリューションサイクルの流体メッシュおよび現在のソリューションサイクルの前記調整された流体メッシュに基づいて流体領域の幾何学的パラメータを計算するステップがさらに含まれている、方法。
  6. 請求項５の方法であって、さらに、流体領域変数および対応する空間導関数を計算するステップを備える、方法。
  7. 請求項６に記載の方法であって、前記流体領域変数には、流体密度、圧力、速度、粘度等が含まれる、方法。
  8. 請求項１に記載の方法であって、流体領域が二つの空間的次元を含んでいる場合、前記空間−時間流体メッシュは三次元である、方法。
  9. 移動流体を用いる空間−時間保存要素／ソリューション要素（ＣＥ／ＳＥ）ソルバーを用いて、流体構造体相互作用（ＦＳＩ）を数値的にシミュレートするシステムであって、
    有限要素解析（ＦＥＡ）アプリケーションモジュールとＣＥ／ＳＥアプリケーションモジュールとに関するコンピュータ可読コードを記憶しているメインメモリと、
    前記メインメモリに接続される少なくとも１つのプロセッサであって、該少なくとも１つのプロセッサが前記メインメモリ内の前記コンピュータ可読コードを実行して、アプリケーションモジュールに、方法に基づいてオペレーションを実行させるシステムであって、該方法が、
    （ａ）ＣＥ／ＳＥソルバーのために構成された空間−時間流体メッシュを備えている流体領域定義と、ＦＥＡモデルを備えている構造体定義と、を受け取るステップと、
    （ｂ）前記空間−時間流体メッシュと前記ＦＥＡモデルとを用いてＦＳＩ境界面を決定するステップと、
    （ｃ）前記ＦＳＩ境界面上で前記ＦＥＡモデルに作用する流体力を初期化するステップと、
    （ｄ）前記受け取った流体力に応じて前記ＦＥＡモデルを用いるＦＥＡを実行することにより、数値的にシミュレートされた構造体挙動を取得するステップであって、前記数値的にシミュレートされた構造体挙動が構造体の外部境界上のＦＥＡモデルのノード位置を含んでいるステップと、
    （ｅ）ＦＳＩ境界面における流体メッシュを、前記構造体の外部境界のノード位置と整合させるよう、更新するステップと、
    （ｆ）ユーザに選択されたメッシュ調整法を用いて、前記ＦＳＩ境界面におけるノード位置に応じて前記流体メッシュの内部ノード位置を調整するステップと、
    （ｇ）前記調整された流体メッシュで、前記ＣＥ／ＳＥソルバーを用いて流体ソリューションを実行するにより、数値的にシミュレートされた流体挙動を取得するステップであって、該数値的にシミュレートされた流体挙動にはＦＳＩ境界面における流体力が含まれているステップと、
    （ｈ）現在のソリューション時間をインクリメントするステップと、
    （ｉ）終了条件が満たされるまでステップ（ｄ）〜（ｈ）を繰り返すステップと、
    を備えている、システム。
  10. 請求項９に記載のシステムであって、流体領域が三つの空間的次元を含んでいる場合、前記空間−時間流体メッシュは四次元である、システム。
  11. 請求項９に記載のシステムであって、流体領域が二つの空間的次元を含んでいる場合、前記空間−時間流体メッシュは三次元である、システム。
  12. 請求項９に記載のシステムであって、前記流体ソリューションの実行には、直前のソリューションサイクルの流体メッシュおよび現在のソリューションサイクルの前記調整された流体メッシュに基づいて流体領域の幾何学的パラメータを計算するステップがさらに含まれている、システム。
  13. 請求項１２に記載の前記システムであって、記方法がさらに、流体領域変数および対応する空間導関数を計算するステップを備えている、システム。
  14. 請求項１３に記載のシステムであって、前記流体領域変数には、流体密度、圧力、速度、粘度等が含まれる、システム。
  15. コンピュータシステムにおいて実行されたとき、方法に基づいて、移動流体メッシュを用いる空間−時間保存要素／ソリューション要素（ＣＥ／ＳＥ）ソルバーを用いて、流体構造体相互作用（ＦＳＩ）を数値的にシミュレートする命令を備えるコンピュータが読み取り可能な非一時的記憶媒体であって、該方法が、
    （ａ）ＦＥＡモジュールおよびＣＥ／ＳＥソルバーアプリケーションモジュールがインストールコンピュータシステムに、ＣＥ／ＳＥソルバーのために構成された空間−時間流体メッシュを備えている流体領域定義と、ＦＥＡモデルを備えている構造体定義と、を受け取るステップと、
    （ｂ）前記空間−時間流体メッシュと前記ＦＥＡモデルとを用いてＦＳＩ境界面を決定するステップと、
    （ｃ）前記ＦＳＩ境界面上で前記ＦＥＡモデルに作用する流体力を初期化するステップと、
    （ｄ）前記受け取った流体力に応じて前記ＦＥＡモデルを用いるＦＥＡを実行することにより、数値的にシミュレートされた構造体挙動を取得するステップであって、前記数値的にシミュレートされた構造体挙動が構造体の外部境界上のＦＥＡモデルのノード位置を含んでいるステップと、
    （ｅ）ＦＳＩ境界面における流体メッシュを、前記構造体の外部境界のノード位置と整合させるよう、更新するステップと、
    （ｆ）ユーザに選択されたメッシュ調整法を用いて、前記ＦＳＩ境界面におけるノード位置に応じて前記流体メッシュの内部ノード位置を調整するステップと、
    （ｇ）前記調整された流体メッシュで、前記ＣＥ／ＳＥソルバーを用いて流体ソリューションを実行するにより、数値的にシミュレートされた流体挙動を取得するステップであって、該数値的にシミュレートされた流体挙動にはＦＳＩ境界面における流体力が含まれているステップと、
    （ｈ）現在のソリューション時間をインクリメントするステップと、
    （ｉ）終了条件が満たされるまでステップ（ｄ）〜（ｈ）を繰り返すステップと、
    を備えている、コンピュータが読み取り可能な非一時的記憶媒体。
  16. 請求項１５に記載のコンピュータが読み取り可能な非一時的記憶媒体であって、流体領域が三つの空間的次元を含んでいる場合、前記空間−時間流体メッシュは四次元であ、るコンピュータが読み取り可能な非一時的記憶媒体。
  17. 請求項１５に記載のコンピュータが読み取り可能な非一時的記憶媒体であって、流体領域が二つの空間的次元を含んでいる場合、前記空間−時間流体メッシュは三次元である、コンピュータが読み取り可能な非一時的記憶媒体。
  18. 請求項１５に記載のコンピュータが読み取り可能な非一時的記憶媒体であって、前記流体ソリューションの実行には、直前のソリューションサイクルの流体メッシュおよび現在のソリューションサイクルの前記調整された流体メッシュに基づいて流体領域の幾何学的パラメータを計算するステップがさらに含まれている、コンピュータが読み取り可能な非一時的記憶媒体。
  19. 請求項１８に記載のコンピュータが読み取り可能な非一時的記憶媒体であって、前記方法がさらに、流体領域変数および対応する空間導関数を計算するステップを備えている、コンピュータが読み取り可能な非一時的記憶媒体。
  20. 請求項１９に記載のコンピュータが読み取り可能な非一時的記憶媒体であって、前記流体領域変数には、流体密度、圧力、速度、粘度等が含まれる、コンピュータが読み取り可能な非一時的記憶媒体。
JP2014036188A 2013-03-05 2014-02-27 流体領域に対して移動メッシュを用いる空間−時間ｃｅ／ｓｅソルバーを利用したｆｓｉの数値的シミュレーション Pending JP2014174990A (ja)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
US13/786,366 2013-03-05
US13/786,366 US20140257765A1 (en) 2013-03-05 2013-03-05 Numerical Simulation of FSI Using The Space-Time CE/SE Solver With A Moving Mesh For The Fluid Domain

Publications (1)

Publication Number Publication Date
JP2014174990A true JP2014174990A (ja) 2014-09-22

Family

ID=51466834

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2014036188A Pending JP2014174990A (ja) 2013-03-05 2014-02-27 流体領域に対して移動メッシュを用いる空間−時間ｃｅ／ｓｅソルバーを利用したｆｓｉの数値的シミュレーション

Country Status (4)

Country Link
US (1) US20140257765A1 (ja)
JP (1) JP2014174990A (ja)
KR (1) KR20140109264A (ja)
CN (1) CN104036062A (ja)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20160051323A (ko) * 2014-11-03 2016-05-11 한국항공우주산업 주식회사 3차원 압력 사상 방법, 이를 구현하기 위한 프로그램이 저장된 기록매체 및 이를 구현하기 위해 매체에 저장된 컴퓨터프로그램
KR20160051349A (ko) * 2014-11-03 2016-05-11 한국항공우주산업 주식회사 항공기 연료탱크의 내부 압력 사상 방법, 이를 구현하기 위한 프로그램이 저장된 기록매체 및 이를 구현하기 위해 매체에 저장된 컴퓨터프로그램

Families Citing this family (18)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP6300244B2 (ja) * 2014-07-03 2018-03-28 富士通株式会社 生体シミュレーション装置、生体シミュレーション装置の制御方法、および生体シミュレーション装置の制御プログラム
US10311181B2 (en) * 2015-05-09 2019-06-04 Livermore Software Technology Corp. Methods and systems for creating computerized mesh model for layered shell-like structure
CN105205225B (zh) * 2015-08-28 2018-06-19 大连理工大学 一种用于浮式海洋结构物静水平衡状态计算的有限元方法
US9798841B2 (en) * 2015-09-15 2017-10-24 Livermore Software Technology Corp. Systems and methods of conducting numerical simulation of an underwater explosion
CN105205269B (zh) * 2015-09-25 2018-05-15 宝鸡石油钢管有限责任公司 一种预测特殊螺纹接头油套管流固耦合特性的分析方法
CN105677969B (zh) * 2016-01-06 2019-02-12 北京航空航天大学 面向单机实时仿真的复杂流体网络稳定性建模方法
JP6983823B2 (ja) * 2016-06-16 2021-12-17 タタ・コンサルタンシー・サーヴィシズ・リミテッド 調節される空間の熱流体管理のためのシステムおよび方法
CN107895069B (zh) * 2017-10-30 2020-10-27 北京理工大学 一种基于复合材料结构的流固耦合数值预测方法
US10891788B2 (en) * 2017-12-13 2021-01-12 Dassault Systemes Simulia Corp. Systems and methods for finite element mesh repair
CN109446471B (zh) * 2018-09-07 2021-07-23 北京航空航天大学 一种考虑载荷不确定性的流固耦合界面数据传递方法
CN110110479B (zh) * 2019-05-20 2023-04-07 合肥美的电冰箱有限公司 制冷设备的温度场仿真方法、终端及可读存储介质
CN111159959B (zh) * 2019-12-12 2023-05-02 智慧航海（青岛）科技有限公司 一种基于Simulink和Fluent的智能船舶运动仿真方法
CN111881600B (zh) * 2020-06-29 2024-03-19 郑州大学 基于MpCCI多物理场作用的管道耦合分析方法
CN112163381B (zh) * 2020-09-25 2024-05-14 西安热工研究院有限公司 一种适用于复杂地形风场流动数值模拟的侧向边界条件设置方法
CN114595647B (zh) * 2022-05-10 2022-08-02 中国空气动力研究与发展中心设备设计与测试技术研究所 一种磁浮飞行风洞气动结构耦合仿真评估方法
CN116341419B (zh) * 2023-05-17 2023-08-01 中国科学院、水利部成都山地灾害与环境研究所 一种流固耦合的数值确定方法、系统及电子设备
CN116484771B (zh) * 2023-06-21 2023-08-25 陕西空天信息技术有限公司 一种轴流压气机cfd网格生成的方法及装置
CN117829026B (zh) * 2024-01-05 2024-06-18 江苏力磁医疗设备有限公司 磁共振兼容转运培养箱流体多物理场耦合仿真计算方法

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20160051323A (ko) * 2014-11-03 2016-05-11 한국항공우주산업 주식회사 3차원 압력 사상 방법, 이를 구현하기 위한 프로그램이 저장된 기록매체 및 이를 구현하기 위해 매체에 저장된 컴퓨터프로그램
KR20160051349A (ko) * 2014-11-03 2016-05-11 한국항공우주산업 주식회사 항공기 연료탱크의 내부 압력 사상 방법, 이를 구현하기 위한 프로그램이 저장된 기록매체 및 이를 구현하기 위해 매체에 저장된 컴퓨터프로그램
KR101682598B1 (ko) 2014-11-03 2016-12-06 한국항공우주산업 주식회사 항공기 연료탱크의 내부 압력 사상 방법, 이를 구현하기 위한 프로그램이 저장된 기록매체 및 이를 구현하기 위해 매체에 저장된 컴퓨터프로그램
KR101682379B1 (ko) 2014-11-03 2016-12-07 한국항공우주산업 주식회사 3차원 압력 사상 방법, 이를 구현하기 위한 프로그램이 저장된 기록매체 및 이를 구현하기 위해 매체에 저장된 컴퓨터프로그램

Also Published As

Publication number Publication date
KR20140109264A (ko) 2014-09-15
US20140257765A1 (en) 2014-09-11
CN104036062A (zh) 2014-09-10

Similar Documents

Publication Publication Date Title
JP2014174990A (ja) 流体領域に対して移動メッシュを用いる空間−時間ｃｅ／ｓｅソルバーを利用したｆｓｉの数値的シミュレーション
JP6043198B2 (ja) 任意ラグランジュ・オイラー法（ａｌｅ）における要素細分化法およびシステム
JP5330300B2 (ja) 有限要素解析法における爆発シミュレーション
Michler Aircraft control surface deflection using RBF‐based mesh deformation
Wang et al. Algorithms for interface treatment and load computation in embedded boundary methods for fluid and fluid–structure interaction problems
US8612186B2 (en) Numerical simulation of structural behaviors using a meshfree-enriched finite element method
US8306793B2 (en) Systems and methods of performing vibro-acoustic analysis of a structure
US8855977B2 (en) Numerically simulating structural behaviors of a product using explicit finite element analysis with a combined technique of mass scaling and subcycling
KR102544597B1 (ko) 유한 요소 분석법으로 보강된 콘크리트의 구조적 거동을 시뮬레이션하는 방법 및 시스템
Sharma et al. Overset meshes for incompressible flows: On preserving accuracy of underlying discretizations
US20180239848A1 (en) Numerical Blast Simulation Methods and Systems Thereof
US8855976B2 (en) Numerically simulating structural behaviors of a product using explicit finite element analysis with a mass scaling enhanced subcycling technique
JP5589198B2 (ja) 有限要素法において８ノード六面体エレメントの剪断ロッキングを低減する方法
US20170255724A1 (en) Enhanced Global Design Variables Used In Structural Topology Optimization Of A Product In An Impact Event
Zheng et al. An improved local remeshing algorithm for moving boundary problems
US20110191068A1 (en) Multiscale substructures in finite element analysis
KR102459848B1 (ko) 입자에 기반하여 대상 객체를 고속으로 모델링하는 방법 및 장치
CN106611078B (zh) 产品的具有时间步长控制方案的高效显式有限元分析
Rendall et al. Conservative unsteady aerodynamic simulation of arbitrary boundary motion using structured and unstructured meshes in time
McDaniel et al. Efficient mesh deformation for computational stability and control analyses on unstructured viscous meshes
US9507892B2 (en) Methods and systems for using bi-directional level sets to partition an undirected graph representing a matrix to be used in CAE
US8676551B1 (en) Multi-solver simulation of dynamic systems in a modeling environment
Martineau et al. A mesh movement algorithm for high quality generalised meshes
JP2016212862A (ja) 層状の殻状構造に関する演算メッシュモデルを生成する方法およびシステム
KR101110342B1 (ko) 유체 시뮬레이션 형상 제어 시스템 및 방법