(57)【要約】【課題】反復計算により解析解を安定して迅速に収束させることができる解析装置等を提供する。【解決手段】入力パラメータを更新しながら漸化式を用いて反復計算を行うことで解析解を収束させる流体解析モデルを用いた流体解析装置において、更新後の入力パラメータＩκ＋１は、不足緩和係数βＩを含む所定の算出式で与えられ、不足緩和係数βＩは、更新後の入力パラメータＩκ＋１と更新前の入力パラメータＩκとの差分が大きい差分大条件において値が小さく、差分が差分大条件に比して小さい差分小条件において値が大きく、且つ、０より大きく１より小さい値であって、微分可能に連続する関数となっており、制御部は、入力パラメータＩκ＋１と入力パラメータＩκとの差分に基づいて、関数から、不足緩和係数βＩの値を導出し、導出した不足緩和係数βＩに基づいて、所定の算出式から、更新後の入力パラメータＩκ＋１を導出する。【選択図】図４
2019191905.tif 000002
【特許請求の範囲】
【請求項１】
入力パラメータに基づき解析モデルを用いて解析することで解析解を導出する演算部を備える解析装置において、
前記解析モデルは、前記入力パラメータを更新しながら漸化式を用いて反復計算を行うことで前記解析解を収束させるモデルとなっており、
更新後の前記入力パラメータをＩκ＋１とし、更新前の前記入力パラメータをＩκとし、前記解析解を収束させるための不足緩和係数をβＩとすると、
更新後の前記入力パラメータＩκ＋１は、前記不足緩和係数βＩを含む算出式で与えられ、
前記入力パラメータＩκ＋１と前記入力パラメータＩκとの差分は、差分大条件と、差分小条件とを含み、
前記不足緩和係数βＩは、
前記入力パラメータＩκ＋１と前記入力パラメータＩκとの差分が前記差分小条件に比して大きい前記差分大条件において値が小さく、
前記入力パラメータＩκ＋１と前記入力パラメータＩκとの差分が前記差分大条件に比して小さい前記差分小条件において値が大きく、
且つ、０より大きく１より小さい値であって、微分可能に連続する関数となっており、
前記演算部は、
前記入力パラメータＩκ＋１と前記入力パラメータＩκとの差分に基づいて、前記関数から、前記不足緩和係数βＩの値を導出し、導出した前記不足緩和係数βＩに基づいて、前記算出式から、更新後の前記入力パラメータＩκ＋１を導出することを特徴とする解析装置。
【請求項２】
前記関数は、下記する式（１）となるロジスティック関数であり、
Ａ，Ｋ，Ｃは、係数であることを特徴とする請求項１に記載の解析装置。
【数１】
2019191905.tif 000006
【請求項３】
前記解析モデルは、気相の流体と液相の流体とからなる二流体の流体解析モデルであり、
前記入力パラメータは、前記流体の圧力と前記流体の速度とを含み、
前記流体解析モデルは、前記流体の圧力から前記流体の速度を導出し、前記流体の速度から前記流体の圧力を導出して、前記流体の圧力と前記流体の速度を更新しながら反復計算を行うモデルであることを特徴とする請求項１または２に記載の解析装置。
【請求項４】
入力パラメータに基づき解析モデルを用いて解析することで解析解を導出する解析装置により実行される解析方法において、
前記解析モデルは、前記入力パラメータを更新しながら漸化式を用いて反復計算を行うことで前記解析解を収束させるモデルとなっており、
更新後の前記入力パラメータをＩκ＋１とし、更新前の前記入力パラメータをＩκとし、前記解析解を収束させるための不足緩和係数をβＩとすると、
更新後の前記入力パラメータＩκ＋１は、前記不足緩和係数βＩを含む算出式で与えられ、
前記入力パラメータＩκ＋１と前記入力パラメータＩκとの差分は、差分大条件と、差分小条件とを含み、
前記不足緩和係数βＩは、
前記入力パラメータＩκ＋１と前記入力パラメータＩκとの差分が前記差分小条件に比して大きい前記差分大条件において値が小さく、
前記入力パラメータＩκ＋１と前記入力パラメータＩκとの差分が前記差分大条件に比して小さい前記差分小条件において値が大きく、
且つ、０より大きく１より小さい値であって、微分可能に連続する関数となっており、
前記入力パラメータＩκ＋１と前記入力パラメータＩκとの差分を導出する差分導出ステップと、
導出した前記差分に基づいて、前記関数から、前記不足緩和係数βＩの値を導出する不足緩和係数導出ステップと、
導出した前記不足緩和係数βＩに基づいて、前記算出式から、更新後の前記入力パラメータＩκ＋１を導出する入力パラメータ導出ステップと、を前記解析装置に実行させることを特徴とする解析方法。
【請求項５】
入力パラメータに基づき解析モデルを用いて解析することで解析解を導出する解析装置により実行される解析プログラムにおいて、
前記解析モデルは、前記入力パラメータを更新しながら漸化式を用いて反復計算を行うことで前記解析解を収束させるモデルとなっており、
更新後の前記入力パラメータをＩκ＋１とし、更新前の前記入力パラメータをＩκとし、前記解析解を収束させるための不足緩和係数をβＩとすると、
更新後の前記入力パラメータＩκ＋１は、前記不足緩和係数βＩを含む算出式で与えられ、
前記入力パラメータＩκ＋１と前記入力パラメータＩκとの差分は、差分大条件と、差分小条件とを含み、
前記不足緩和係数βＩは、
前記入力パラメータＩκ＋１と前記入力パラメータＩκとの差分が前記差分小条件に比して大きい前記差分大条件において値が小さく、
前記入力パラメータＩκ＋１と前記入力パラメータＩκとの差分が前記差分大条件に比して小さい前記差分小条件において値が大きく、
且つ、０より大きく１より小さい値であって、微分可能に連続する関数となっており、
前記入力パラメータＩκ＋１と前記入力パラメータＩκとの差分を導出する差分導出ステップと、
導出した前記差分に基づいて、前記関数から、前記不足緩和係数βＩの値を導出する不足緩和係数導出ステップと、
導出した前記不足緩和係数βＩに基づいて、前記算出式から、更新後の前記入力パラメータＩκ＋１を導出する入力パラメータ導出ステップと、を前記解析装置に実行させることを特徴とする解析プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、解析装置、解析方法及び解析プログラムに関するものである。
【背景技術】
【０００２】
従来、解析装置として、気相の流体と液相の流体とからなる二流体の流体解析モデルを用いた流体解析装置がある。このような流体解析装置として、入力されたパラメータに基づいて気液二相流の流動挙動を解析する気液二相流解析装置が知られている（例えば、特許文献１）。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】特開２０１１−１５０６８０号公報
【発明の概要】
【発明が解決しようとする課題】
【０００４】
ところで、二流体の流体解析モデルでは、入力される入力パラメータとして、圧力と流速と含み、圧力と流速とを反復して計算することで、解析解を算出する計算手法が利用されている。二流体モデルにおいて、圧力と流速を反復して計算する場合、解析解は、発散する可能性がある。このため、一般的に、解析解を収束させて安定的に解析するために、二流体の流体解析モデルには、不足緩和を導入している。この不足緩和は、従来において、その値を経験則で与えていた。しかしながら、経験則で与えられる不足緩和の値は最適値ではない。つまり、不足緩和の値が小さいと、詳細な進展を計算することになるため、計算回数が多くなり、計算時間が長くなってしまう。一方で、計算時間を短くしようとして、不足緩和の値を大きくすると、解析解が発散してしまう可能性がある。
【０００５】
そこで、本発明は、反復計算により解析解を安定して迅速に収束させることができる解析装置、解析方法及び解析プログラムを提供することを課題とする。
【課題を解決するための手段】
【０００６】
本発明の解析装置は、入力パラメータに基づき解析モデルを用いて解析することで解析解を導出する演算部を備える解析装置において、前記解析モデルは、前記入力パラメータを更新しながら漸化式を用いて反復計算を行うことで前記解析解を収束させるモデルとなっており、更新後の前記入力パラメータをＩκ＋１とし、更新前の前記入力パラメータをＩκとし、前記解析解を収束させるための不足緩和係数をβＩとすると、更新後の前記入力パラメータＩκ＋１は、前記不足緩和係数βＩを含む算出式で与えられ、前記入力パラメータＩκ＋１と前記入力パラメータＩκとの差分は、差分大条件と、差分小条件とを含み、前記不足緩和係数βＩは、前記入力パラメータＩκ＋１と前記入力パラメータＩκとの差分が前記差分小条件に比して大きい前記差分大条件において値が小さく、前記入力パラメータＩκ＋１と前記入力パラメータＩκとの差分が前記差分大条件に比して小さい前記差分小条件において値が大きく、且つ、０より大きく１より小さい値であって、微分可能に連続する関数となっており、前記演算部は、前記入力パラメータＩκ＋１と前記入力パラメータＩκとの差分に基づいて、前記関数から、前記不足緩和係数βＩの値を導出し、導出した前記不足緩和係数βＩに基づいて、前記算出式から、更新後の前記入力パラメータＩκ＋１を導出することを特徴とする。
【０００７】
また、本発明の解析方法は、入力パラメータに基づき解析モデルを用いて解析することで解析解を導出する解析装置により実行される解析方法において、前記解析モデルは、前記入力パラメータを更新しながら漸化式を用いて反復計算を行うことで前記解析解を収束させるモデルとなっており、更新後の前記入力パラメータをＩκ＋１とし、更新前の前記入力パラメータをＩκとし、前記解析解を収束させるための不足緩和係数をβＩとすると、更新後の前記入力パラメータＩκ＋１は、前記不足緩和係数βＩを含む算出式で与えられ、前記入力パラメータＩκ＋１と前記入力パラメータＩκとの差分は、差分大条件と、差分小条件とを含み、前記不足緩和係数βＩは、前記入力パラメータＩκ＋１と前記入力パラメータＩκとの差分が前記差分小条件に比して大きい前記差分大条件において値が小さく、前記入力パラメータＩκ＋１と前記入力パラメータＩκとの差分が前記差分大条件に比して小さい前記差分小条件において値が大きく、且つ、０より大きく１より小さい値であって、微分可能に連続する関数となっており、前記入力パラメータＩκ＋１と前記入力パラメータＩκとの差分を導出する差分導出ステップと、導出した前記差分に基づいて、前記関数から、前記不足緩和係数βＩの値を導出する不足緩和係数導出ステップと、導出した前記不足緩和係数βＩに基づいて、前記算出式から、更新後の前記入力パラメータＩκ＋１を導出する入力パラメータ導出ステップと、を前記解析装置に実行させることを特徴とする。
【０００８】
また、本発明の解析プログラムは、入力パラメータに基づき解析モデルを用いて解析することで解析解を導出する解析装置により実行される解析プログラムにおいて、前記解析モデルは、前記入力パラメータを更新しながら漸化式を用いて反復計算を行うことで前記解析解を収束させるモデルとなっており、更新後の前記入力パラメータをＩκ＋１とし、更新前の前記入力パラメータをＩκとし、前記解析解を収束させるための不足緩和係数をβＩとすると、更新後の前記入力パラメータＩκ＋１は、前記不足緩和係数βＩを含む算出式で与えられ、前記入力パラメータＩκ＋１と前記入力パラメータＩκとの差分は、差分大条件と、差分小条件とを含み、前記不足緩和係数βＩは、前記入力パラメータＩκ＋１と前記入力パラメータＩκとの差分が前記差分小条件に比して大きい前記差分大条件において値が小さく、前記入力パラメータＩκ＋１と前記入力パラメータＩκとの差分が前記差分大条件に比して小さい前記差分小条件において値が大きく、且つ、０より大きく１より小さい値であって、微分可能に連続する関数となっており、前記入力パラメータＩκ＋１と前記入力パラメータＩκとの差分を導出する差分導出ステップと、導出した前記差分に基づいて、前記関数から、前記不足緩和係数βＩの値を導出する不足緩和係数導出ステップと、導出した前記不足緩和係数βＩに基づいて、前記算出式から、更新後の前記入力パラメータＩκ＋１を導出する入力パラメータ導出ステップと、を前記解析装置に実行させることを特徴とする。
【０００９】
これらの構成によれば、入力パラメータＩκ＋１と入力パラメータＩκとの差分が大きく、解析解が発散し易い差分大条件において、不足緩和係数βＩの値を小さくすることができる。このため、修正値Ｉ’を小さくできることから、解析解の発散を抑制することができ、安定的な解析を実行することができる。一方で、入力パラメータＩκ＋１と入力パラメータＩκとの差分が小さく、解析解が発散し難い差分小条件において、不足緩和係数βＩの値が大きくすることができる。このため、修正値Ｉ’を大きくできることから、解析解の収束を早めることができ、迅速な解析を実行することができる。このため、反復計算により解析解を安定して迅速に収束させることができる。
【００１０】
また、前記関数は、下記する式（１）となるロジスティック関数であり、Ａ，Ｋ，Ｃは、係数であることが、好ましい。
【００１１】
【数１】
2019191905.tif 000003
【００１２】
この構成によれば、不足緩和係数βＩを、取り扱いやすい関数を用いて表現することができる。つまり、係数Ａ，Ｋ，Ｃを適宜設定することで、関数を容易に変化させることができるため、例えば、入力パラメータの種類に応じた適切な関数を容易に設定することが可能となる。
【００１３】
また、前記解析モデルは、気相の流体と液相の流体とからなる二流体の流体解析モデルであり、前記入力パラメータは、前記流体の圧力と前記流体の速度とを含み、前記流体解析モデルは、前記流体の圧力から前記流体の速度を導出し、前記流体の速度から前記流体の圧力を導出して、前記流体の圧力と前記流体の速度を更新しながら反復計算を行うモデルであることが、好ましい。
【００１４】
この構成によれば、二流体の流体解析モデルに適用することができ、二流体の流体解析の解析解を安定して迅速に収束させることができる。
【図面の簡単な説明】
【００１５】
【図１】図１は、本実施形態の流体解析装置を模式的に表した概略構成図である。
【図２】図２は、不足緩和係数の変化を示すグラフである。
【図３】図３は、二流体の流体解析に関するフローチャートである。
【図４】図４は、不足緩和係数の導出に関するフローチャートである。
【発明を実施するための形態】
【００１６】
以下に、本発明に係る実施形態を図面に基づいて詳細に説明する。なお、この実施形態によりこの発明が限定されるものではない。また、下記実施形態における構成要素には、当業者が置換可能かつ容易なもの、あるいは実質的に同一のものが含まれる。さらに、以下に記載した構成要素は適宜組み合わせることが可能であり、また、実施形態が複数ある場合には、各実施形態を組み合わせることも可能である。
【００１７】
［本実施形態］
本実施形態の解析方法及び解析装置は、原子炉内の冷却材の挙動を解析する流体解析方法及び流体解析装置１０である。原子炉内の冷却材は、事故時において液相と気相とが混在する二流体となる場合がある。このため、本実施形態の流体解析では、二流体となる冷却材の流動を解析している。なお、本実施形態では、流体解析方法及び流体解析装置１０に適用して説明するが、下記する解析を実行可能な方法及び装置であれば、いずれの解析方法及び解析装置であってもよく、特に限定されない。
【００１８】
図１を参照して、気相の流体と液相の流体とからなる二流体の流体解析を実行する流体解析装置１０について説明する。図１は、本実施形態の流体解析装置を模式的に表した概略構成図である。流体解析装置１０は、二流体の流速、圧力、ボイド率、エンタルピ及び密度等の入力パラメータに基づいて、流体解析モデルを用いて解析することで、解析解を導出している。
【００１９】
流体解析装置１０は、各種プログラムを実行して各種処理を実行可能な制御部（演算部）１１と、各種プログラムおよびデータを記憶する記憶部１２と、キーボード等の入力デバイスで構成された入力部１３と、モニタ等の表示デバイスで構成された表示部１４とを有している。なお、流体解析装置１０は、単体の装置で構成してもよいし、他の装置と一体に構成してもよいし、演算装置及びデータサーバ等の各種装置を組み合わせたシステムとして構成してもよく、特に限定されない。
【００２０】
記憶部１２には、各種プログラムとして、例えば、流体解析に用いられる流体解析モデルを含む流体解析プログラムＰが記憶されている。また、記憶部１２には、データとして、例えば、上記した入力パラメータを含む各種データが記憶されている。
【００２１】
流体解析装置１０は、流体解析プログラムＰを実行すると、入力部１３から入力された入力データ及び記憶部１２に記憶されたデータ等を制御部１１が取得する。制御部１１は、取得した入力パラメータを含む各種データを入力パラメータとし、この入力パラメータに基づいて流体解析モデルを用いて計算を行うことで、解析解を算出する。
【００２２】
次に、流体解析モデルについて具体的に説明する。本実施形態の流体解析モデルは、入力パラメータを更新しながら漸化式を用いて反復計算を行うことで解析解を収束させるモデルとなっており、いわゆるＩＰＳＡ法（Inter  Phase  Slip  Algorithm）と呼ばれる流体解析法を用いている。ＩＰＳＡ法は、気相及び液相が混在する二流体の流体解析を行う手法であり、反復法により気相及び液相の連続の式を満足する運動方程式を解いている。ここで、図３を参照し、流体解析の演算フローについて説明する。図３は、二流体の流体解析に関するフローチャートである。
【００２３】
ここで、解析対象となる流体が存在する解析空間には、解析空間を離散化するために、多数の格子点が設定されている。先ず、制御部１１は、各格子点に対して、入力パラメータの初期値を設定する（ステップＳ１１）。具体的に、ステップＳ１１では、各格子点における流体の状態量及び物性値、流速、圧力に対して初期値を設定する。
【００２４】
次に、制御部１１は、気相及び液相の連続の式に、仮想値となる流速、密度等の入力パラメータを設定する（ステップＳ１２）。ステップＳ１２において、制御部１１は、所定のボイド率となる流体に対して、つまり、ボイド率が０となる気相状態の流体、及びボイド率が１となる液相状態の流体に対して、仮想値となる流速、密度等の入力パラメータを与える。
【００２５】
次に、制御部１１は、流速の推定値を算出する（ステップＳ１３）。ステップＳ１３において、制御部１１は、最新となる入力パラメータに基づいて、所定の算出式から、各格子点における流速の推定値を算出する。続いて、制御部１１は、算出した流速の推定値を用いて圧力場を修正する（ステップＳ１４）。ステップＳ１４において、制御部１１は、算出した流速の推定値を含む入力パラメータに基づいて、所定の算出式から、圧力Ｐ＊を算出すると共に、圧力の修正量Ｐ’  を算出し、下記する（１）式から、修正後の圧力を算出する。ここで、Ｐκ＋１は、修正後（更新後）の圧力であり、Ｐ＊は、流速の推定値から算出される圧力であり、β（Ｉ＝Ｐ）は、解析解を収束させるための不足緩和係数であり、Ｐ’  は、流速の推定値から算出される圧力の修正量である。
【００２６】
【数２】
2019191905.tif 000004
【００２７】
また、制御部１１は、ボイド率を更新する（ステップＳ１５）。ステップＳ１５において、制御部１１は、所定の算出式から、気相の流体と液相の流体との体積率を算出し、算出した体積率からボイド率を算出すると共に、不足緩和係数β（Ｉ＝α）を考慮して、修正後（更新後）のボイド率を算出する。
【００２８】
さらに、制御部１１は、エンタルピ及び密度を算出する（ステップＳ１６）。ステップＳ１６において、制御部１１は、不足緩和係数β（Ｉ＝ｈ）を考慮した所定の算出式から、修正後（更新後）のエンタルピを算出する。また、制御部１１は、算出したエンタルピに基づいて、所定の算出式から、密度を算出する。
【００２９】
続いて、制御部１１は、解析解が収束したか否かを判定する（ステップＳ１７）。制御部１１は、解析解が収束していないと判定する（ステップＳ１７：Ｎｏ）と、上記ステップにおいて更新された入力パラメータを、反復計算のための入力パラメータとして設定する（ステップＳ１８）。そして、制御部１１は、更新後の入力パラメータに基づいて、ステップＳ１２からステップＳ１６を実行し、解析解が収束するまで反復計算を実行する（ステップＳ１９）。一方で、制御部１１は、解析解が収束したと判定する（ステップＳ１７：Ｙｅｓ）と、流体解析の実行を終了する。
【００３０】
ここで、上記した圧力、ボイド率及びエンタルピの不足緩和係数βＩ（＝ｐ，α，ｈ）は、下記する（２）式で示すロジスティック関数で表現されている。ここで、Ｉ（＝Ｐ，α，ｈ）κ＋１は、更新後の入力パラメータとなっており、Ｉ（＝Ｐ，α，ｈ）κは、更新前の入力パラメータとなっている。また、Ａ，Ｃ，Ｋは、係数となっている。
【００３１】
【数３】
2019191905.tif 000005
【００３２】
ロジスティック関数となる不足緩和係数βＩは、例えば、入力パラメータが圧力Ｐである場合、図２に示すような関数となっている。ここで、図２は、その横軸が、入力パラメータＰκ＋１と入力パラメータＰκとの差分（｜Ｐκ＋１−Ｐκ｜）となっており、その縦軸が、不足緩和係数βｐとなっている。図２に示すように、不足緩和係数βｐは、入力パラメータＰκ＋１と入力パラメータＰκとの差分が大きい差分大条件において値が小さく、入力パラメータＰκ＋１と入力パラメータＰκとの差分が差分大条件に比して小さい差分小条件において値が大きくなっている。また、不足緩和係数βｐは、０より大きく１より小さい値であって、微分可能に連続する関数となっている。このため、差分大条件では、不足緩和係数βｐの値を小さくすることで、更新前の入力パラメータＰκに対する入力パラメータの修正量Ｐ’を小さくして、更新後の入力パラメータＰκ＋１を算出している。一方で、差分小条件では、不足緩和係数βｐの値を大きくすることで、更新前の入力パラメータＰκに対する入力パラメータの修正量Ｐ’を大きくして、更新後の入力パラメータＰκ＋１を算出している。
【００３３】
なお、図２では、圧力に関する不足緩和係数βｐについて説明したが、ボイド率及びエンタルピについても同様のロジスティック関数となっているため、説明を省略する。このとき、係数Ａ，Ｃ，Ｋは、任意の値となっており、入力パラメータの種類、つまり、圧力、ボイド率及びエンタルピに応じた適切な修正量となるように設定されている。
【００３４】
次に、図４を参照して、不足緩和係数βｐを用いた更新後の入力パラメータの算出について説明する。制御部１１は、入力パラメータＩκ＋１と入力パラメータＩκとの差分を算出する（ステップＳ２１：差分導出ステップ）。続いて、制御部１１は、算出した差分に基づいて、（２）式のロジスティック関数から、不足緩和係数βＩの値を算出する（ステップＳ２２：不足緩和係数導出ステップ）。制御部１１は、不足緩和係数βＩを含む所定の算出式から、更新後の入力パラメータＩκ＋１を算出する（ステップＳ２３：入力パラメータ導出ステップ）。
【００３５】
以上のように、本実施形態によれば、解析解が発散し易い差分大条件において、不足緩和係数βＩの値を小さくすることで、修正値を小さくできることから、解析解の発散を抑制することができ、安定的な解析を実行することができる。一方で、解析解が発散し難い差分小条件において、不足緩和係数βＩの値が大きくすることで、修正値を大きくできることから、解析解の収束を早めることができ、迅速な解析を実行することができる。このため、反復計算により解析解を安定して迅速に収束させることができる。
【００３６】
また、本実施形態によれば、不足緩和係数βＩを、取り扱いやすい関数を用いて表現することができる。つまり、係数Ａ，Ｋ，Ｃを適宜設定することで、関数を容易に変化させることができるため、例えば、入力パラメータの種類に応じた適切な関数を容易に設定することが可能となる。
【００３７】
また、本実施形態によれば、二流体の流体解析モデルに、ロジスティック関数となる不足緩和係数を適用することができるため、二流体の流体解析の解析解を安定して迅速に収束させることができる。
【００３８】
なお、本実施形態では、不足緩和係数βＩを流体解析モデルに適用して説明したが、入力パラメータを更新する反復計算を用いた解析モデルであれば、いずれに適用してもよい。このため、例えば、単相流の流体解析モデルである、いわゆるＳＩＭＰＬＥ法に、本実施形態の不足緩和係数βＩを適用してもよい。
【００３９】
また、本実施形態では、不足緩和係数βＩをロジスティック関数を用いて表現したが、ロジスティック関数に、特に限定されない。不足緩和係数βＩは、差分大条件において値が小さく、差分小条件において値が大きく、また、０より大きく１より小さい値であって、微分可能に連続する関数であれば、いずれの関数であってもよい。
【符号の説明】
【００４０】
１０  流体解析装置
１１  制御部
１２  記憶部
１３  入力部
１４  表示部
Ｐ  流体解析プログラム
【図１】
2019191905.tif 000007
【図２】
2019191905.tif 000008
【図３】
2019191905.tif 000009
【図４】
2019191905.tif 000010