ガス導管内の物性変化推定プログラム、物性変化推定装置、及び物性変化推定方法
【課題】ガスの輸送導管を成分の異なる複数種類のガス輸送に使用する場合における導管内の物性変化推定プログラムであって、ガスの混合比が変化した場合の導管下流側における濃度変化を、低コストでかつ適確に把握することのできる物性変化推定プログラム等を提供する。
【解決手段】ガス導管内を流れる互いに種類の異なるガスの混合比を変更させた場合における、ガスの物性変化挙動を推定する処理をコンピュータに実行させる物性変化推定プログラムが、推定開始後の各時刻において、推定処理の対象範囲を所定長さで分割した第一分割領域毎に、混合比に基づいて、ガスの流速を求めるステップと、前記第一分割領域を更に分割した第二分割領域毎に、流速に基づいて混合比を求めるステップとを、前記コンピュータに実行させる。
【解決手段】ガス導管内を流れる互いに種類の異なるガスの混合比を変更させた場合における、ガスの物性変化挙動を推定する処理をコンピュータに実行させる物性変化推定プログラムが、推定開始後の各時刻において、推定処理の対象範囲を所定長さで分割した第一分割領域毎に、混合比に基づいて、ガスの流速を求めるステップと、前記第一分割領域を更に分割した第二分割領域毎に、流速に基づいて混合比を求めるステップとを、前記コンピュータに実行させる。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、ガスの輸送導管を成分の異なる複数種類のガス輸送に使用する場合における、導管内の物性変化推定プログラム等に関し、特に、ガスの混合比が変化した場合の導管下流側における濃度変化を、低コストでかつ適確に把握することのできる物性変化推定プログラム等に関する。
【背景技術】
【０００２】
都市ガスは、地域に張り巡らせたガス導管網を経由して、ガス製造工場等のガス供給元からガスを使用するガス消費設備等へ供給される。また、供給する都市ガスは、複数の地域から産出され組成のバラツキがある天然ガスに、ＬＰＧ、ブタン、プロパンなどを混合して熱量が一定になるように調整がなされる。しかし、ガスの輸送導管の有効活用や、ガス輸送導管を持たない事業者への公平な参入機会の創出などの要請から、複数の地域から産出した天然ガスをそのままガス導管網に導入するなど、熱量がそれぞれ異なる複数種類のガスがガス導管網に導入されることもある。
【０００３】
また、今後は、益々エネルギーの多様化の要請が高まり、ガス導管網をこれまで以上に性状の異なる異種ガスの輸送に使用することが考えられる。さらに、ガス供給の新規参入者等のガス輸送を代行する託送に、ガス導管網が利用される頻度が多くなることが予想され、これによってもガス導管網において、性状のかなり異なるガスの輸送を行なう機会が増えると考えられる。
【０００４】
このような状況では、ガス導管内を流れるガスの性状が変化する機会が増え、そのような場合には、供給されるガスを消費する設備等においては、ガスの性状が大きく変化する場合があり、その変化挙動を適確に把握しなければ設備の運転に支障をきたしたり、危険が生じてしまう場合もある。
【０００５】
そこで、下記特許文献１では、輸送するガスの置換、混合比の変更による導管内の物性変化挙動を所定の拡散方程式を用いて求めることについて提案がなされている。
【０００６】
また、下記特許文献２には、上記拡散方程式における拡散定数を、導管内のガス流れが属する各領域毎に定められた所定の関係式によって決定される技術が開示されている。
【０００７】
また、下記非特許文献１では、配管破断事故後に生じる自然循環の特性を調べるための、２成分気体の拡散と自然循環に関する研究について示され、下記非特許文献２には、ガスパイプラインの非定常流解析シミュレーションについて示されている。
【先行技術文献】
【特許文献】
【０００８】
【特許文献１】特開２００６−３３６７９５号公報
【特許文献２】特開２００８−２４０８６４号公報
【非特許文献】
【０００９】
【非特許文献１】武田哲明、菱田誠、「２成分気体の拡散と自然循環に関する研究」、日本機械学会論文集、１９８９年８月、５５巻、５１６号、ｐ．２７９５−２７９９
【非特許文献２】山根総一郎、佐藤律夫、「ガスパイプライン異種ガス対応非定常流解析シミュレータの開発」、ＪＦＥ技報、２００６年１１月、Ｎｏ．１４、ｐ．６０−６４
【発明の概要】
【発明が解決しようとする課題】
【００１０】
しかしながら、上記特許文献１、２に記載の方法は、導管内の流速、導管の内径等の仕様が一定である場合、推定開始時の異種ガスの界面が単純な形状をしている場合に適用可能であり、適用範囲が限られるという課題があった。
【００１１】
また、上記非特許文献１に記載の方法は、対象範囲が層流であり導管内の流れが乱流である場合には対応できず、また、管の長さが短い場合を対象としているので、対象が上述したガスの輸送導管（網）のように長い場合には、解析のための手法が計算負荷の面で向いていない。
【００１２】
また、上記非特許文献２に記載の方法では、気体の移流のみによる効果が反映され、気体の拡散が考慮されていないので、推定結果の精度に課題がある。
【００１３】
そこで、本発明の目的は、ガスの輸送導管を成分の異なる複数種類のガス輸送に使用する場合における導管内の物性変化推定プログラムであって、ガスの混合比が変化した場合の導管下流側における濃度変化を、低コストでかつ適確に把握することのできると共に、様々な条件下における実際のケースに使用可能な物性変化推定プログラム、等を提供することである。
【課題を解決するための手段】
【００１４】
上記の目的を達成するために、本発明の一つの側面は、ガス導管内を流れる互いに種類の異なるガスの混合比を変更させた場合における、前記ガス導管内を流れるガスの物性変化挙動を推定する処理をコンピュータに実行させる物性変化推定プログラムが、前記推定開始後の各時刻において、前記推定処理の対象範囲を所定長さで分割した第一分割領域毎に、前記混合比に基づいて、前記ガスの流速を求める流速決定ステップと、前記第一分割領域を更に分割した第二分割領域毎に、前記流速に基づいて前記混合比を求める（すなわち、界面が存在する領域については界面形状を求める）混合比決定ステップとを、前記コンピュータに実行させる、ことである。
【００１５】
更に、上記の発明において、その好ましい態様は、前記混合比決定ステップは、前記混合比が変化する部分である界面が存在する領域について実行される、ことを特徴とする。
【００１６】
また、上記の発明において、好ましい態様は、更に、少なくとも、前記第二分割領域毎の前記ガス導管の仕様と、前記推定開始時における前記ガス導管の各位置における前記混合比と、前記ガス導管の最上流における前記ガスの流速及び圧力を取得するステップを、前記コンピュータに実行させ、前記流速決定ステップと前記混合比決定ステップにおいて、当該取得された情報が用いられる、ことを特徴とする。
【００１７】
更に、上記の発明において、一つの態様は、前記ガス導管の最上流におけるガスの流速及び圧力は、時間を変数とする関数として取得される、ことを特徴とする。
【００１８】
更にまた、上記の発明において、好ましい態様は、前記流速決定ステップは、運動方程式に基づいて実行され、前記第混合比決定は、移流拡散方程式に基づいて実行される、ことを特徴とする。
【００１９】
更に、上記の発明において、その好ましい態様は、前記流速決定ステップは、モル密度保存式に基づいて実行される、ことを特徴とする。
【００２０】
上記の目的を達成するために、本発明の別の側面は、ガス導管内を流れる互いに種類の異なるガスの混合比を変更させた場合における、前記ガス導管内を流れるガスの物性変化挙動を推定する物性変化推定装置が、前記推定のための条件を受け付ける入力部と、前記推定のためのデータを格納する格納部と、前記入力部が受け付ける条件と前記格納部が格納するデータに基づいて、前記推定開始後の各時刻において、前記推定の対象範囲を所定長さで分割した第一分割領域毎に、前記混合比に基づいて、前記ガスの流速を求め、前記第一分割領域を更に分割した第二分割領域毎に、前記流速に基づいて前記混合比を求める推定処理部と、を有する、ことである。
【００２１】
上記の目的を達成するために、本発明の更に別の側面は、ガス導管内を流れる互いに種類の異なるガスの混合比を変更させた場合における、前記ガス導管内を流れるガスの物性変化挙動を推定する物性変化推定装置における物性変化推定方法において、前記物性変化推定装置は、前記推定のための条件を受け付ける入力部と、前記推定のためのデータを格納する格納部と、推定処理部とを備え、前記推定開始後の各時刻において、前記入力部が受け付ける条件と前記格納部が格納するデータに基づいて、前記推定処理部が、前記推定の対象範囲を所定長さで分割した第一分割領域毎に、前時刻で求めた前記混合比に基づいて、前記ガスの流速を求める工程と、前記推定処理部が、前記第一分割領域を更に分割した第二分割領域毎に、前記求めた流速に基づいて前記混合比を求める工程を有する、ことである。
【００２２】
本発明の更なる目的及び、特徴は、以下に説明する発明の実施の形態から明らかになる。
【図面の簡単な説明】
【００２３】
【図１】本発明を適用したガス導管内の物性変化推定プログラム等の対象となるガス導管網を例示した図である。
【図２】本実施の形態例に係る手法を説明するための図である。
【図３】図２に示した場合における、混合比の変更後の濃度変化挙動を例示した図である。
【図４】本実施の形態例に係る物性変化推定装置の構成図である。
【図５】本物性変化推定装置１０による処理の流れを例示したフローチャートである。
【図６】本推定計算で用いるメッシュを説明するための図である。
【図７】初期の界面形状Ｃ（ｘ）について説明するための図である。
【図８】本物性変化推定装置１０による推定計算の結果の一例を示した図である。
【発明を実施するための形態】
【００２４】
以下、図面を参照して本発明の実施の形態例を説明する。しかしながら、かかる実施の形態例が、本発明の技術的範囲を限定するものではない。なお、図において、同一又は類似のものには同一の参照番号又は参照記号を付して説明する。
【００２５】
図１は、本発明を適用したガス導管内の物性変化推定プログラム等の対象となるガス導管網を例示した図である。本発明を適用した物性変化推定プログラム等の実施の形態例は、図１に示すガス導管網１のような複数種類のガス輸送に用いられるガス導管網において、輸送するガスの置換、混合比の変更を行なった際の、導管内の物性変化の挙動を、導管をサイズの異なる２種類のメッシュに分割して数値解析し、従来法よりも柔軟でかつ適確な推定を、効率良く実行しようとするものである。
【００２６】
まず、図１に例示した推定の対象となる設備について説明する。図１に示すガス供給元２は、ガス導管網１内に設置されるガス製造工場（２ａ、２ｃ）、または、ガス導管網１外のガス製造工場やその他のガス供給源（２ｂ）であり、当該ガス導管網１について複数存在する。また、ガス消費設備３は、前記ガス供給元２から供給されるガスを使用する燃焼装置等の設備とその設置箇所を表しており、これも当該ガス導管網１に複数存在する（図１の例では、３ａ〜３ｅ）。
【００２７】
これらガス供給元２から送られるガスは、ガス導管網１を経由して各ガス消費設備３へ供給されるが、各ガス供給元２から送られるガスの性状及び量、各ガス消費設備３におけるガス消費量などにより、ガス導管網１内の各所において輸送されるガスの性状及び量は変化する。例えば、図１において、点線で囲われた部分イに着目すると、ガス導管網１のＡで示す部分では、ガス供給元２ａからのガス（ガスａ）とガス供給元２ｂからのガス（ガスａとは種類（性状）の異なるガスｂ）の両方が流れる構造となっている。従って、当該ガス導管１の使用方法によって、Ａの部分を流れるガスの性状及び量は変化することになる。
【００２８】
例えば、通常はＡの部分をガス供給元２ａからのガスａのみの輸送に使用しているが、時々、ガス供給元２ｂからのガスｂの託送に用いるという場合が考えられる。この場合、基本的にはガスａとガスｂを混合して輸送するが、何らかの理由によってその混合比を変更する場合や、輸送されるガスをガスａからガスｂに置換する場合や、ガスｂの成分が変化する場合等も考えられる。
【００２９】
このような、ガス導管１のＡの部分に流される２種類のガスの混合比の変更は、当該２種類のガスが合流する図１のＢで示す地点で発生することになるが、Ａの部分に流されるガスを使用する例えばガス消費設備３ａの地点（図１のＣ）においては、前記ガスの混合比の変更後、供給されるガスの性状が時間の経過とともに変化していくことになる。本実施の形態例に係る物性変化推定プログラム等では、このように、種類の異なるガスの混合比を変更した場合に、その下流におけるガス導管１内のガスの性状変化を推定しようとするものである。
【００３０】
図２は、本実施の形態例に係る手法を説明するための図である。図２に示すのは、２種理の異なるガスａ、ｂをそれぞれ輸送する二つの導管が合流し、合流後のガス導管１がガス消費設備３ａにつながっている場合、すなわち、図１のイ部分に示した場合の例である。本実施の形態例に係る物性変化推定は、２種類のガスが合流する地点でそれらガスの混合比を変更させた際に、合流後のガス導管１内における各ガスの濃度変化を求めるものであり、より具体的には、合流後のガス導管１内に流すガスａ＋ガスｂの中のガスｂの濃度Ｃ（以下、ガスｂの濃度とは合流後のガス導管１内におけるガスａ＋ガスｂの中のガスｂの濃度を指す）をαからβに変更した場合に、混合発生点（ガス導管１の軸方向の距離をｘとするとｘ＝０の地点、前記合流地点）から下流に任意の地点の、濃度の変更から任意の時間経過後における、ガスｂの濃度Ｃを求めるものである。例えば、この手法により、混合したガスを使用するガス消費設備３ａにおける（図２の消費点（ｘ＝Ｌの地点）における）ガス性状の時間的変化を求めることができる。具体的な推定手法については後述するが、当該推定では、上記距離ｘ及び経過時間ｔの関数であるガス流速ｕ（ｘ,ｔ）及びガス圧力Ｐ（ｘ,ｔ）を求めるステップと上記濃度Ｃを求めるステップが交互に実行される。
【００３１】
図３は、図２に示した場合における、混合比の変更後の濃度変化挙動を例示した図である。図３の（ａ）の左側は、前記混合発生点（ｘ＝０）で、ガスｂの濃度Ｃをαからβに変更した直後（ｔ＝ｔ０）における、混合発生点前後のガスｂの濃度Ｃを表している。この濃度変更後、時間の経過とともに混合ガスはガス導管１内を下流に移動していくが、その間にガスの拡散現象により、混合比が変更された前後の混合ガス間の境目（以下、界面と呼ぶ）はぼやけたものになっていく。すなわち、当該界面はガス導管１の軸方向にある幅を有するようになる。従って、ガスｂの濃度変更後、ある時間（Δｔ）経過後には、前記界面部分は混合発生点からΔｘだけ下流に移動し、その前後におけるガスｂの濃度は、例えば、図３の（ａ）の右側に示すような状態となる。図に示されるように、図３の（ａ）の左側においてロで示される部分（界面）は、前記拡散現象に伴い傾斜をもつようになる。
【００３２】
図３の（ｂ）は、図２における消費点（ｘ＝Ｌ）におけるガスｂの濃度の時間的変化を例示している。図に示すように、前記界面部分の通過に伴い、消費点ではガスｂの濃度がある時間をかけてαからβに変化していくと考えられる。この濃度変化の挙動を適確に捉えることは、ガス消費設備３にとって重要なことである。なお、図３の例では、混合発生点における濃度変化が瞬時にαからβに変化するという、界面の形状が単純なケースであるが、後述するように本推定方法では、初期時の界面形状がより複雑なケースであっても対応可能である。
【００３３】
本実施の形態例に係る物性変化推定プログラム等は、このようなガス混合比の変更後の挙動を、すなわち、ガス導管内の任意位置、任意時間での濃度Ｃを、ガス導管を分割したメッシュ毎の数値計算により、適確に把握しようとするものである。
【００３４】
なお、ガス濃度が１００％、すなわち、そのガスだけの状態、の場合の濃度Ｃを１で表すものとすると、前記ガスｂの濃度変更において、αが０、βが１である場合には、ガスａからガスｂへ置換することを意味することになる。ここでは、かかる置換も含めてガスの混合比の変更と称することとする。
【００３５】
なお、本発明におけるガス（例えば、ガスａ、ガスｂ）とは、単成分で構成されているガス及び複数成分で構成されているガスのいずれをも意味しており、本実施の形態例におけるガスａ、ｂについては、単成分又は複数成分のどちらであっても良い。さらに、ガスｂが複数成分で構成されている場合で、ガスｂ中の成分変化により、ガスｂからガスｂ’へ変化する場合には、そのガスｂを構成している単成分毎に、ガスａとの合流後のガスａ＋ガスｂ’中の当該単成分毎に濃度の変化を本物性変化推定方法にて推定し、かかるのちにその成分毎の濃度から物性を求めることで、ガスｂからガスｂ’への変化に対しての、合流後のガス導管内のガス物性変化を推定できる。
【００３６】
図４は、本実施の形態例に係る物性変化推定装置の構成図である。図４に示す物性変化推定装置１０が、本実施の形態例に係る物性変化推定装置であり、具体的には、パーソナルコンピューターなどのコンピュータシステムで構成することができる。入力部１１は、前述した混合比の変更後の挙動を推定するために必要な各種推定条件の入力を行なうための部分である。
【００３７】
例えば、混合させるガスａ、ｂ、初期状態での界面の形状、すなわち、ガス導管の各位置ｘにおける濃度Ｃ（ｘ）、及び、ｘ＝０におけるガス流速ｕ及びガス圧力Ｐなどが当該入力部１１から物性変化推定装置１０に入力される。なお、この入力部１１は、当該装置の操作者が入力操作を行なうためのキーボード、マウス、ディスプレイ等、及び前記入力操作を促すための画面表示や入力されたデータの取得の処理を記述したプログラム、当該プログラムに従って処理を実行する制御装置（ＣＰＵ）等で構成される。
【００３８】
次に、推定処理部１２は、前記入力部１１から入力された推定条件等に基づいて、混合比の変更後のガス導管１内におけるガス濃度変化を推定する部分である。推定処理の具体的な内容については後述するが、当該推定処理部１２が実行する数値解析の手法に大きな特徴を有している。なお、当該推定処理部１２は、処理内容を記述したプログラム、当該プログラムに従って処理を実行する制御装置（ＣＰＵ）等で構成される。
【００３９】
データ格納部１３は、前記推定処理部１２が処理を実行するために必要な各種データを格納するデータ記憶手段であり、推定処理の度に操作者によって入力部１１から入力しなくてはならないデータ以外のデータを格納する。例えば、本物性変化推定装置１０が処理の対象とするガス導管網１の各所における導管の寸法（内径、長さ等）、輸送する各ガスの物性値（動粘度係数又は静粘度係数、等）などが予め格納される。なお、当該データ格納部１３は、ＨＤＤ等で構成される。
【００４０】
また、出力部１４は、前記推定処理部１２が行なった推定処理の結果を出力する部分であり、ガス導管１の所定の箇所又は所定の時間におけるガス濃度変化の推定結果を、表形式やグラフ形式で出力する。なお、この出力部１４は、ディスプレイ、プリンタへの出力インターフェース等、及び前記ディスプレイへの画面表示や前記プリンタへの出力の処理を記述したプログラム、当該プログラムに従って処理を実行する制御装置（ＣＰＵ）等で構成され、前記推定結果は、画面表示や紙等への印刷という形で出力される。
【００４１】
本物性変化推定装置１０は、以上説明したような構成を有し、以下に示すような内容で処理を実行する。図５は、本物性変化推定装置１０による処理の流れを例示したフローチャートである。
【００４２】
当該装置を用いてガス導管１内のガス混合比を変更させた場合のその後の挙動を推定しようとする操作者は、まず、前述した入力部１１を用いて推定処理に必要な各種の条件を入力する（ステップＳ１）。なお、前述のとおり、本物性変化推定装置１０では、図２に示した状態を前提とした推定処理を行なうものとする。
【００４３】
この推定条件の入力では、ガス導管１で輸送されるガスａ、ｂの種類、混合比変更前後のガスｂの濃度（α、β）及び初期状態における界面の形状（Ｃ（ｘ））、初期位置（ｘ＝０）における混合ガスのガス流速ｕ及びガス圧力Ｐ、ガス導管１の位置（ガス導管１を特定する情報）、ガス導管１内を流れるガスの温度等が入力され、ガスの濃度変化挙動を把握したい特定の位置（例えば、図２の消費点等）がある場合にはその位置を特定する情報を入力する。なお、データ格納部１３に、ガス導管網１の各所における導管の寸法及び物性値、各ガスの物性値などが格納されていない場合には、これらのうち必要なデータについては操作者によって入力される。また、本推定処理における数値計算では、後述するように、ガス導管１をメッシュに分割してメッシュ毎に逐次計算を行っていくが、当該メッシュのサイズを変更可能とする場合には、操作者がメッシュサイズも合わせて入力する。
【００４４】
次に、推定処理部１２は、当該推定処理を行うために必要な各種データを前記データ格納部１３から取得する（ステップＳ２）。具体的には、推定対象とするガス導管１について、その長さＬ、各部の内径Ｄ、及び内壁の摩擦係数λ、また、ガスａ及びガスｂの気体分子量Ｍａ、Ｍｂ、普遍気体定数Ｒ０、ガスの静粘度係数μ等のデータを取得する。なお、ガス導管１の各部（各メッシュ）の内径Ｄ及び摩擦係数λは、下記推定計算の過程で必要な際に取得するようにしてもよい。
【００４５】
次に、推定処理部１２は、以上入力されたデータと取得したデータを用いて、推定計算を実行する（ステップＳ３）。まず、当該推定計算で用いる式について説明する。当該推定計算では、以下に示す７つの式を用いて、各時間ｔ及び各位置ｘにおけるガス濃度Ｃを求める。
【００４６】
【数１】
【００４７】
【数２】
【００４８】
【数３】
【００４９】
なお、上式においてＫは、実験から以下のように与えられる。
Ｋ＝Ａ３ｕＤＲｅ−０．１２５ （４−１）
（分子量の式）
Ｍ＝Ｍａ＋（Ｍｂ−Ｍａ）Ｃ （５）
（気体定数の式）
Ｒ＝Ｒ０／Ｍ （６）
（モル密度と密度の関係式）
Ｘ＝ρ／Ｍ （７）
なお、上記（１）−（７）式において、各記号が表わすものは以下の通りである。
ｔ：時間 ［ｓ］
ｘ：流れ方向の座標（位置） ［ｍ］
ｕ：ガス流速 ［ｍ／ｓ］
Ｐ：ガス圧力 ［Ｐａ］
ρ：ガスの密度 ［ｋｇ／ｍ３］
λ：ガス導管内壁の摩擦係数 ［‐］
Ｄ：ガス導管の内径 ［ｍ］
Ｘ：ガスｂのモル密度 ［ｍｏｌ／ｍ３］
Ｒ：気体定数 ［Ｊ／（ｋｇ・Ｋ）］
θ：ガスの温度 ［Ｋ］
Ｃ：ガスｂの濃度 ［‐］
Ｋ：軸方向（流れ方向）の拡散係数 ［ｍ２／ｓ］
Ｒｅ：レイノズル数 ［‐］
Ｍ：ガス分子量 ［ｋｇ／ｍｏl］
Ｍａ：ガスａの分子量 ［ｋｇ／ｍｏl］
Ｍｂ：ガスｂの分子量 ［ｋｇ／ｍｏl］
Ｒ０：普遍気体定数 ８．１３４［Ｊ／（ｋｇ・Ｋ）］
Ａ３：比例定数（例えば、２．５９）
なお、レイノズル数Ｒｅは、上述したガスの静粘度係数μ、ガス導管の内径Ｄ、及びガスの密度ρから求めることができる。
【００５０】
以上の式において、時間ｔ及び位置ｘに対して未知数である変数は、Ｃ、Ｍ、Ｐ、ｕ、ρ、Ｒ、Ｘの７つであり、推定処理部１２は、この７つの未知数を上記７つの式（（１）−（７）式）を解くことによって求める。すなわち、推定計算を実行する。
【００５１】
なお、上記（２）式において、一般的な質量保存側に係る式を用いずに、モル密度保存式を用いたことが本装置の一つの特徴である。一般的な質量保存側に係る式を用いた場合、ガス密度が用いられることになるが、この密度の変化の要因としては、ガス圧力の変化に起因するものと、ガス成分（濃度）の変化に起因するものの二つが存在する。ここで対象としている混合ガスは、上述の通り、界面付近で濃度が大きく変化するので、上記ガス成分（濃度）の変化に起因する密度の変化が大きく、一般的な質量保存側に係る式を用いた場合、数値解析で用いる差分法において発散が生じてしまう可能性が高い。そこで、本装置では、このような発散を防止するため、モル密度保存式を用いている。
【００５２】
次に、推定処理部１２が行う上記７つの式を用いた推定計算の具体的な内容について説明する。当該推定計算はいわゆる数値解析によって行うが、その際に用いるメッシュ（解析対象を分割したもの）について、まず、説明する。
【００５３】
図６は、本推定計算で用いるメッシュを説明するための図である。図６において、ガス導管１の長さＬの部分（ｘ＝０からｘ＝Ｌの部分）が推定対象を表している。そして、図中の矢印がガスの流れ方向を示しており、従って、推定対象において左側が上流となる。
【００５４】
当該推定対象をｎ個に等分割した各領域（第一分割領域）が、本実施の形態例における１種類目のメッシュ（以下、粗メッシュと呼ぶ）であり、図中のＥ０〜Ｅｎ−１がそれに相当する。かかる粗メッシュは、ガス流速ｕ及びガス圧力Ｐを求める際に用いるものであり、例えば、推定対象の長さＬが１０ｋｍ程度である場合に、これに限定されるものではないが、一つの粗メッシュの長さを１００ｍ程度とすることができる。図６に示す例では、ある時間ｔにおいて、最上流の粗メッシュＥ０に対して、その左端のガス流速ｕ０及びガス圧力Ｐ０を用いて計算がなされ、Ｅ０の右端のガス流速ｕ１及びガス圧力Ｐ１が求められる。以下、同様に下流に向かって各粗メッシュに対して計算がなされる。
【００５５】
次に、上述した一つの粗メッシュをさらにｋ個に等分割した各領域（第二分割領域）が、本実施の形態例における２種類目のメッシュ（以下、細メッシュと呼ぶ）であり、図中のｅｉ，０〜ｅｉ，ｋ−１がそれに相当する。かかる細メッシュは、ガスｂの濃度Ｃを求める際に、すなわち界面形状を求める際に、用いるものであり、例えば、一つの粗メッシュの長さが１００ｍ程度である場合に、これに限定されるものではないが、一つの細メッシュの長さを１〜１０ｍ程度とすることができる。推定計算では、ある時間ｔにおいて、各細メッシュについて、それぞれ、ガスｂの濃度Ｃ（図中のＣｉ，０〜Ｃｉ，ｋ−１）が求められることになる。なお、後述するが、この細メッシュを用いた計算は全ての粗メッシュについて行う必要はなく、界面が存在する粗メッシュ及びその近傍について行われる。また、当該細メッシュ毎に、前述したガス導管１の各仕様が格納されていることが望ましい。
【００５６】
このように、本推定方法では、２種類のメッシュを使用して計算を行う。これは、気体の特性に基づくものであり、ガス成分が変化する前記界面付近においても、ガスの流速ｕ及び圧力Ｐは大きな変化をしない一方、ガス濃度Ｃは急激に変化するので、前者については、比較的大きなメッシュでの計算でも精度上問題がなく、後者については、精度良く状態を把握するために細かいメッシュでの計算が必要となる。
【００５７】
以下、具体的な計算手順について説明する。図５に戻って、まず、初期時におけるガス濃度Ｃを決定する。（ステップＳ３−１）。図７は、初期の界面形状Ｃ（ｘ）について説明するための図である。通常の推定においては、ガス混合地点からのガス混合時点以降を対象とするので、当該初期時においては、通常、界面が粗メッシュＥ０に存在する。従って、上記界面が存在する箇所でのガス濃度Ｃの決定は、図７に例示するように、入力され界面の形状Ｃ（ｘ）（図中の曲線）に従って、各細メッシュのＣ値を決定する。図７の例では、ｅ０，０、ｅ０，１等の細メッシュについて、対応するＣ０，０、Ｃ０，１等の値が、図中のプロット点×のように決定される。
【００５８】
そして、上記界面が存在する箇所については、上記各細メッシュのＣの値から、対応する粗メッシュの濃度Ｃを決定する。具体的には、それらの平均値から粗メッシュの値を求める。また、界面が存在しない箇所については、界面の下流側ではＣの値がαであり、その上流側ではβであるので、粗メッシュ単位でそのように決定する。
【００５９】
次に、時間変数ｔをｔ０に設定後（ステップＳ３−２）、初期位置（ｘ＝０）における圧力Ｐ０及び流速ｕ０を、前記入力値に従って設定する（ステップＳ３−３）。そして、この時点におけるガス流速ｕ及びガス圧力Ｐを算出する（ステップＳ３−４）。ここで、前述した（５）式において、上記決定された濃度ＣによりＣの値が既知となるので、Ｍの値が決定され、そのＭの値により、上記（６）式からＲの値が求められる。この時点で未知数は、ｕ、Ｐ、Ｘ、ρの４つであり、それらの関係式が上記（１）−（３）及び（７）式として４つ与えられるので、これらの関係式を用いて上記４つの未知数の値が決定される。
【００６０】
具体的には、スタガート格子差分法などを用いて各粗メッシュ毎に数値計算を行って求める。図６に示した例では、前述のようにｕ０、Ｐ０の値が与えられているので、Ｅ０について数値計算を行いｕ１、Ｐ１の値等を算出する。以下、同様に、下流方向に計算を行って、ｕ２、Ｐ２〜ｕｎ、Ｐｎの値を求める。
【００６１】
次に、時間変数ｔをｔ＋Δｔとし（ステップＳ３−５）、次の時間（Δｔ経過後）におけるガス濃度界面の中心位置を計算する（ステップＳ３−６）。ガス濃度界面の中心位置ｘ界面は、前時刻における界面の中心位置ｘ前界面と界面が存在した粗メッシュの流速ｕ界面が存在する粗メッシュから、ｘ界面＝ｘ前界面＋Δｔ×ｕ界面が存在する粗メッシュとして計算される。
【００６２】
その結果から、界面が推定対象の最下流の位置を通過したかを判断し（ステップＳ３−７）、通過していない場合には（ステップＳ３−７のＮｏ）、この時刻におけるガス濃度Ｃと圧力Ｐ及び流速ｕの算出を行う。
【００６３】
具体的には、まず、上記界面の中心位置近傍の細メッシュについて濃度Ｃの値を算出する。すなわち、界面形状を求める（ステップＳ３−８）。ここでは、前時刻におけるガス流速ｕが求められているので、その値を前述した（４）及び（４−１）式に代入してＣの値を求めることができる。ここでも、具体的にはスタガート格子差分法などを用いて各細メッシュ毎に数値計算を行って求める。その後、上記ステップＳ３−１と同様にして、粗メッシュ毎のガス濃度Ｃを決定する（ステップＳ３−９）。
【００６４】
次に、処理が前記ステップＳ３−４に戻って、上述の通り、この時刻における圧力Ｐ及び流速ｕが求められる。
【００６５】
以降、上記ステップＳ３−７において、界面が推定対象の最下流の位置を通過していると判断されるまで、同様の処理が繰り返される。そして、界面が通過していると判断された場合には（ステップＳ３−７のＹｅｓ）、推定計算（Ｓ３）を終了し、推定結果の出力を行う（ステップＳ４）。
【００６６】
このように、本推定計算では、各時間（各時刻）において、前時刻で求められた粗メッシュ毎のガス流速ｕを用いて界面近傍について細メッシュ毎にガス濃度Ｃを求め（すなわち、界面形状を求め）、その求められた濃度Ｃに基づいて粗メッシュ毎に流速ｕを求める。そして、当該処理を時間経過に従って順次繰り返して実行することにより、任意の時刻ｔ、任意の位置ｘでのガス濃度Ｃを求めることができる。なお、上記説明では、各時刻において先に濃度Ｃを求めその後に流速ｕを求めたが、前時刻における濃度Ｃを用いて先に流速ｕを求め、その後に濃度Ｃを求めるようにしてもよい。
【００６７】
推定結果の出力では、用途に応じて様々な形式での出力が可能であるが、例えば、ガス消費設備３の位置での、界面の通過時における濃度Ｃの時間変化をグラフ形式で出力する、ことが可能である。
【００６８】
図８は、本物性変化推定装置１０による推定計算の結果の一例を示した図である。図８に示す例は、一つの粗メッシュにおける各時間での界面（におけるＣの値）をプロットしたものである。実線は界面が当該メッシュに入ってから１０秒後の状態を示しており、点線は、２０秒後の状態を示している。１０秒後の状態と２０秒後の状態を比較すると、後者では界面の傾斜がなだらかになり、ｘ方向に拡散現象が起こっていることがわかる。また、３０秒後では、当該領域に界面がなくなり、界面が通過していることがわかる。
【００６９】
なお、上述の説明では、初期位置（ｘ＝０）における混合ガスのガス流速ｕ及びガス圧力Ｐを推定条件として入力したが、これらの初期値を時間ｔによって変化する変数ｕ（ｔ）及びＰ（ｔ）として与えるようにしてもよい。この場合には、各時間ｔにおけるｕ０及びＰ０をｕ（ｔ）及びＰ（ｔ）によって得られる値とすることにより、上述の場合と同様に計算を行うことができる。従って、この場合には、図５に示した処理手順において、ステップＳ３−９の後に、処理がステップＳ３−３に移行することになる。
【００７０】
また、図２等においては、ｘ＝０を混合発生点、ｘ＝Ｌを消費点として、これらの間を推定対象としたが、これに限定されることはなく、本推定計算は任意の領域に適用することができる。
【００７１】
なお、本物性変化推定装置１０は、スタンドアローンのシステムとしたが、ネットワークを介して、操作者による入力や結果の出力が行なわれるようにしてもよい。また、ガス導管網１の各所で測定されるガス流速などのリアルタイムデータを推定の入力値として用いるべく、それら測定装置等と通信可能にして用いてもよい。
【００７２】
以上説明したように、本実施の形態例に係る物性変化推定装置１０では、推定対象範囲をメッシュに分割して数値計算するのでメッシュ毎にガス導管１の仕様を設定可能であり、また、初期時の界面形状をＣ（ｘ）で与えられるので複雑な界面形状の場合にも対応可能である。さらに、上述のように、初期位置のガス流速ｕ及びガス圧力Ｐを時間ｔの関数として与えることもできる。従って、本装置１０により、従来の装置よりも適用範囲が広くなって、実際に存在する様々なケースのシミュレーションを適確に行うことが可能となる。
【００７３】
また、上記（４−１）式ではガス導管１内の流れが乱流である場合の式を採用しており、流れが層流である場合に限らず適用することができる。
【００７４】
また、前述の通り、連続性を担保する関係式としてモル密度保存式を用い、数値計算上の不具合を防止している。
【００７５】
また、本物性変化推定装置１０では、推定対象を粗メッシュと細メッシュに分割し、気体の特性に基づいて、細メッシュ毎の計算が必要な場合のみ細メッシュでの計算を実行する。従って、ガス導管網のように推定対象が長い場合に解析上の負荷を大幅に軽減することが可能である。
【００７６】
このように、本物性変化推定装置１０によって、ガスの物性が変化した場合の挙動を、適確にかつ低コストで把握することが可能となり、ガス供給の信頼性をより向上させることができる。
【００７７】
本発明の保護範囲は、上記の実施の形態に限定されず、特許請求の範囲に記載された発明とその均等物に及ぶものである。
【符号の説明】
【００７８】
１ ガス導管（網）
２ ガス供給元
３ ガス消費設備
１０ 物性変化推定装置
１１ 入力部
１２ 推定処理部
１３ データ格納部
１４ 出力部
【技術分野】
【０００１】
本発明は、ガスの輸送導管を成分の異なる複数種類のガス輸送に使用する場合における、導管内の物性変化推定プログラム等に関し、特に、ガスの混合比が変化した場合の導管下流側における濃度変化を、低コストでかつ適確に把握することのできる物性変化推定プログラム等に関する。
【背景技術】
【０００２】
都市ガスは、地域に張り巡らせたガス導管網を経由して、ガス製造工場等のガス供給元からガスを使用するガス消費設備等へ供給される。また、供給する都市ガスは、複数の地域から産出され組成のバラツキがある天然ガスに、ＬＰＧ、ブタン、プロパンなどを混合して熱量が一定になるように調整がなされる。しかし、ガスの輸送導管の有効活用や、ガス輸送導管を持たない事業者への公平な参入機会の創出などの要請から、複数の地域から産出した天然ガスをそのままガス導管網に導入するなど、熱量がそれぞれ異なる複数種類のガスがガス導管網に導入されることもある。
【０００３】
また、今後は、益々エネルギーの多様化の要請が高まり、ガス導管網をこれまで以上に性状の異なる異種ガスの輸送に使用することが考えられる。さらに、ガス供給の新規参入者等のガス輸送を代行する託送に、ガス導管網が利用される頻度が多くなることが予想され、これによってもガス導管網において、性状のかなり異なるガスの輸送を行なう機会が増えると考えられる。
【０００４】
このような状況では、ガス導管内を流れるガスの性状が変化する機会が増え、そのような場合には、供給されるガスを消費する設備等においては、ガスの性状が大きく変化する場合があり、その変化挙動を適確に把握しなければ設備の運転に支障をきたしたり、危険が生じてしまう場合もある。
【０００５】
そこで、下記特許文献１では、輸送するガスの置換、混合比の変更による導管内の物性変化挙動を所定の拡散方程式を用いて求めることについて提案がなされている。
【０００６】
また、下記特許文献２には、上記拡散方程式における拡散定数を、導管内のガス流れが属する各領域毎に定められた所定の関係式によって決定される技術が開示されている。
【０００７】
また、下記非特許文献１では、配管破断事故後に生じる自然循環の特性を調べるための、２成分気体の拡散と自然循環に関する研究について示され、下記非特許文献２には、ガスパイプラインの非定常流解析シミュレーションについて示されている。
【先行技術文献】
【特許文献】
【０００８】
【特許文献１】特開２００６−３３６７９５号公報
【特許文献２】特開２００８−２４０８６４号公報
【非特許文献】
【０００９】
【非特許文献１】武田哲明、菱田誠、「２成分気体の拡散と自然循環に関する研究」、日本機械学会論文集、１９８９年８月、５５巻、５１６号、ｐ．２７９５−２７９９
【非特許文献２】山根総一郎、佐藤律夫、「ガスパイプライン異種ガス対応非定常流解析シミュレータの開発」、ＪＦＥ技報、２００６年１１月、Ｎｏ．１４、ｐ．６０−６４
【発明の概要】
【発明が解決しようとする課題】
【００１０】
しかしながら、上記特許文献１、２に記載の方法は、導管内の流速、導管の内径等の仕様が一定である場合、推定開始時の異種ガスの界面が単純な形状をしている場合に適用可能であり、適用範囲が限られるという課題があった。
【００１１】
また、上記非特許文献１に記載の方法は、対象範囲が層流であり導管内の流れが乱流である場合には対応できず、また、管の長さが短い場合を対象としているので、対象が上述したガスの輸送導管（網）のように長い場合には、解析のための手法が計算負荷の面で向いていない。
【００１２】
また、上記非特許文献２に記載の方法では、気体の移流のみによる効果が反映され、気体の拡散が考慮されていないので、推定結果の精度に課題がある。
【００１３】
そこで、本発明の目的は、ガスの輸送導管を成分の異なる複数種類のガス輸送に使用する場合における導管内の物性変化推定プログラムであって、ガスの混合比が変化した場合の導管下流側における濃度変化を、低コストでかつ適確に把握することのできると共に、様々な条件下における実際のケースに使用可能な物性変化推定プログラム、等を提供することである。
【課題を解決するための手段】
【００１４】
上記の目的を達成するために、本発明の一つの側面は、ガス導管内を流れる互いに種類の異なるガスの混合比を変更させた場合における、前記ガス導管内を流れるガスの物性変化挙動を推定する処理をコンピュータに実行させる物性変化推定プログラムが、前記推定開始後の各時刻において、前記推定処理の対象範囲を所定長さで分割した第一分割領域毎に、前記混合比に基づいて、前記ガスの流速を求める流速決定ステップと、前記第一分割領域を更に分割した第二分割領域毎に、前記流速に基づいて前記混合比を求める（すなわち、界面が存在する領域については界面形状を求める）混合比決定ステップとを、前記コンピュータに実行させる、ことである。
【００１５】
更に、上記の発明において、その好ましい態様は、前記混合比決定ステップは、前記混合比が変化する部分である界面が存在する領域について実行される、ことを特徴とする。
【００１６】
また、上記の発明において、好ましい態様は、更に、少なくとも、前記第二分割領域毎の前記ガス導管の仕様と、前記推定開始時における前記ガス導管の各位置における前記混合比と、前記ガス導管の最上流における前記ガスの流速及び圧力を取得するステップを、前記コンピュータに実行させ、前記流速決定ステップと前記混合比決定ステップにおいて、当該取得された情報が用いられる、ことを特徴とする。
【００１７】
更に、上記の発明において、一つの態様は、前記ガス導管の最上流におけるガスの流速及び圧力は、時間を変数とする関数として取得される、ことを特徴とする。
【００１８】
更にまた、上記の発明において、好ましい態様は、前記流速決定ステップは、運動方程式に基づいて実行され、前記第混合比決定は、移流拡散方程式に基づいて実行される、ことを特徴とする。
【００１９】
更に、上記の発明において、その好ましい態様は、前記流速決定ステップは、モル密度保存式に基づいて実行される、ことを特徴とする。
【００２０】
上記の目的を達成するために、本発明の別の側面は、ガス導管内を流れる互いに種類の異なるガスの混合比を変更させた場合における、前記ガス導管内を流れるガスの物性変化挙動を推定する物性変化推定装置が、前記推定のための条件を受け付ける入力部と、前記推定のためのデータを格納する格納部と、前記入力部が受け付ける条件と前記格納部が格納するデータに基づいて、前記推定開始後の各時刻において、前記推定の対象範囲を所定長さで分割した第一分割領域毎に、前記混合比に基づいて、前記ガスの流速を求め、前記第一分割領域を更に分割した第二分割領域毎に、前記流速に基づいて前記混合比を求める推定処理部と、を有する、ことである。
【００２１】
上記の目的を達成するために、本発明の更に別の側面は、ガス導管内を流れる互いに種類の異なるガスの混合比を変更させた場合における、前記ガス導管内を流れるガスの物性変化挙動を推定する物性変化推定装置における物性変化推定方法において、前記物性変化推定装置は、前記推定のための条件を受け付ける入力部と、前記推定のためのデータを格納する格納部と、推定処理部とを備え、前記推定開始後の各時刻において、前記入力部が受け付ける条件と前記格納部が格納するデータに基づいて、前記推定処理部が、前記推定の対象範囲を所定長さで分割した第一分割領域毎に、前時刻で求めた前記混合比に基づいて、前記ガスの流速を求める工程と、前記推定処理部が、前記第一分割領域を更に分割した第二分割領域毎に、前記求めた流速に基づいて前記混合比を求める工程を有する、ことである。
【００２２】
本発明の更なる目的及び、特徴は、以下に説明する発明の実施の形態から明らかになる。
【図面の簡単な説明】
【００２３】
【図１】本発明を適用したガス導管内の物性変化推定プログラム等の対象となるガス導管網を例示した図である。
【図２】本実施の形態例に係る手法を説明するための図である。
【図３】図２に示した場合における、混合比の変更後の濃度変化挙動を例示した図である。
【図４】本実施の形態例に係る物性変化推定装置の構成図である。
【図５】本物性変化推定装置１０による処理の流れを例示したフローチャートである。
【図６】本推定計算で用いるメッシュを説明するための図である。
【図７】初期の界面形状Ｃ（ｘ）について説明するための図である。
【図８】本物性変化推定装置１０による推定計算の結果の一例を示した図である。
【発明を実施するための形態】
【００２４】
以下、図面を参照して本発明の実施の形態例を説明する。しかしながら、かかる実施の形態例が、本発明の技術的範囲を限定するものではない。なお、図において、同一又は類似のものには同一の参照番号又は参照記号を付して説明する。
【００２５】
図１は、本発明を適用したガス導管内の物性変化推定プログラム等の対象となるガス導管網を例示した図である。本発明を適用した物性変化推定プログラム等の実施の形態例は、図１に示すガス導管網１のような複数種類のガス輸送に用いられるガス導管網において、輸送するガスの置換、混合比の変更を行なった際の、導管内の物性変化の挙動を、導管をサイズの異なる２種類のメッシュに分割して数値解析し、従来法よりも柔軟でかつ適確な推定を、効率良く実行しようとするものである。
【００２６】
まず、図１に例示した推定の対象となる設備について説明する。図１に示すガス供給元２は、ガス導管網１内に設置されるガス製造工場（２ａ、２ｃ）、または、ガス導管網１外のガス製造工場やその他のガス供給源（２ｂ）であり、当該ガス導管網１について複数存在する。また、ガス消費設備３は、前記ガス供給元２から供給されるガスを使用する燃焼装置等の設備とその設置箇所を表しており、これも当該ガス導管網１に複数存在する（図１の例では、３ａ〜３ｅ）。
【００２７】
これらガス供給元２から送られるガスは、ガス導管網１を経由して各ガス消費設備３へ供給されるが、各ガス供給元２から送られるガスの性状及び量、各ガス消費設備３におけるガス消費量などにより、ガス導管網１内の各所において輸送されるガスの性状及び量は変化する。例えば、図１において、点線で囲われた部分イに着目すると、ガス導管網１のＡで示す部分では、ガス供給元２ａからのガス（ガスａ）とガス供給元２ｂからのガス（ガスａとは種類（性状）の異なるガスｂ）の両方が流れる構造となっている。従って、当該ガス導管１の使用方法によって、Ａの部分を流れるガスの性状及び量は変化することになる。
【００２８】
例えば、通常はＡの部分をガス供給元２ａからのガスａのみの輸送に使用しているが、時々、ガス供給元２ｂからのガスｂの託送に用いるという場合が考えられる。この場合、基本的にはガスａとガスｂを混合して輸送するが、何らかの理由によってその混合比を変更する場合や、輸送されるガスをガスａからガスｂに置換する場合や、ガスｂの成分が変化する場合等も考えられる。
【００２９】
このような、ガス導管１のＡの部分に流される２種類のガスの混合比の変更は、当該２種類のガスが合流する図１のＢで示す地点で発生することになるが、Ａの部分に流されるガスを使用する例えばガス消費設備３ａの地点（図１のＣ）においては、前記ガスの混合比の変更後、供給されるガスの性状が時間の経過とともに変化していくことになる。本実施の形態例に係る物性変化推定プログラム等では、このように、種類の異なるガスの混合比を変更した場合に、その下流におけるガス導管１内のガスの性状変化を推定しようとするものである。
【００３０】
図２は、本実施の形態例に係る手法を説明するための図である。図２に示すのは、２種理の異なるガスａ、ｂをそれぞれ輸送する二つの導管が合流し、合流後のガス導管１がガス消費設備３ａにつながっている場合、すなわち、図１のイ部分に示した場合の例である。本実施の形態例に係る物性変化推定は、２種類のガスが合流する地点でそれらガスの混合比を変更させた際に、合流後のガス導管１内における各ガスの濃度変化を求めるものであり、より具体的には、合流後のガス導管１内に流すガスａ＋ガスｂの中のガスｂの濃度Ｃ（以下、ガスｂの濃度とは合流後のガス導管１内におけるガスａ＋ガスｂの中のガスｂの濃度を指す）をαからβに変更した場合に、混合発生点（ガス導管１の軸方向の距離をｘとするとｘ＝０の地点、前記合流地点）から下流に任意の地点の、濃度の変更から任意の時間経過後における、ガスｂの濃度Ｃを求めるものである。例えば、この手法により、混合したガスを使用するガス消費設備３ａにおける（図２の消費点（ｘ＝Ｌの地点）における）ガス性状の時間的変化を求めることができる。具体的な推定手法については後述するが、当該推定では、上記距離ｘ及び経過時間ｔの関数であるガス流速ｕ（ｘ,ｔ）及びガス圧力Ｐ（ｘ,ｔ）を求めるステップと上記濃度Ｃを求めるステップが交互に実行される。
【００３１】
図３は、図２に示した場合における、混合比の変更後の濃度変化挙動を例示した図である。図３の（ａ）の左側は、前記混合発生点（ｘ＝０）で、ガスｂの濃度Ｃをαからβに変更した直後（ｔ＝ｔ０）における、混合発生点前後のガスｂの濃度Ｃを表している。この濃度変更後、時間の経過とともに混合ガスはガス導管１内を下流に移動していくが、その間にガスの拡散現象により、混合比が変更された前後の混合ガス間の境目（以下、界面と呼ぶ）はぼやけたものになっていく。すなわち、当該界面はガス導管１の軸方向にある幅を有するようになる。従って、ガスｂの濃度変更後、ある時間（Δｔ）経過後には、前記界面部分は混合発生点からΔｘだけ下流に移動し、その前後におけるガスｂの濃度は、例えば、図３の（ａ）の右側に示すような状態となる。図に示されるように、図３の（ａ）の左側においてロで示される部分（界面）は、前記拡散現象に伴い傾斜をもつようになる。
【００３２】
図３の（ｂ）は、図２における消費点（ｘ＝Ｌ）におけるガスｂの濃度の時間的変化を例示している。図に示すように、前記界面部分の通過に伴い、消費点ではガスｂの濃度がある時間をかけてαからβに変化していくと考えられる。この濃度変化の挙動を適確に捉えることは、ガス消費設備３にとって重要なことである。なお、図３の例では、混合発生点における濃度変化が瞬時にαからβに変化するという、界面の形状が単純なケースであるが、後述するように本推定方法では、初期時の界面形状がより複雑なケースであっても対応可能である。
【００３３】
本実施の形態例に係る物性変化推定プログラム等は、このようなガス混合比の変更後の挙動を、すなわち、ガス導管内の任意位置、任意時間での濃度Ｃを、ガス導管を分割したメッシュ毎の数値計算により、適確に把握しようとするものである。
【００３４】
なお、ガス濃度が１００％、すなわち、そのガスだけの状態、の場合の濃度Ｃを１で表すものとすると、前記ガスｂの濃度変更において、αが０、βが１である場合には、ガスａからガスｂへ置換することを意味することになる。ここでは、かかる置換も含めてガスの混合比の変更と称することとする。
【００３５】
なお、本発明におけるガス（例えば、ガスａ、ガスｂ）とは、単成分で構成されているガス及び複数成分で構成されているガスのいずれをも意味しており、本実施の形態例におけるガスａ、ｂについては、単成分又は複数成分のどちらであっても良い。さらに、ガスｂが複数成分で構成されている場合で、ガスｂ中の成分変化により、ガスｂからガスｂ’へ変化する場合には、そのガスｂを構成している単成分毎に、ガスａとの合流後のガスａ＋ガスｂ’中の当該単成分毎に濃度の変化を本物性変化推定方法にて推定し、かかるのちにその成分毎の濃度から物性を求めることで、ガスｂからガスｂ’への変化に対しての、合流後のガス導管内のガス物性変化を推定できる。
【００３６】
図４は、本実施の形態例に係る物性変化推定装置の構成図である。図４に示す物性変化推定装置１０が、本実施の形態例に係る物性変化推定装置であり、具体的には、パーソナルコンピューターなどのコンピュータシステムで構成することができる。入力部１１は、前述した混合比の変更後の挙動を推定するために必要な各種推定条件の入力を行なうための部分である。
【００３７】
例えば、混合させるガスａ、ｂ、初期状態での界面の形状、すなわち、ガス導管の各位置ｘにおける濃度Ｃ（ｘ）、及び、ｘ＝０におけるガス流速ｕ及びガス圧力Ｐなどが当該入力部１１から物性変化推定装置１０に入力される。なお、この入力部１１は、当該装置の操作者が入力操作を行なうためのキーボード、マウス、ディスプレイ等、及び前記入力操作を促すための画面表示や入力されたデータの取得の処理を記述したプログラム、当該プログラムに従って処理を実行する制御装置（ＣＰＵ）等で構成される。
【００３８】
次に、推定処理部１２は、前記入力部１１から入力された推定条件等に基づいて、混合比の変更後のガス導管１内におけるガス濃度変化を推定する部分である。推定処理の具体的な内容については後述するが、当該推定処理部１２が実行する数値解析の手法に大きな特徴を有している。なお、当該推定処理部１２は、処理内容を記述したプログラム、当該プログラムに従って処理を実行する制御装置（ＣＰＵ）等で構成される。
【００３９】
データ格納部１３は、前記推定処理部１２が処理を実行するために必要な各種データを格納するデータ記憶手段であり、推定処理の度に操作者によって入力部１１から入力しなくてはならないデータ以外のデータを格納する。例えば、本物性変化推定装置１０が処理の対象とするガス導管網１の各所における導管の寸法（内径、長さ等）、輸送する各ガスの物性値（動粘度係数又は静粘度係数、等）などが予め格納される。なお、当該データ格納部１３は、ＨＤＤ等で構成される。
【００４０】
また、出力部１４は、前記推定処理部１２が行なった推定処理の結果を出力する部分であり、ガス導管１の所定の箇所又は所定の時間におけるガス濃度変化の推定結果を、表形式やグラフ形式で出力する。なお、この出力部１４は、ディスプレイ、プリンタへの出力インターフェース等、及び前記ディスプレイへの画面表示や前記プリンタへの出力の処理を記述したプログラム、当該プログラムに従って処理を実行する制御装置（ＣＰＵ）等で構成され、前記推定結果は、画面表示や紙等への印刷という形で出力される。
【００４１】
本物性変化推定装置１０は、以上説明したような構成を有し、以下に示すような内容で処理を実行する。図５は、本物性変化推定装置１０による処理の流れを例示したフローチャートである。
【００４２】
当該装置を用いてガス導管１内のガス混合比を変更させた場合のその後の挙動を推定しようとする操作者は、まず、前述した入力部１１を用いて推定処理に必要な各種の条件を入力する（ステップＳ１）。なお、前述のとおり、本物性変化推定装置１０では、図２に示した状態を前提とした推定処理を行なうものとする。
【００４３】
この推定条件の入力では、ガス導管１で輸送されるガスａ、ｂの種類、混合比変更前後のガスｂの濃度（α、β）及び初期状態における界面の形状（Ｃ（ｘ））、初期位置（ｘ＝０）における混合ガスのガス流速ｕ及びガス圧力Ｐ、ガス導管１の位置（ガス導管１を特定する情報）、ガス導管１内を流れるガスの温度等が入力され、ガスの濃度変化挙動を把握したい特定の位置（例えば、図２の消費点等）がある場合にはその位置を特定する情報を入力する。なお、データ格納部１３に、ガス導管網１の各所における導管の寸法及び物性値、各ガスの物性値などが格納されていない場合には、これらのうち必要なデータについては操作者によって入力される。また、本推定処理における数値計算では、後述するように、ガス導管１をメッシュに分割してメッシュ毎に逐次計算を行っていくが、当該メッシュのサイズを変更可能とする場合には、操作者がメッシュサイズも合わせて入力する。
【００４４】
次に、推定処理部１２は、当該推定処理を行うために必要な各種データを前記データ格納部１３から取得する（ステップＳ２）。具体的には、推定対象とするガス導管１について、その長さＬ、各部の内径Ｄ、及び内壁の摩擦係数λ、また、ガスａ及びガスｂの気体分子量Ｍａ、Ｍｂ、普遍気体定数Ｒ０、ガスの静粘度係数μ等のデータを取得する。なお、ガス導管１の各部（各メッシュ）の内径Ｄ及び摩擦係数λは、下記推定計算の過程で必要な際に取得するようにしてもよい。
【００４５】
次に、推定処理部１２は、以上入力されたデータと取得したデータを用いて、推定計算を実行する（ステップＳ３）。まず、当該推定計算で用いる式について説明する。当該推定計算では、以下に示す７つの式を用いて、各時間ｔ及び各位置ｘにおけるガス濃度Ｃを求める。
【００４６】
【数１】
【００４７】
【数２】
【００４８】
【数３】
【００４９】
なお、上式においてＫは、実験から以下のように与えられる。
Ｋ＝Ａ３ｕＤＲｅ−０．１２５ （４−１）
（分子量の式）
Ｍ＝Ｍａ＋（Ｍｂ−Ｍａ）Ｃ （５）
（気体定数の式）
Ｒ＝Ｒ０／Ｍ （６）
（モル密度と密度の関係式）
Ｘ＝ρ／Ｍ （７）
なお、上記（１）−（７）式において、各記号が表わすものは以下の通りである。
ｔ：時間 ［ｓ］
ｘ：流れ方向の座標（位置） ［ｍ］
ｕ：ガス流速 ［ｍ／ｓ］
Ｐ：ガス圧力 ［Ｐａ］
ρ：ガスの密度 ［ｋｇ／ｍ３］
λ：ガス導管内壁の摩擦係数 ［‐］
Ｄ：ガス導管の内径 ［ｍ］
Ｘ：ガスｂのモル密度 ［ｍｏｌ／ｍ３］
Ｒ：気体定数 ［Ｊ／（ｋｇ・Ｋ）］
θ：ガスの温度 ［Ｋ］
Ｃ：ガスｂの濃度 ［‐］
Ｋ：軸方向（流れ方向）の拡散係数 ［ｍ２／ｓ］
Ｒｅ：レイノズル数 ［‐］
Ｍ：ガス分子量 ［ｋｇ／ｍｏl］
Ｍａ：ガスａの分子量 ［ｋｇ／ｍｏl］
Ｍｂ：ガスｂの分子量 ［ｋｇ／ｍｏl］
Ｒ０：普遍気体定数 ８．１３４［Ｊ／（ｋｇ・Ｋ）］
Ａ３：比例定数（例えば、２．５９）
なお、レイノズル数Ｒｅは、上述したガスの静粘度係数μ、ガス導管の内径Ｄ、及びガスの密度ρから求めることができる。
【００５０】
以上の式において、時間ｔ及び位置ｘに対して未知数である変数は、Ｃ、Ｍ、Ｐ、ｕ、ρ、Ｒ、Ｘの７つであり、推定処理部１２は、この７つの未知数を上記７つの式（（１）−（７）式）を解くことによって求める。すなわち、推定計算を実行する。
【００５１】
なお、上記（２）式において、一般的な質量保存側に係る式を用いずに、モル密度保存式を用いたことが本装置の一つの特徴である。一般的な質量保存側に係る式を用いた場合、ガス密度が用いられることになるが、この密度の変化の要因としては、ガス圧力の変化に起因するものと、ガス成分（濃度）の変化に起因するものの二つが存在する。ここで対象としている混合ガスは、上述の通り、界面付近で濃度が大きく変化するので、上記ガス成分（濃度）の変化に起因する密度の変化が大きく、一般的な質量保存側に係る式を用いた場合、数値解析で用いる差分法において発散が生じてしまう可能性が高い。そこで、本装置では、このような発散を防止するため、モル密度保存式を用いている。
【００５２】
次に、推定処理部１２が行う上記７つの式を用いた推定計算の具体的な内容について説明する。当該推定計算はいわゆる数値解析によって行うが、その際に用いるメッシュ（解析対象を分割したもの）について、まず、説明する。
【００５３】
図６は、本推定計算で用いるメッシュを説明するための図である。図６において、ガス導管１の長さＬの部分（ｘ＝０からｘ＝Ｌの部分）が推定対象を表している。そして、図中の矢印がガスの流れ方向を示しており、従って、推定対象において左側が上流となる。
【００５４】
当該推定対象をｎ個に等分割した各領域（第一分割領域）が、本実施の形態例における１種類目のメッシュ（以下、粗メッシュと呼ぶ）であり、図中のＥ０〜Ｅｎ−１がそれに相当する。かかる粗メッシュは、ガス流速ｕ及びガス圧力Ｐを求める際に用いるものであり、例えば、推定対象の長さＬが１０ｋｍ程度である場合に、これに限定されるものではないが、一つの粗メッシュの長さを１００ｍ程度とすることができる。図６に示す例では、ある時間ｔにおいて、最上流の粗メッシュＥ０に対して、その左端のガス流速ｕ０及びガス圧力Ｐ０を用いて計算がなされ、Ｅ０の右端のガス流速ｕ１及びガス圧力Ｐ１が求められる。以下、同様に下流に向かって各粗メッシュに対して計算がなされる。
【００５５】
次に、上述した一つの粗メッシュをさらにｋ個に等分割した各領域（第二分割領域）が、本実施の形態例における２種類目のメッシュ（以下、細メッシュと呼ぶ）であり、図中のｅｉ，０〜ｅｉ，ｋ−１がそれに相当する。かかる細メッシュは、ガスｂの濃度Ｃを求める際に、すなわち界面形状を求める際に、用いるものであり、例えば、一つの粗メッシュの長さが１００ｍ程度である場合に、これに限定されるものではないが、一つの細メッシュの長さを１〜１０ｍ程度とすることができる。推定計算では、ある時間ｔにおいて、各細メッシュについて、それぞれ、ガスｂの濃度Ｃ（図中のＣｉ，０〜Ｃｉ，ｋ−１）が求められることになる。なお、後述するが、この細メッシュを用いた計算は全ての粗メッシュについて行う必要はなく、界面が存在する粗メッシュ及びその近傍について行われる。また、当該細メッシュ毎に、前述したガス導管１の各仕様が格納されていることが望ましい。
【００５６】
このように、本推定方法では、２種類のメッシュを使用して計算を行う。これは、気体の特性に基づくものであり、ガス成分が変化する前記界面付近においても、ガスの流速ｕ及び圧力Ｐは大きな変化をしない一方、ガス濃度Ｃは急激に変化するので、前者については、比較的大きなメッシュでの計算でも精度上問題がなく、後者については、精度良く状態を把握するために細かいメッシュでの計算が必要となる。
【００５７】
以下、具体的な計算手順について説明する。図５に戻って、まず、初期時におけるガス濃度Ｃを決定する。（ステップＳ３−１）。図７は、初期の界面形状Ｃ（ｘ）について説明するための図である。通常の推定においては、ガス混合地点からのガス混合時点以降を対象とするので、当該初期時においては、通常、界面が粗メッシュＥ０に存在する。従って、上記界面が存在する箇所でのガス濃度Ｃの決定は、図７に例示するように、入力され界面の形状Ｃ（ｘ）（図中の曲線）に従って、各細メッシュのＣ値を決定する。図７の例では、ｅ０，０、ｅ０，１等の細メッシュについて、対応するＣ０，０、Ｃ０，１等の値が、図中のプロット点×のように決定される。
【００５８】
そして、上記界面が存在する箇所については、上記各細メッシュのＣの値から、対応する粗メッシュの濃度Ｃを決定する。具体的には、それらの平均値から粗メッシュの値を求める。また、界面が存在しない箇所については、界面の下流側ではＣの値がαであり、その上流側ではβであるので、粗メッシュ単位でそのように決定する。
【００５９】
次に、時間変数ｔをｔ０に設定後（ステップＳ３−２）、初期位置（ｘ＝０）における圧力Ｐ０及び流速ｕ０を、前記入力値に従って設定する（ステップＳ３−３）。そして、この時点におけるガス流速ｕ及びガス圧力Ｐを算出する（ステップＳ３−４）。ここで、前述した（５）式において、上記決定された濃度ＣによりＣの値が既知となるので、Ｍの値が決定され、そのＭの値により、上記（６）式からＲの値が求められる。この時点で未知数は、ｕ、Ｐ、Ｘ、ρの４つであり、それらの関係式が上記（１）−（３）及び（７）式として４つ与えられるので、これらの関係式を用いて上記４つの未知数の値が決定される。
【００６０】
具体的には、スタガート格子差分法などを用いて各粗メッシュ毎に数値計算を行って求める。図６に示した例では、前述のようにｕ０、Ｐ０の値が与えられているので、Ｅ０について数値計算を行いｕ１、Ｐ１の値等を算出する。以下、同様に、下流方向に計算を行って、ｕ２、Ｐ２〜ｕｎ、Ｐｎの値を求める。
【００６１】
次に、時間変数ｔをｔ＋Δｔとし（ステップＳ３−５）、次の時間（Δｔ経過後）におけるガス濃度界面の中心位置を計算する（ステップＳ３−６）。ガス濃度界面の中心位置ｘ界面は、前時刻における界面の中心位置ｘ前界面と界面が存在した粗メッシュの流速ｕ界面が存在する粗メッシュから、ｘ界面＝ｘ前界面＋Δｔ×ｕ界面が存在する粗メッシュとして計算される。
【００６２】
その結果から、界面が推定対象の最下流の位置を通過したかを判断し（ステップＳ３−７）、通過していない場合には（ステップＳ３−７のＮｏ）、この時刻におけるガス濃度Ｃと圧力Ｐ及び流速ｕの算出を行う。
【００６３】
具体的には、まず、上記界面の中心位置近傍の細メッシュについて濃度Ｃの値を算出する。すなわち、界面形状を求める（ステップＳ３−８）。ここでは、前時刻におけるガス流速ｕが求められているので、その値を前述した（４）及び（４−１）式に代入してＣの値を求めることができる。ここでも、具体的にはスタガート格子差分法などを用いて各細メッシュ毎に数値計算を行って求める。その後、上記ステップＳ３−１と同様にして、粗メッシュ毎のガス濃度Ｃを決定する（ステップＳ３−９）。
【００６４】
次に、処理が前記ステップＳ３−４に戻って、上述の通り、この時刻における圧力Ｐ及び流速ｕが求められる。
【００６５】
以降、上記ステップＳ３−７において、界面が推定対象の最下流の位置を通過していると判断されるまで、同様の処理が繰り返される。そして、界面が通過していると判断された場合には（ステップＳ３−７のＹｅｓ）、推定計算（Ｓ３）を終了し、推定結果の出力を行う（ステップＳ４）。
【００６６】
このように、本推定計算では、各時間（各時刻）において、前時刻で求められた粗メッシュ毎のガス流速ｕを用いて界面近傍について細メッシュ毎にガス濃度Ｃを求め（すなわち、界面形状を求め）、その求められた濃度Ｃに基づいて粗メッシュ毎に流速ｕを求める。そして、当該処理を時間経過に従って順次繰り返して実行することにより、任意の時刻ｔ、任意の位置ｘでのガス濃度Ｃを求めることができる。なお、上記説明では、各時刻において先に濃度Ｃを求めその後に流速ｕを求めたが、前時刻における濃度Ｃを用いて先に流速ｕを求め、その後に濃度Ｃを求めるようにしてもよい。
【００６７】
推定結果の出力では、用途に応じて様々な形式での出力が可能であるが、例えば、ガス消費設備３の位置での、界面の通過時における濃度Ｃの時間変化をグラフ形式で出力する、ことが可能である。
【００６８】
図８は、本物性変化推定装置１０による推定計算の結果の一例を示した図である。図８に示す例は、一つの粗メッシュにおける各時間での界面（におけるＣの値）をプロットしたものである。実線は界面が当該メッシュに入ってから１０秒後の状態を示しており、点線は、２０秒後の状態を示している。１０秒後の状態と２０秒後の状態を比較すると、後者では界面の傾斜がなだらかになり、ｘ方向に拡散現象が起こっていることがわかる。また、３０秒後では、当該領域に界面がなくなり、界面が通過していることがわかる。
【００６９】
なお、上述の説明では、初期位置（ｘ＝０）における混合ガスのガス流速ｕ及びガス圧力Ｐを推定条件として入力したが、これらの初期値を時間ｔによって変化する変数ｕ（ｔ）及びＰ（ｔ）として与えるようにしてもよい。この場合には、各時間ｔにおけるｕ０及びＰ０をｕ（ｔ）及びＰ（ｔ）によって得られる値とすることにより、上述の場合と同様に計算を行うことができる。従って、この場合には、図５に示した処理手順において、ステップＳ３−９の後に、処理がステップＳ３−３に移行することになる。
【００７０】
また、図２等においては、ｘ＝０を混合発生点、ｘ＝Ｌを消費点として、これらの間を推定対象としたが、これに限定されることはなく、本推定計算は任意の領域に適用することができる。
【００７１】
なお、本物性変化推定装置１０は、スタンドアローンのシステムとしたが、ネットワークを介して、操作者による入力や結果の出力が行なわれるようにしてもよい。また、ガス導管網１の各所で測定されるガス流速などのリアルタイムデータを推定の入力値として用いるべく、それら測定装置等と通信可能にして用いてもよい。
【００７２】
以上説明したように、本実施の形態例に係る物性変化推定装置１０では、推定対象範囲をメッシュに分割して数値計算するのでメッシュ毎にガス導管１の仕様を設定可能であり、また、初期時の界面形状をＣ（ｘ）で与えられるので複雑な界面形状の場合にも対応可能である。さらに、上述のように、初期位置のガス流速ｕ及びガス圧力Ｐを時間ｔの関数として与えることもできる。従って、本装置１０により、従来の装置よりも適用範囲が広くなって、実際に存在する様々なケースのシミュレーションを適確に行うことが可能となる。
【００７３】
また、上記（４−１）式ではガス導管１内の流れが乱流である場合の式を採用しており、流れが層流である場合に限らず適用することができる。
【００７４】
また、前述の通り、連続性を担保する関係式としてモル密度保存式を用い、数値計算上の不具合を防止している。
【００７５】
また、本物性変化推定装置１０では、推定対象を粗メッシュと細メッシュに分割し、気体の特性に基づいて、細メッシュ毎の計算が必要な場合のみ細メッシュでの計算を実行する。従って、ガス導管網のように推定対象が長い場合に解析上の負荷を大幅に軽減することが可能である。
【００７６】
このように、本物性変化推定装置１０によって、ガスの物性が変化した場合の挙動を、適確にかつ低コストで把握することが可能となり、ガス供給の信頼性をより向上させることができる。
【００７７】
本発明の保護範囲は、上記の実施の形態に限定されず、特許請求の範囲に記載された発明とその均等物に及ぶものである。
【符号の説明】
【００７８】
１ ガス導管（網）
２ ガス供給元
３ ガス消費設備
１０ 物性変化推定装置
１１ 入力部
１２ 推定処理部
１３ データ格納部
１４ 出力部
【特許請求の範囲】
【請求項１】
ガス導管内を流れる互いに種類の異なるガスの混合比を変更させた場合における、前記ガス導管内を流れるガスの物性変化挙動を推定する処理をコンピュータに実行させる物性変化推定プログラムであって、
前記推定開始後の各時刻において、
前記推定処理の対象範囲を所定長さで分割した第一分割領域毎に、前記混合比に基づいて、前記ガスの流速を求める流速決定ステップと、
前記第一分割領域を更に分割した第二分割領域毎に、前記流速に基づいて前記混合比を求める混合比決定ステップとを、
前記コンピュータに実行させる
ことを特徴とする物性変化推定プログラム。
【請求項２】
請求項１において、
前記混合比決定ステップは、前記混合比が変化する部分である界面が存在する領域について実行される
ことを特徴とする物性変化推定プログラム。
【請求項３】
請求項１あるいは請求項２において、更に、
少なくとも、前記第二分割領域毎の前記ガス導管の仕様と、前記推定開始時における前記ガス導管の各位置における前記混合比と、前記ガス導管の最上流における前記ガスの流速及び圧力を取得するステップを、前記コンピュータに実行させ、
前記流速決定ステップと前記混合比決定ステップにおいて、当該取得された情報が用いられる
ことを特徴とする物性変化推定プログラム。
【請求項４】
請求項３において、
前記ガス導管の最上流におけるガスの流速及び圧力は、時間を変数とする関数として取得される
ことを特徴とする物性変化推定プログラム。
【請求項５】
請求項１乃至請求項４のいずれかにおいて、
前記流速決定ステップは、運動方程式に基づいて実行され、
前記混合比決定ステップは、移流拡散方程式に基づいて実行される
ことを特徴とする物性変化推定プログラム。
【請求項６】
請求項１乃至請求項５のいずれかにおいて、
前記流速決定ステップは、モル密度保存式に基づいて実行される
ことを特徴とする物性変化推定プログラム。
【請求項７】
ガス導管内を流れる互いに種類の異なるガスの混合比を変更させた場合における、前記ガス導管内を流れるガスの物性変化挙動を推定する物性変化推定装置であって、
前記推定のための条件を受け付ける入力部と、
前記推定のためのデータを格納する格納部と、
前記入力部が受け付ける条件と前記格納部が格納するデータに基づいて、前記推定開始後の各時刻において、
前記推定の対象範囲を所定長さで分割した第一分割領域毎に、前記混合比に基づいて、前記ガスの流速を求め、
前記第一分割領域を更に分割した第二分割領域毎に、前記流速に基づいて前記混合比を求める
推定処理部と、を有する
ことを特徴とする物性変化推定装置。
【請求項８】
ガス導管内を流れる互いに種類の異なるガスの混合比を変更させた場合における、前記ガス導管内を流れるガスの物性変化挙動を推定する物性変化推定装置における物性変化推定方法であって、
前記物性変化推定装置は、前記推定のための条件を受け付ける入力部と、前記推定のためのデータを格納する格納部と、推定処理部とを備え、
前記推定開始後の各時刻において、
前記入力部が受け付ける条件と前記格納部が格納するデータに基づいて、
前記推定処理部が、前記推定の対象範囲を所定長さで分割した第一分割領域毎に、前時刻で求めた前記混合比に基づいて、前記ガスの流速を求める工程と、
前記推定処理部が、前記第一分割領域を更に分割した第二分割領域毎に、前記求めた流速に基づいて前記混合比を求める工程を有する
ことを特徴とする物性変化推定方法。
【請求項１】
ガス導管内を流れる互いに種類の異なるガスの混合比を変更させた場合における、前記ガス導管内を流れるガスの物性変化挙動を推定する処理をコンピュータに実行させる物性変化推定プログラムであって、
前記推定開始後の各時刻において、
前記推定処理の対象範囲を所定長さで分割した第一分割領域毎に、前記混合比に基づいて、前記ガスの流速を求める流速決定ステップと、
前記第一分割領域を更に分割した第二分割領域毎に、前記流速に基づいて前記混合比を求める混合比決定ステップとを、
前記コンピュータに実行させる
ことを特徴とする物性変化推定プログラム。
【請求項２】
請求項１において、
前記混合比決定ステップは、前記混合比が変化する部分である界面が存在する領域について実行される
ことを特徴とする物性変化推定プログラム。
【請求項３】
請求項１あるいは請求項２において、更に、
少なくとも、前記第二分割領域毎の前記ガス導管の仕様と、前記推定開始時における前記ガス導管の各位置における前記混合比と、前記ガス導管の最上流における前記ガスの流速及び圧力を取得するステップを、前記コンピュータに実行させ、
前記流速決定ステップと前記混合比決定ステップにおいて、当該取得された情報が用いられる
ことを特徴とする物性変化推定プログラム。
【請求項４】
請求項３において、
前記ガス導管の最上流におけるガスの流速及び圧力は、時間を変数とする関数として取得される
ことを特徴とする物性変化推定プログラム。
【請求項５】
請求項１乃至請求項４のいずれかにおいて、
前記流速決定ステップは、運動方程式に基づいて実行され、
前記混合比決定ステップは、移流拡散方程式に基づいて実行される
ことを特徴とする物性変化推定プログラム。
【請求項６】
請求項１乃至請求項５のいずれかにおいて、
前記流速決定ステップは、モル密度保存式に基づいて実行される
ことを特徴とする物性変化推定プログラム。
【請求項７】
ガス導管内を流れる互いに種類の異なるガスの混合比を変更させた場合における、前記ガス導管内を流れるガスの物性変化挙動を推定する物性変化推定装置であって、
前記推定のための条件を受け付ける入力部と、
前記推定のためのデータを格納する格納部と、
前記入力部が受け付ける条件と前記格納部が格納するデータに基づいて、前記推定開始後の各時刻において、
前記推定の対象範囲を所定長さで分割した第一分割領域毎に、前記混合比に基づいて、前記ガスの流速を求め、
前記第一分割領域を更に分割した第二分割領域毎に、前記流速に基づいて前記混合比を求める
推定処理部と、を有する
ことを特徴とする物性変化推定装置。
【請求項８】
ガス導管内を流れる互いに種類の異なるガスの混合比を変更させた場合における、前記ガス導管内を流れるガスの物性変化挙動を推定する物性変化推定装置における物性変化推定方法であって、
前記物性変化推定装置は、前記推定のための条件を受け付ける入力部と、前記推定のためのデータを格納する格納部と、推定処理部とを備え、
前記推定開始後の各時刻において、
前記入力部が受け付ける条件と前記格納部が格納するデータに基づいて、
前記推定処理部が、前記推定の対象範囲を所定長さで分割した第一分割領域毎に、前時刻で求めた前記混合比に基づいて、前記ガスの流速を求める工程と、
前記推定処理部が、前記第一分割領域を更に分割した第二分割領域毎に、前記求めた流速に基づいて前記混合比を求める工程を有する
ことを特徴とする物性変化推定方法。
【図１】
【図２】
【図３】
【図４】
【図５】
【図６】
【図７】
【図８】
【図２】
【図３】
【図４】
【図５】
【図６】
【図７】
【図８】
【公開番号】特開２０１０−２７６０５３（Ｐ２０１０−２７６０５３Ａ）
【公開日】平成２２年１２月９日（２０１０．１２．９）
【国際特許分類】
【出願番号】特願２００９−１２６７８７（Ｐ２００９−１２６７８７）
【出願日】平成２１年５月２６日（２００９．５．２６）
【出願人】（０００２２０２６２）東京瓦斯株式会社 (1,166)
【Ｆターム（参考）】
【公開日】平成２２年１２月９日（２０１０．１２．９）
【国際特許分類】
【出願日】平成２１年５月２６日（２００９．５．２６）
【出願人】（０００２２０２６２）東京瓦斯株式会社 (1,166)
【Ｆターム（参考）】
[ Back to top ]