JP5582878B2 - 数値解析装置及び要素生成プログラム - Google Patents

数値解析装置及び要素生成プログラム Download PDF

Info

Publication number
JP5582878B2
JP5582878B2 JP2010132114A JP2010132114A JP5582878B2 JP 5582878 B2 JP5582878 B2 JP 5582878B2 JP 2010132114 A JP2010132114 A JP 2010132114A JP 2010132114 A JP2010132114 A JP 2010132114A JP 5582878 B2 JP5582878 B2 JP 5582878B2
Authority
JP
Japan
Prior art keywords
volume
model
fluid
piping
junction
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.)
Active
Application number
JP2010132114A
Other languages
English (en)
Other versions
JP2011257999A (ja
Inventor
仁意 小野
利春 中林
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.)
Mitsubishi Heavy Industries Ltd
Original Assignee
Mitsubishi Heavy Industries Ltd
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 Mitsubishi Heavy Industries Ltd filed Critical Mitsubishi Heavy Industries Ltd
Priority to JP2010132114A priority Critical patent/JP5582878B2/ja
Priority to PCT/JP2011/063167 priority patent/WO2011155539A1/ja
Priority to US13/695,087 priority patent/US9372944B2/en
Publication of JP2011257999A publication Critical patent/JP2011257999A/ja
Application granted granted Critical
Publication of JP5582878B2 publication Critical patent/JP5582878B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00—Details relating to CAD techniques
    • G06F2111/10—Numerical modelling
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00—Details relating to the application field
    • G06F2113/08—Fluids
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00—Details relating to the application field
    • G06F2113/14—Pipes
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/10—Geometric CAD
    • G06F30/18—Network design, e.g. design based on topological or interconnect aspects of utility systems, piping, heating ventilation air conditioning [HVAC] or cabling
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06G—ANALOGUE COMPUTERS
    • G06G7/00—Devices in which the computing operation is performed by varying electric or magnetic quantities
    • G06G7/48—Analogue computers for specific processes, systems or devices, e.g. simulators
    • G06G7/57—Analogue computers for specific processes, systems or devices, e.g. simulators for fluid flow ; for distribution networks

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Geometry (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Complex Calculations (AREA)

Description

本発明は、数値解析装置及び要素生成プログラムに関するものである。
プラント等で用いられる流体機器及び配管で構成される管路網に対する種々の解析を行うためにＣＡＤ（Computer Aided Design）ツールやＣＡＥ（Computer Aided Engineering）ツール等が用いられている。
そして、構造解析や熱流動解析の分野では、ＣＡＤツールによって構築された３次元モデルから、ＣＡＥを行うための入力データを自動生成する種々のプリプロセッサが開発されている。そのため、ＣＡＤツールで構築されたモデルがあれば、ＣＡＥツールによる構造解析や熱流動解析が比較的容易に実現でき、ＣＡＤツールを用いた設計と、ＣＡＥツールによる解析による設計妥当性の評価、評価結果をもとにしたＣＡＤツールによる再設計のサイクルを、短時間で実現することが実現可能となっている。
また、特許文献１には、配管系統を、分岐合流点および機器に対応する単体要素としてのポイント・オブジェクトと、各ポイント・オブジェクト間に存在する複数の配管や機器を一括して表す複合要素としてのインターバル・オブジェクトとに分割し、各オブジェクトに、オブシェクト相互の接続関係を表す属性項目と、そのオブジェクトを構成する要素の属性を示す属性項目と設けるとともに、上記インターバル・オブジェクトの属性値の配列を物理的配管要素の配列と同順とし、このポイント・オブジェクトとインターバル・オブジェクトによって、各種モデルに共通の汎用モデルを構築する配管系統のモデリング方法が記載されている。
特開平０７−１０３４００号公報
ところで、従来、管路内の流体の温度や圧力、流速等の物理量の時間的な変化を解析（以下、「管路網非定常解析」という。）するための入力データを、ＣＡＤツールで構築した管路網モデルから自動生成することはできなかった。そのため、ＣＡＤツールで構築された管路網モデルから、人手を介して管路網非定常解析を行うための入力データを読み取る必要があった。また、特許文献１に記載の技術は、配管系統をポイントとインターバルに分割してモデル化する手法は記載されているが、ＣＡＤツール等を用いて算出された配管系統のデータのモデルへの反映方法については言及されていない。
本発明は、このような事情に鑑みてなされたものであって、３次元モデルとして構築された管路網モデルから管路網における流体の非定常解析を行うために必要な情報を読み出すことができる数値解析装置及び要素生成プログラムを提供することを目的とする。
上記課題を解決するために、本発明の数値解析装置は以下の手段を採用する。
すなわち、本発明に係る数値解析装置は、流体機器及び配管を含む管路網の３次元モデルとして構築された管路網モデルに基づいて、ボリューム・ジャンクション法により流体の物理量の時間的な変化の解析である非定常解析を行う数値解析装置であって、予め構築された前記管路網モデルにおいて設定された始点と設定された終点との間に位置する流体機器モデル及び配管モデルのうち、前記非定常解析に用いる該流体機器モデル及び該配管モデルを選択する選択手段と、前記選択手段によって選択された前記流体機器モデル及び前記配管モデルを、該流体機器モデル及び該配管モデルの容積を示す要素である複数のボリューム要素、並びに前記ボリューム要素と他の前記ボリューム要素との境界を示すと共に流体の運動エネルギーの増減を示す要素である複数のジャンクション要素に分割する分割手段と、前記分割手段によって分割された前記ボリューム要素毎の容積及び前記ジャンクション要素毎に対応した圧力損失の大きさを示す係数を、前記流体機器モデルの形状、前記配管モデルの形状、並びに前記管路網モデルを流れる流体の物理量に基づいて導出し、導出した該容積を前記ボリューム要素に対応付け、導出した該係数を前記ジャンクション要素に対応付ける導出手段と、流体の流量と各前記ボリューム要素の容積とから該流体の各前記ボリューム要素毎の通過時間を算出する算出手段と、前記算出手段によって算出された通過時間が各前記ボリューム要素において所定の差以内となるように、前記ボリューム要素及び前記ジャンクション要素を再び分割する再分割手段と、を備える。
本発明によれば、選択手段によって、流体機器及び配管を含む管路網の３次元モデルとして予め構築された管路網モデルにおいて設定された始点と設定された終点との間に位置する流体機器モデル及び配管モデルのうち、非定常解析に用いる該流体機器モデル及び該配管モデルが選択される。
そして、分割手段によって、選択手段で選択された流体機器モデル及び配管モデルが、該流体機器モデル及び該配管モデルの容積を示す要素である複数のボリューム要素、並びにボリューム要素と他のボリューム要素との境界を示すと共に流体の運動エネルギーの増減を示す要素である複数のジャンクション要素に分割される。さらに、導出手段によって、分割手段で分割されたボリューム要素毎の容積、及びジャンクション要素毎に対応した圧力損失の大きさを示す係数が流体機器モデルの形状及び配管モデルの形状に基づいて導出され、導出された該容積がボリューム要素に対応付けられ、導出された該係数がジャンクション要素に対応付けられる。
なお、上記ボリューム要素及びジャンクション要素は、ボリューム・ジャンクション法により流体の物理量の時間的な変化の解析である非定常解析に用いられる。
このように、本発明は、３次元モデルである管路網モデルを構成する流体機器モデル及び配管モデルが、ボリューム要素及びジャンクション要素に分割され、流体機器モデルの形状、配管モデルの形状、及び流体の物理量に基づいて、ボリューム要素毎の容積及びジャンクション要素毎の圧損係数が導出されるので、３次元モデルとして構築された管路網モデルから管路網における流体の非定常解析を行うために必要な情報を読み出すことができる。
また、本発明は、流体の流量と各前記ボリューム要素の容積とから該流体の各前記ボリューム要素毎の通過時間を算出する算出手段と、前記算出手段によって算出された通過時間が各前記ボリューム要素において所定の差以内となるように、前記ボリューム要素及び前記ジャンクション要素を再び分割する再分割手段と、をさらに備えてもよい。
本発明によれば、算出手段によって流体の流量と各ボリューム要素の容積とから該流体の各ボリューム要素毎の通過時間が算出され、再分割手段によって、算出手段で算出された通過時間が各ボリューム要素において所定の差以内となるように、ボリューム要素及びジャンクション要素が再び分割される。これにより、本発明は、各ボリューム要素を通過する流体の通過時間を同等にすることができ、ボリューム・ジャンクション法による解析結果をより安定にすることができる。
また、本発明は、前記ジャンクション要素を介して隣接する前記ボリューム要素の位置情報から、該ボリューム要素間の高度差のみを抽出し登録する登録手段をさらに備えてもよい。
本発明によれば、登録手段によって、ジャンクション要素を介して隣接するボリューム要素間の高度差が登録されるので、該高度差を用いることで、ボリューム・ジャンクション法による非定常解析において、３次元の機器配置情報から高度差のみの１次元情報にモデルパラメータの数を少なくすることができ、計算負荷を抑制することができる。
また、本発明は、前記導出手段が、前記係数の大きさが同じ場合は、前記係数を示す情報を統合させてもよい。
本発明によれば、圧力損失の大きさを示す係数が同じ場合は、該係数を示す情報を統合させるので、ボリューム・ジャンクション法で非定常解析を行うために用いる情報の量を少なくすることができる。
また、上記課題を解決するために、本発明の要素生成プログラムは以下の手段を採用する。
すなわち、本発明に係る要素生成プログラムは、流体機器及び配管を含む管路網の３次元モデルとして構築された管路網モデルに基づいて、ボリューム・ジャンクション法により流体の物理量の時間的な変化の解析である非定常解析を行う数値解析装置の要素生成プログラムであって、コンピュータを、予め構築された前記管路網モデルにおいて設定された始点と設定された終点との間に位置する流体機器モデル及び配管モデルのうち、前記非定常解析に用いる該流体機器モデル及び該配管モデルを選択する選択手段と、前記選択手段によって選択された前記流体機器モデル及び前記配管モデルを、該流体機器モデル及び該配管モデルの容積を示す要素である複数のボリューム要素、並びに前記ボリューム要素と他の前記ボリューム要素との境界を示すと共に流体の運動エネルギーの増減を示す要素である複数のジャンクション要素に分割する分割手段と、前記分割手段によって分割された前記ボリューム要素毎の容積及び前記ジャンクション要素毎に対応した圧力損失の大きさを示す係数を、前記流体機器モデルの形状、前記配管モデルの形状、並びに前記管路網モデルを流れる流体の物理量に基づいて導出し、導出した該容積を前記ボリューム要素に対応付け、導出した該係数を前記ジャンクション要素に対応付ける導出手段と、流体の流量と各前記ボリューム要素の容積とから該流体の各前記ボリューム要素毎の通過時間を算出する算出手段と、前記算出手段によって算出された通過時間が各前記ボリューム要素において所定の差以内となるように、前記ボリューム要素及び前記ジャンクション要素を再び分割する再分割手段と、して機能させるための要素生成プログラムである。
本発明によれば、３次元モデルである管路網モデルを構成する流体機器モデル及び配管モデルが、ボリューム要素及びジャンクション要素に分割され、流体機器モデルの形状、配管モデルの形状、及び流体の物理量に基づいて、ボリューム要素毎の容積及びジャンクション要素毎の圧損係数が導出されるので、３次元モデルとして構築された管路網モデルから管路網における流体の非定常解析を行うために必要な情報を読み出すことができる。
３次元モデルとして構築された管路網モデルから管路網における流体の非定常解析を行うために必要な情報を読み出す。
本発明の第１実施形態に係る数値解析装置の電気的構成を示すブロック図である。 本発明の第１実施形態に係る３次元モデルとして構築された管路網モデルの一例を示す模式図である。 本発明の第１実施形態に係る非定常解析に用いる管路網ネットワークの一例を示す模式図である。 本発明の第１実施形態に係る要素生成プログラムの処理の流れを示すフローチャートである。 本発明の第２実施形態に係るエレベーション差を登録する処理の流れを示すフローチャートである。
以下に、本発明に係る数値解析装置及び要素生成プログラムの一実施形態について、図面を参照して説明する。
〔第１実施形態〕
以下、本発明の第１実施形態について説明する。図１に、本第１実施形態に係る数値解析装置１０の電気的構成を示す。
数値解析装置１０は、数値解析装置１０全体の動作を司るＣＰＵ（Central Processing Unit）１２、各種プログラムや各種パラメータ等が予め記憶されたＲＯＭ（Read Only Memory）１４、ＣＰＵ１２による各種プログラムの実行時のワークエリア等として用いられるＲＡＭ（Random Access Memory）１６、ＣＡＥによる数値解析プログラム及び詳細を後述する要素生成プログラム等の各種プログラム及び各種情報を記憶する記憶手段としてのＨＤＤ（Hard Disk Drive）１８を備えている。
さらに、数値解析装置１０は、キーボード及びマウス等から構成され、各種操作の入力を受け付ける操作入力部２０、ＣＡＤツールによって構築された管路網モデルや数値解析の結果等の各種情報を表示する、例えばＬＣＤ（Liquid Crystal Display）等の画像表示部２２、並びに他の情報処理装置及びプリンタ等の外部装置と接続され、該外部装置への各種情報の送受信を行う外部インタフェース２４を備えている。
これらＣＰＵ１２、ＲＯＭ１４、ＲＡＭ１６、ＨＤＤ１８、操作入力部２０、画像表示部２２、及び外部インタフェース２４は、システムバス３０を介して相互に電気的に接続されている。従って、ＣＰＵ１２は、ＲＯＭ１４、ＲＡＭ１６、及びＨＤＤ１８へのアクセス、操作入力部２０に対する操作状態の把握、画像表示部２２に対する各種の画像の表示、並びに外部インタフェース２４を介した上記外部装置との各種情報の送受信等を各々行なうことができる。
なお、ＨＤＤ１８には、流体機器及び配管で構成される管路網をＣＡＤツールによって３次元モデルとして構築された管路網モデルを示す管路網モデルデータが予め記憶されている。また、ＣＡＤツールは、数値解析装置１０にインストールされ、数値解析装置１０によって管路網モデルが構築されていてもよいし、例えば、外部インタフェース２４を介して接続されている他の情報処理装置を用いて管路網モデルが構築され、該管路網モデルがＨＤＤ１８に記憶されてもよい。
図２にＣＡＤツールで作成された管路網モデルの一例（配管アイソメ図）を示す。
同図に示す管路網モデルは、流体機器としてタンク４０Ａ，４０Ｂ、弁４２Ａ〜４２Ｆ、ポンプ４４を含んでおり、各流体機器の間は、配管５０で結合されている。また、配管５０においても、口径の異なる配管５０同士はフランジを介して結合されている。
なお、管路網モデルに含まれる流体機器モデル及び配管モデルは、各々機器番号が付されており、各流体機器モデル及び配管モデルは、該機器番号によってその種類（例えば、配管であればその種別（直管、拡大、縮小、ベント、分岐、合流））が識別される。また、流体機器モデル及び配管モデルは、その形状（例えば、配管の形状データとして、等価粗さ、直管内径、長さ、入口内径、出口内径、拡大角、ベント内径、曲率半径、及び曲げ角度等）が入力されている。
さらに、管路網モデルには、ＣＡＥツールによる非定常解析において流体の物理量（温度、圧力、流速、及び流量等）を検出する個所を示す検出端が、予め設置されている。
一方、ＣＡＥツールによる管路網を流れる流体に対する非定常解析は、質量保存式、運動量保存式、及びエネルギー保存式の３種類の保存式（連続の式）からなる基礎方程式を用いて、流体の物理量を数値シミュレーションで解析する。そのため、基礎方程式をコンピュータで計算可能な式に変換する必要がある。すなわち、流体の流れ方向に連続的に変化する温度、圧力、流速、及び流量といった物理量を、空間で分割して離散化する必要がある。
このため、管路網を構成する流路をボリューム要素で分割し、ボリューム要素とボリューム要素の境界をジャンクション要素で接続するボリューム・ジャンクション法により、上記基礎方程式を空間的に離散化する。
なお、ボリューム要素とは、流体が流れる流路の容積を示す要素であり、具体的には、配管、タンク、熱交換器、ヘッダー配管等をボリューム要素である。また、ジャンクション要素とは、流体の運動エネルギーの増減（流体の圧力を減少させる抵抗や流体に圧力を与えるもの）を示す要素であり、具体的には、弁、オリフィス、及びフランジ等の流路面積が急に拡大又は急に縮小する部位、並びにポンプ等である。
図３に、図２に示した管路網モデルに基づいて生成したボリューム要素及びジャンクション要素を模式化した管路網ネットワーク図を示す。
図３に示す管路網ネットワーク図では、各配管５０、タンク４０Ａ，４０Ｂをボリューム要素としている。一方、弁４２Ａ〜４２Ｆ、ポンプ４４、フランジ５２をジャンクション要素としている。そして、同図に示すように、ボリューム要素とボリューム要素との境界には、ジャンクション要素が配置されている。なお、各ボリューム要素は、フランジ５２を介して隣り合う配管５０の容積の平均値が対応する容積となる。
そして、本第１実施形態に係る数値解析装置１０は、ボリューム・ジャンクション法による非定常解析を行うために、管路網モデルに基づいて、該管路網モデルに対応したボリューム要素及びジャンクション要素を生成する要素生成処理を行う。すなわち、本第１実施形態に係る数値解析装置１０は、上記ボリューム要素及びジャンクション要素を生成する要素生成装置としても作用する。
次に、本第１実施形態に係る数値解析装置１０の作用を説明する。
図４は、要素生成処理を行う場合に、ＣＰＵ１２によって実行される要素生成プログラムの処理の流れを示すフローチャートであり、該要素生成プログラムはＨＤＤ１８の所定領域に予め記憶されている。
まず、ステップ１００では、ＨＤＤ１８に記憶されている管路網モデルを読み出し、該管路網モデルにおいて、非定常解析の対象となる始点及び終点を設定する。なお、該始点及び終点は、操作入力部２０を介してユーザによって設定される。
次のステップ１０２では、ステップ１００において設定された始点と終点との間に配置されている流体機器モデル及び配管モデルのうち、非定常解析に用いる流体機器モデル及び配管モデルを選択する。
具体的には、本第１実施形態では、管路網モデルにおいて、各流体機器モデル毎に付されている機器番号により、各々を区別しているため、該機器番号に基づいて非定常解析に用いる流体機器モデル及び配管モデルを選択する。また、流体機器モデル同士を結合させる配管モデルのうち、上述した検出端が設けられている配管モデル、及び予めユーザにより設定された口径以上の配管モデルが非定常解析に用いるものとして選択される。さらに、流体機器モデル同士を結合させる配管モデルでなくても、管路網モデル系外への流体の流れを示す配管であるベント管及びドレン管を示すモデル等を、非定常解析に用いるものとしてユーザにより選択できるとしてもよい。
次のステップ１０４では、ステップ１０２において選択された流体機器モデル及び配管モデルをボリューム要素及びジャンクション要素に分割する。
本第１実施形態では、上述したように、配管、タンク、熱交換器、ヘッダー配管等の配管の分岐部分は、ボリューム要素として分割され、ポンプ、並びに流路面積が拡大又は縮小する部分である流体機器（弁、オリフィス、及びフランジ等）は、ジャンクション要素として分割される。また、温度、圧力、及び濃度等を検出するための検出端が設置されているモデルは、ボリューム要素として分割され、流速及び流量等の流体の流れの大きさを示す物理量を検出するための検出端が設置されているモデルは、ジャンクション要素として分割される。なお、上述したように、ジャンクション要素は、ボリューム要素とボリューム要素との間に配置される。
次のステップ１０６では、管路網モデルのボリューム要素として分割された流体機器モデル及び配管モデルから、流体が流れる管の口径及び管の長さを抽出することによって、ボリューム要素毎の容積を導出する。また、流体機器モデルの形状、配管モデルの形状、及び管路網モデルを流れる流体の物理量に基づいて圧力損失の大きさを示す係数(以下、「圧損係数」という。)を導出する。そして、導出した容積をボリューム要素毎に対応付け、導出した圧損係数をジャンクション要素毎に対応付ける。
なお、上記管の口径及び管の長さとは、配管モデルの口径及び長さ、並びに流体機器モデル内の配管の口径及び長さをいう。また、流体の物理量とは、例えば、始点に入力される流体の成分（濃度）、基準圧力、基準温度、及び基準流量である。
なお、本第１実施形態では、連続している同一口径の管は、予め定められた容積以下となるように分割する。そして、該分割した容積をボリューム要素として対応付ける。
また、圧損係数を導出するためには、圧損係数を導出するための他のプログラムを用いて導出してもよい。この場合、流体機器モデルの形状及び配管モデルの形状を示す情報を該プログラムに入力する。そして、該プログラムによって導出された圧損係数を読み込み、各ジャンクション要素毎に対応付ける。
なお、導出した圧損係数が同じ場合は、該圧損係数を示す情報を統合させる。具体的には、例えば、同一の圧損係数を示す情報を所定のメモリ領域に記憶させ、ジャンクション要素には、該所定のメモリ領域のアドレスを示すアドレス情報を対応付ける。そして、非定常解析を行う場合に、上記アドレス情報に基づいて、記憶させた圧損係数を読み出す。これにより、ボリューム・ジャンクション法で非定常解析を行うために用いる情報の量を少なくすることができる。
次のステップ１０８では、流体の流量と各ボリューム要素の容積とから該流体の各ボリューム要素毎の通過時間を算出する。なお、ここでいう流体の流量とは、たとえば、管路網モデルが、定格運転している場合における体積流量である。
次のステップ１１０では、ステップ１０８で算出した通過時間が各ボリューム要素において所定の差以内となっているか否かを判定し、肯定判定の場合は、本プログラムを終了する一方、否定判定の場合は、ステップ１１２へ移行する。具体的には、本ステップ１００のおける処理は、ステップ１０８で算出した通過時間が最も短いボリューム要素を特定し、最も短い通過時間よりも所定倍（一例として２倍）の通過時間のボリューム要素がある場合を、上記所定の差を越えたボリューム要素があるとし、ステップ１１０で否定判定とする。
ステップ１１２では、上記所定の差を越えたボリューム要素をさらに予め定められた数で再分割し、ステップ１０６へ戻る。これにより、各ボリューム要素を通過する流体の通過時間を同等にすることができ、ボリューム・ジャンクション法による解析結果をより安定にすることができる。なお、分割したボリューム要素間には新たにジャンクション要素が配置される。
そして、本第１実施形態にかかる数値解析装置１０は、本プログラムが終了した後に、本プログラムによってボリューム要素及びジャンクション要素が配置された管路網ネットワークに基づいて、ボリューム・ジャンクション法を用いた数値解析プログラムによって非定常解析を行う。
以上説明したように、本第１実施形態に係る数値解析装置１０は、３次元モデルである管路網モデルを構成する流体機器モデル及び配管モデルがボリューム要素及びジャンクション要素に分割され、流体機器モデルの形状、配管モデルの形状、及び流体の物理量に基づいて、ボリューム要素毎の容積及びジャンクション要素毎の圧損係数が導出されるので、３次元モデルとして構築された管路網モデルから管路網における流体の非定常解析を行うために必要な情報を読み出すことができる。
〔第２実施形態〕
以下、本発明の第２実施形態について説明する。
なお、本第２実施形態に係る数値解析装置１０の電気的構成は、図１と同様であるので説明を省略する。
次に、本第２実施形態に係る数値解析装置１０の作用を説明する。
なお、第２実施形態に係る要素生成プログラムの処理の流れは、図４に示す処理と同様であるが、図４のステップ１０６において、ボリューム要素の位置を示す情報から、ジャンクション要素を介して隣接するボリューム要素の高度差（以下、「エレベーション差」という。）を算出する。
図５に、図４のステップ１０６で行われるエレベーション差を導出する処理の流れを示す。
まずステップ２００では、各ボリューム要素の中心点の座標（ｘ，ｙ，ｚ）を管路網モデルから抽出する。
次のステップ２０２では、ステップ２００で抽出した中心点の座標（ｘ，ｙ，ｚ）のうち、高さ方向の情報（ｚ）の値を抽出する。
次のステップ２０４では、ジャンクション要素を介して隣接するボリューム要素のエレベーション差を、ステップ２０２で抽出した該ボリューム要素の高さ方向の情報（ｚ）から算出する。
次のステップ２０６では、ステップ２０４で算出したエレベーション差をボリューム要素又はジャンクション要素に登録する。
以上説明したように、本第２実施形態に係る数値解析装置１０では、ジャンクション要素を介して隣接するボリューム要素の位置情報として、該ボリューム要素間のエレベーション差のみを登録するので、登録したエレベーション差を用いることで、ボリューム・ジャンクション法による非定常解析で用いるパラメータの数を少なくすることができ、計算負荷を抑制することができる。
以上、本発明を、上記各実施形態を用いて説明したが、本発明の技術的範囲は上記実施形態に記載の範囲には限定されない。発明の要旨を逸脱しない範囲で上記各実施形態に多様な変更または改良を加えることができ、該変更または改良を加えた形態も本発明の技術的範囲に含まれる。
例えば、上記各実施形態では、ＨＤＤ１８に要素生成プログラムを記憶させていたが、本発明は、これに限定されるものではなく、要素生成プログラムを、磁気ディスク、ＣＤ（Compact Disc）、ＤＶＤ（Digital Versatile Disc）等の光ディスク、メモＩＣ（Integrated Circuit）カード、及びメモリカード等の可搬型記憶媒体に記憶させ、該可搬型記憶媒体に記憶されている要素生成プログラムを数値解析装置１０が読み取り、実行する形態としてもよい。
１０ 数値解析装置
１２ ＣＰＵ
１８ ＨＤＤ
２０ 操作入力部

Claims (4)

  1. 流体機器及び配管を含む管路網の３次元モデルとして構築された管路網モデルに基づいて、ボリューム・ジャンクション法により流体の物理量の時間的な変化の解析である非定常解析を行う数値解析装置であって、
    予め構築された前記管路網モデルにおいて設定された始点と設定された終点との間に位置する流体機器モデル及び配管モデルのうち、前記非定常解析に用いる該流体機器モデル及び該配管モデルを選択する選択手段と、
    前記選択手段によって選択された前記流体機器モデル及び前記配管モデルを、該流体機器モデル及び該配管モデルの容積を示す要素である複数のボリューム要素、並びに前記ボリューム要素と他の前記ボリューム要素との境界を示すと共に流体の運動エネルギーの増減を示す要素である複数のジャンクション要素に分割する分割手段と、
    前記分割手段によって分割された前記ボリューム要素毎の容積及び前記ジャンクション要素毎に対応した圧力損失の大きさを示す係数を、前記流体機器モデルの形状、前記配管モデルの形状、並びに前記管路網モデルを流れる流体の物理量に基づいて導出し、導出した該容積を前記ボリューム要素に対応付け、導出した該係数を前記ジャンクション要素に対応付ける導出手段と、
    流体の流量と各前記ボリューム要素の容積とから該流体の各前記ボリューム要素毎の通過時間を算出する算出手段と、
    前記算出手段によって算出された通過時間が各前記ボリューム要素において所定の差以内となるように、前記ボリューム要素及び前記ジャンクション要素を再び分割する再分割手段と、
    を備えた数値解析装置。
  2. 前記ジャンクション要素を介して隣接する前記ボリューム要素の位置情報から、該ボリューム要素間の高度差のみを抽出し登録する登録手段をさらに備えた請求項１記載の数値解析装置。
  3. 前記導出手段は、前記係数の大きさが同じ場合は、前記係数を示す情報を統合させる請求項１または請求項２記載の数値解析装置。
  4. 流体機器及び配管を含む管路網の３次元モデルとして構築された管路網モデルに基づいて、ボリューム・ジャンクション法により流体の物理量の時間的な変化の解析である非定常解析を行う数値解析装置の要素生成プログラムであって、
    コンピュータを、
    予め構築された前記管路網モデルにおいて設定された始点と設定された終点との間に位置する流体機器モデル及び配管モデルのうち、前記非定常解析に用いる該流体機器モデル及び該配管モデルを選択する選択手段と、
    前記選択手段によって選択された前記流体機器モデル及び前記配管モデルを、該流体機器モデル及び該配管モデルの容積を示す要素である複数のボリューム要素、並びに前記ボリューム要素と他の前記ボリューム要素との境界を示すと共に流体の運動エネルギーの増減を示す要素である複数のジャンクション要素に分割する分割手段と、
    前記分割手段によって分割された前記ボリューム要素毎の容積及び前記ジャンクション要素毎に対応した圧力損失の大きさを示す係数を、前記流体機器モデルの形状、前記配管モデルの形状、並びに前記管路網モデルを流れる流体の物理量に基づいて導出し、導出した該容積を前記ボリューム要素に対応付け、導出した該係数を前記ジャンクション要素に対応付ける導出手段と、
    流体の流量と各前記ボリューム要素の容積とから該流体の各前記ボリューム要素毎の通過時間を算出する算出手段と、
    前記算出手段によって算出された通過時間が各前記ボリューム要素において所定の差以内となるように、前記ボリューム要素及び前記ジャンクション要素を再び分割する再分割手段と、
    して機能させるための要素生成プログラム。
JP2010132114A 2010-06-09 2010-06-09 数値解析装置及び要素生成プログラム Active JP5582878B2 (ja)

Priority Applications (3)

Application Number Priority Date Filing Date Title
JP2010132114A JP5582878B2 (ja) 2010-06-09 2010-06-09 数値解析装置及び要素生成プログラム
PCT/JP2011/063167 WO2011155539A1 (ja) 2010-06-09 2011-06-08 数値解析装置、要素生成プログラム、及び数値解析方法
US13/695,087 US9372944B2 (en) 2010-06-09 2011-06-08 Numerical analysis device, element generation program, and numerical analysis method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2010132114A JP5582878B2 (ja) 2010-06-09 2010-06-09 数値解析装置及び要素生成プログラム

Publications (2)

Publication Number Publication Date
JP2011257999A JP2011257999A (ja) 2011-12-22
JP5582878B2 true JP5582878B2 (ja) 2014-09-03

Family

ID=45098143

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2010132114A Active JP5582878B2 (ja) 2010-06-09 2010-06-09 数値解析装置及び要素生成プログラム

Country Status (3)

Country Link
US (1) US9372944B2 (ja)
JP (1) JP5582878B2 (ja)
WO (1) WO2011155539A1 (ja)

Families Citing this family (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
IN2014DN05788A (ja) * 2012-01-13 2015-05-15 Process Systems Entpr Ltd
JP5900951B2 (ja) * 2012-01-17 2016-04-06 国立大学法人広島大学 モデル化方法、解析方法及びプログラム
CN103955186B (zh) * 2014-04-22 2016-08-24 中国石油大学(北京) 天然气管网管流状态参数确定方法及装置
JP2020144409A (ja) * 2019-03-04 2020-09-10 株式会社日立製作所 プラント診断装置
JP7280123B2 (ja) * 2019-06-26 2023-05-23 株式会社日立製作所 ３ｄモデル作成支援システムおよび３ｄモデル作成支援方法
CN113434993B (zh) * 2020-03-23 2023-01-31 株洲中车时代电气股份有限公司 一种管网流量计算系统及方法
JP6908758B1 (ja) * 2020-05-08 2021-07-28 清 菅井 ガス流量解析方法
CN112257014B (zh) * 2020-10-22 2024-01-05 国药集团重庆医药设计院有限公司 一种低压饱和蒸汽管道压力的计算方法及系统
JP7257472B2 (ja) * 2020-10-29 2023-04-13 みずほリサーチ＆テクノロジーズ株式会社 設計支援システム、設計支援方法及び設計支援プログラム

Family Cites Families (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH07103400A (ja) 1993-10-07 1995-04-18 Toshiba Corp 配管系統のモデリング方法
JP3556384B2 (ja) * 1996-03-28 2004-08-18 三菱電機株式会社 水路系模擬装置
JPH11305646A (ja) * 1998-04-23 1999-11-05 Mitsubishi Electric Corp 初期値生成方法および装置
JP3236835B2 (ja) * 1999-06-17 2001-12-10 核燃料サイクル開発機構 一次元プラントモデル非定常流動評価方法、評価システム、及び記憶媒体
EP2472391B1 (en) * 2002-05-20 2016-10-26 Tyco Fire Products LP System for analyzing models of dry pipe systems
JP2005050027A (ja) * 2003-07-31 2005-02-24 Toshiba Corp 流動解析システムおよび流動解析方法
JP4585323B2 (ja) * 2005-01-25 2010-11-24 株式会社東芝 プラント配管エンジニアリング計算用リスト出力装置およびプラント配管エンジニアリング計算装置
EP1934709A4 (en) * 2005-10-03 2013-03-13 Central Sprinkler Company SYSTEM AND METHOD FOR EVALUATING A FLUID FLOW IN A PIPE SYSTEM

Also Published As

Publication number Publication date
JP2011257999A (ja) 2011-12-22
US20130046527A1 (en) 2013-02-21
US9372944B2 (en) 2016-06-21
WO2011155539A1 (ja) 2011-12-15

Similar Documents

Publication Publication Date Title
JP5582878B2 (ja) 数値解析装置及び要素生成プログラム
JP3522408B2 (ja) 数値流体解析結果の誤差見積方法、数値流体解析結果の誤差見積装置、数値流体解析方法、及び数値流体解析装置
SaemI et al. Computation of two-and three-dimensional water hammer flows
Stel et al. Investigation of the motion of bubbles in a centrifugal pump impeller
Gorman et al. Evaluation of the efficacy of turbulence models for swirling flows and the effect of turbulence intensity on heat transfer
JP6429682B2 (ja) ネットワークシミュレーション装置、ネットワークシミュレーション方法及びネットワークシミュレーションプログラム
JP2017211987A (ja) 配管の経路生成方法、複数の配管の経路生成方法、配管の経路生成装置、複数の配管の経路生成装置及びプログラム
Meziou et al. Low-dimensional modeling of transient two-phase flow in pipelines
Sondermann et al. Numerical simulation of non-isothermal two-phase flow in pipelines using a two-fluid model
Gooding et al. The impact of various modeling decisions on flow field predictions in a centrifugal compressor
JP2006258391A (ja) 浮力を伴う乱流の流体的及び熱的諸特性の推定方法及び推定プログラム
Sharp et al. Optimizing the ASME Venturi recovery cone angle to minimize head loss
Surendra et al. Investigation of swirl flows applied to the oil and gas industry
WO2024051724A1 (zh) 流致振动试验装置、方法、计算机设备、存储介质和产品
Aly Computational fluid dynamics and optimization of flow and heat transfer in coiled tube-in-tube heat exchangers under turbulent flow conditions
JP6089755B2 (ja) 熱対策機器導入省エネ試算評価システム及び評価方法
Zanino et al. CFD modeling of ITER cable‐in‐conduit superconductors. Part I: Friction in the central channel
Coussirat et al. Study of available turbulence and cavitation models to reproduce flow patterns in confined flows
Himr Investigation and numerical simulation of a water hammer with column separation
Figueiredo et al. A straightforward strategy for leak localization in two-phase gas pipelines
Karbon et al. Large-eddy simulation of the flow in Z-Shape duct
US11113436B1 (en) Method for simulation of flow in fluid flow network having one-dimensional and multi-dimensional flow components
Schab et al. Parameter study of geometrically induced flow maldistribution in shell and tube heat exchangers
Ishii et al. Numerical Approach to Study Nonuniform Gas–Liquid Distribution in the Refrigerant Distributor in an Air Conditioner
Guyot et al. A Vortex modeling with 3D CFD

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20130529

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20140401

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20140602

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20140624

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20140715

R151 Written notification of patent or utility model registration

Ref document number: 5582878

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R151

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250