JP6167523B2

JP6167523B2 - 計算装置

Info

Publication number: JP6167523B2
Application number: JP2013002500A
Authority: JP
Inventors: 内海　雅彦; 雅彦内海
Original assignee: IHI Corp
Current assignee: IHI Corp
Priority date: 2013-01-10
Filing date: 2013-01-10
Publication date: 2017-07-26
Anticipated expiration: 2033-01-10
Also published as: JP2014134957A

Description

本発明は、管群に直交して流れる流体により発生する管群の振動を予測する管群振動予測方法およびこのような振動を計算する計算装置に関する。

従来、加圧水型原子炉（ＰＷＲ）に備えられる蒸気発生器（ＳＧ）においては、原子炉から供給される一次冷却水を流す伝熱管を並列配置して管群とし、この管群を伝熱面としてこの管群の軸方向と直交する向きに二次冷却水を流して熱交換している。このような蒸気発生器においては、熱交換の効率化を図るため、管群に直交した流れの流速を高めることが求められている。

このような管群において、この直交流の流速がある値に達すると、不安定振動が発生することが知られている（下記非特許文献１参照）。原子炉の安全性を確保するためにはこのような不安定振動を制御することが必要であり、不安定振動を引き起こす臨界流速を予測する研究がなされてきた。

この臨界流速の予測には、ＣＦＤ（数値流体力学）による数値シミュレーション（下記非特許文献２参照）が提供されている。また、１本の管の振動によって隣接する複数の管に作用する力を実験により求め、この実験データに基づいて不安定振動解析モデルを作成する方法が提供されている（下記非特許文献３参照）。

日本機械学会編、事例に学ぶ流体関連振動、ｐ．６６市岡丈彦、川田裕、中村友道、泉元、小林敏雄、高松洋、「流動解析を用いた管群の流力弾性振動に関する研究（第１報、並列２円柱および１列管群の解析）」、日本機械学会論文集（Ｂ編）、１９９５年、６１巻、５８２号、ｐ．１４５−１５１田中博喜、「円柱群の流力弾性振動に関する研究（流水中の円柱群）」、昭和５６年、日本機械学会論文集（Ｂ編）、４７巻、４１９号、ｐ．１２２４−１２３３

しかしながら、ＣＦＤによるシミュレーションは、モデルの入力、計算処理等に多くの時間とコストを要し、多数のパラメータについて臨界流速を予測することは困難であった。また、実験データに基づいて不安定振動解析モデルを作成する方法についても、多数のパラメータについて臨界流速を予測することは困難であった。

この出願に係る発明は、上述の実情に鑑みて提案されるものであって、多くの時間とコストを要することなく、多数のパラメータについて臨界流速を予測することができるような管群振動予測方法及び計算装置を提供することを目的とする。

この発明に係る管群振動予測方法は、所定径の管が所定間隔にて周期的に配置されてなる管群において、この管群の軸方向に直交して流れる流体により当該管群に発生する振動を予測する管群振動予測方法であって、前記管群の軸方向に直交する面内において、前記管群を含む領域について各管をそれぞれ含む正多角形のセルに分割するステップと、前記正多角形のセルにおいて、この正多角形のセルの外周と、この正多角形のセルに含まれる管との境界をなす内周とについて、未知量をモード展開するステップと、前記正多角形のセルにおいて、管の中心を始点として前記外周と前記内周を結ぶ動径上で未知量をモード展開するステップとを有するものである。

また、この発明に係る管群振動予測方法は、所定径の管が所定間隔にて周期的に配置されてなる管群において、この管群の軸方向に直交して流れる流体により当該管群に発生する振動を予測する管群振動予測方法であって、前記管群の軸方向に直交する面内において、前記管群を含むセルについて各管をそれぞれ含む正多角形のセルに分割するステップと、前記正多角形のセルにおいて、この正多角形のセルの外周と、この正多角形のセルに含まれる管との境界をなす内周とについて、未知量をモード展開するステップと、前記正多角形のセルにおいて、管の中心を始点として前記外周と前記内周を結ぶ動径上で未知量をモード展開するステップと、前記モード展開した未知量についてガレルキン法を適用するステップとを有するものである。

前記未知量は前記流体の流速及び圧力であり、これら流速及び圧力についてナビエ・ストークスの式及び連続条件に基づいて、前記流速及び圧力の一般化座標の時間に関する常微分方程式を導き、この常微分方程式に基づいて前記流体の定常流による流れ場を決定することが好ましい。前記未知量として前記流速及び圧力にそれぞれ振動外乱を重畳し、これら流速、圧力及び振動外乱についてナビエ・ストークスの式、連続条件及び前記定常流による流れ場に基づいて、前記振動外乱の一般化座標の時間に関する常微分方程式を導き、この常微分方程式に基づいて前記管群の振動を予測することが好ましい。前記振動外乱の一般化座標の時間に関する常微分方程式に基づいて、前記管群の不安定振動が発生する前記流体の臨界流速を決定することが好ましい。

前記管は円柱形状であり、前記管は４回回転対称性を有するように配置され、前記正多角形は正方形であることが好ましい。前記管群は加圧水型原子炉の蒸気発生器における伝熱管群を構成し、前記流体は前記蒸気発生器における二次冷却水であることが好ましい。

この発明に係る計算装置は、前記管群振動予測方法を適用したものであって、前記流速及び圧力の一般化座標の時間に関する常微分方程式と前記振動外乱の一般化座標の時間に関する常微分方程式の表現を含む所定のプログラムを格納した記憶手段と、前記管の径、前記管を配置した間隔、前記管の単位長さ当たりの質量、前記管の固有振動数、前記流体の密度、前記流体の粘性係数の少なくとも１つの値を入力する入力手段と、前記記憶手段に格納された前記プログラムを読み込んで実行し、前記入力手段に入力された前記少なくとも１つの値に基づいて前記臨界流速を計算し、この値が閾値を越えたかどうかを判定する演算手段と、前記演算手段で計算した前記臨界流速と前記判定の結果を出力する出力手段とを含むものである。

この発明は、管を含む正多角形のセルについてのモード展開により解析的に得られた常微分方程式を使用しているので、多くの時間とコストを要することなく、多数のパラメータについて臨界流速を予測することができる。

流れを受ける管群を示す図である。１セルにおける座標設定を示す図である。流体密度の大きい領域における換算流速の臨界値の質量減衰パラメータ依存性を示す図である。レイノルズ数２．２９×１０^４における流れ場を示す図である。流体密度の小さい領域における換算流速の臨界値の質量減衰パラメータ依存性を示す図である。式（２４）における粘性力の導出を示す図である。管群振動を予測する計算装置の構成を示す図である。計算装置における一連の動作の流れを示す図である。

以下、本発明に係る管群振動予測方法及び計算装置について、図面を参照して詳細に説明する。

〔１．まえがき〕
従来、加圧水型原子炉（ＰＷＲ）の蒸気発生器において、管群が管の軸と直交する方向の流れを受け、その流速がある値に達すると不安定振動が発生することがよく知られている（非特許文献１）。このため、不安定振動を引き起こす流速（臨界流速）を予測することは、原子力機器の安全性を確保するために、非常に重要な問題となっている。

従来、臨界流速の予測法として、ＣＦＤ数値シミュレーション（非特許文献２）、１本の管の振動によって隣接する複数の管に作用する力を実験により求め、この実験データを下に不安定振動解析モデルを作成する方法（非特許文献３）が用いられている。しかし、管群に対する流体支配方程式系を解析的手法で解く方法は、研究されていない。本実施の形態は、このような方法を提案したものである。

〔２．方法の概要〕
ＣＦＤ解析では、百万オーダの膨大自由度を要し、膨大な計算時間、コストがかかる。従って、設計のパラメータスタディのために多くのパラメータについて流速を少しずつ上げて臨界流速を決定することは、実際上困難となる。また、隣接管に作用する力係数を実験で求めることによって解析モデルを作成する方法においても、多数のパラメータスタディを行う上での不便が生じる。

このような問題を解決するため本実施の形態は、解析的手法による高速計算法を開発することを目的とする。本実施の形態では、有限要素法など従来の多くの数値解法が領域全体を小要素に多数分割するのと対照的に、比較的大きい領域（各管を含む多角形の領域）に分割し、それらの境界での未知量をモード展開する。そして、これらのモード座標で各領域内での未知量を表わすことによって自由度を低減する、境界モード展開法を提案する。このような準解析的方法による自由度低減、計算効率化の可能性は、従来看過されているようである。提案する高速解析法により、従来データに近い臨界流速が高速算定できることを確かめる。

〔３．解析方法〕
図１のように流れを受ける管群１０について、不安定振動が発生する臨界流速を求める。この管群１０は、それぞれ半径ａの第１の管１１_１から第９の管１１_９までの９本の管１１_１〜１１_９が、軸間の間隔（ピッチ）をｐとして、Ｘ方向及びＹ方向にそれぞれ３列に、各管１１_１〜１１_９の軸が延びるｚ方向について４回回転対称になるように配置されている。以下では簡単のため、各管１１_１〜１１_９を特定する添字１〜９を省略する。

管１１の軸方向に流速、圧力、管の変位の変化しない２次元問題を考え、この軸に直交する２次元平面を想定する。この面内において、最初に各管１１の静止条件下での定常（平均）流れ場を決定し（３．１節）、次に振動外乱を重畳させて安定性解析を行う（３．２節）。

〔３．１流れ場の決定〕
〔３．１．１解の表示〕
各管１１を含む正方形領域であるセル１２の頂点の番号を定義し、頂点ｉでのｘ，ｙ方向の流速、圧力をそれぞれ、

とする。

図２のような１セル１２に着目する。各外周辺に沿って、ｘまたはｙが増す方向に座標ｓを設定し、各辺上でのｘ，ｙ方向の流速、圧力を、辺の両端ｉ，ｊでの値の直線補間項に精度向上のためのモード展開項を付加して、次のように表わす。

ここで

は、次のような直線補間関数およびモード関数

であり、

は、ｉ，ｊを結ぶ辺上で定義された一般化座標であり、ｃは辺の番号であり、

は、辺の両端の周方向角

であり、

とｓの関係は各外周辺で次式によって与えられる。

ここで

流体運動解析は、図２のような、各セル１２に関する極座標を用いて行う。このため、式（１）の流速成分を、次のように極座標に変換する。

ここで、方向余弦は

である。

周方向座標

での内周と外周を結ぶ動径上での流速、圧力を、次のように表す。

ここで

である。

式（６）右辺第１項は、内周ｒ＝ａ、外周ｒ＝ｒ_ｍａｘでそれぞれ０，１となり、外周で解（４）に等しくなるための項である。式（６）右辺第２項は、精度向上のための直交関数展開項で、

である。内周での流速は、境界層が非常に薄いため、法線方向成分のみが０となる境界条件を満たすように設定している。

ここで、数（６）の第１項は、外周ｒ＝ｒ_ｍａｘにおいて

となって外周での解（式（４））に接続、帰着させることができる。内周ｒ＝ａでは消える。

数（６）の第２項は、外周ｒ＝ｒ_ｍａｘでは消え、内周ｒ＝ａで半径方向流速が０、周方向流速と圧力は０でないという条件を満たしながら、動径上

における解をｓｉｎ，ｃｏｓの直交関数系でフーリエ展開したものである。ここで、動径座標は管１１の中心を始点とすることができる。

式（４）に式（１）を代入した結果を、式（６）右辺第１項に代入することにより、各セル１２内の任意位置での流速と圧力が、頂点でのｘ，ｙ方向の流速と圧力、式（１）右辺第３項のモード展開係数、式（６）右辺第２項のモード展開係数とで表される。すなわち、流体の無限自由度が、これらの時間変数に縮小する。そこで、このようにして得た式を自由度縮小式と称する。

〔３．１．２ガレルキン法〕
自由度縮小式を、下記の連続条件とナビエ・ストークス方程式の変分重み付き表示式に代入して、未知の時間変数に関する常微分方程式を導く。

ここで

である。

常微分方程式は、次のマトリックス方程式の形に導く。

ここで

である。

式（１１）を時間積分し、時間変動の小さい定常近似解に達した解によって流れ場を決定する。このとき数値計算の安定化を行った。多用されている人工粘性は流れ場を変えてしまうので使わず、運動方程式の時間微分を陽に含む非定常慣性項について係数ρ_ｆを

とする人工慣性を導入した。粘性項と異なり非定常慣性項は定常状態で０に近づくため、求めるべき定常近似解への影響は小さいことに留意した工夫である。

〔３．２振動解析〕
〔３．２．１解の表示〕
解を次のように設定する。

これは、前節で求めた定常成分

に振動外乱

が重畳した形である。

３．１．１節の自由度低減手法を、振動成分に関する方程式系の導出にも用いる。まず、セル１２の外周辺上でのｘ，ｙ方向の流速、圧力に関する振動成分を、辺の両端ｉ，ｊでの値

の直線補間と、モード展開項の和で表し、式（１）に対応して次式を得る。

流体運動解析は、図２のような、各セル１２に関する極座標を用いて行うため、流速成分を極座標に変換し、周方向座標

での内周と外周を結ぶ動径上での流速、圧力（振動成分）を、式（６）に対応して次のように表す。

このようにして、各セル１２内の任意位置での流速と圧力（振動成分）を、頂点でのｘ，ｙ方向の流速と圧力、式（１３）右辺第３項のモード展開係数、式（１４）右辺のモード展開係数で表す。すなわち、流体の無限自由度を、これらの時間変数に縮小する。

〔３．２．２ガレルキン法〕
ナビエ・ストークス方程式は、変分重み付きの形で次のように表される。

運動学的条件は、管１１表面のｘ，ｙ方向の速度ｑ_ｘ，ｑ_ｙの法線方向成分を管１１表面での湧出しと考え、湧出し項を流体連続条件式に導入することによって、次のようになる。

各管の運動方程式は、次式によって与えられる。

ここで

である。

式（１４）を式（１５），（１６），（１７），（１８）に代入することにより、未知の時間変数に関する常微分方程式系を下記のマトリックス方程式の形に導く。

ここで

は一般化座標によって構成される列ベクトルであり、マトリックス

の固有値を求めることによって安定性を判定する。

〔４．計算結果〕
〔４．１計算条件〕
次のパラメータについて計算を行った。

流れ場の決定では、粘性係数の代わりに乱流粘性係数をμの６０倍の値に設定した（ＣＦＤにおいては２０〜１００倍の値を用いることに基づく）。

管群１０の問題では、管１１が多数あるため、図１で、次の３個の領域

での流れ場は同じである。そこで、１個のセル１２のみについて、流れ場を決定する計算を行った。次に、Ｘ方向についても管１１が多数であるため周期対称な問題である。そこで、出入り口での湧出し、吸込みの影響の少ない次の領域

で求まった流れ場を９本の管１１の周りのセル１２に設定して安定性解析を行った。

〔４．２臨界流速の従来データとの比較〕
質量減衰パラメータと呼ばれる次の量（無次元）

を、減衰比ζ_ｂを変えることによって変え、換算流速（reduced velocity）と称されている下記の量（無次元）の臨界値を求めた。

ここで、Ｖ_ｃは管１１の間の隙間の流速で、流れ場を決めるときの式（１０）中の湧出し、吸込みソースの強さＶ_ｓｒｃのｐ／（ｐ−２ａ）倍である。臨界値を従来データ（非特許文献１）と比較した結果を図３に示す。図中の“□”は従来データであり、図中の折線は本解析により得られた値である。本解析により従来データに近い値が得られていることが確かめられる。

図４に流れ場の計算例を示す。この計算例においてレイノルズ数

である。

本解析のＶ_ｒ＝２．８はＶ_ｃ＝０．８４ｍ／ｓに相当する。図４でＸ＝０．０６ｍでの隣接管との隙間での最大流速は０．９ｍ／ｓであり、数値計算上の精度限界によりＶ_ｃ＝０．８４ｍ／ｓを超えているが、隙間に渡る平均ではＶ_ｃの０．８倍程度になる。従って、計算結果として生ずる隙間での平均的流速によって臨界流速を求めると、図４に示されている臨界流速の０．８倍となり、より保守的な評価となる。

上記は、流体密度の大きい液体の場合を対象とした、質量減衰パラメータの小さい範囲の結果である。しかし工学上は、蒸気発生器（ＳＧ）のように、ボイド率９０％、すなわち流体密度が１／１０に低下し、式（２０）の質量減衰パラメータがより大きい場合が重要になる。

図５は、この場合の計算結果を示す。図中の“□”は従来データ（非特許文献１）であり、図中の折線は本解析による値である。本解析では、減衰を一般的な値ζ_ｂ＝０．０４に固定し、管１１の単位長さ当りの質量を変えることによって、質量減衰パラメータを変えている。この場合も、本解析により、従来データに近い値が得られていることが確かめられる。従来データは、質量減衰パラメータが１を超えると、換算流速の臨界値の質量減衰パラメータに対する変化が不連続的に大きくなる。この傾向が、本解析においても現れている。

〔４．３自由度低減効果〕
本解析では、自由度数がＣＦＤに比べると少なくて済むことを確かめる。解の表示式（１３），（１４）より、流体解析に関する自由度数Ｎは次式によって求められる。

ここで、Ｎ_１は頂点の個数、Ｎ_２は辺の個数、Ｎ_１３は管の個数（セルの個数）であり、係数３は、ｘ，ｙ方向の流速成分と圧力とで３個の変数があることによる。図１の３×３配列の管群１０の場合について自由度数を求めると、Ｎ_１＝１６，Ｎ_２＝２４，Ｎ_３＝９であり、ｎ_ｍａｘ＝２，ｍ_ｍａｘ＝１，ｋ_ｍａｘ＝１で上のように従来データを説明し得る解析結果が得られる。

この場合の自由度数はＮ＝２１９である。ＣＦＤの数百万オーダに比べるとかなり少なくて済む。安定判別の固有値を求めるために要する計算時間は１ケースについて２分程度である（ただし臨界流速を求めるためには、固有値計算を１０回程度繰り返す必要がある）。このような低次元の固有値問題に帰着できることは、ＣＦＤ使用の際に流体と管１１の周波数差により流体解析に合わせて時間きざみを非常に小さくして時間積分する必要が生じる問題の回避に有効である。

〔５．計算装置〕
前述したような臨界流速の予測は、図７に示すような計算装置２０によって実現することができる。この計算装置２０は、ＣＰＵ、ＤＳＰの如き演算部２１、ＲＡＭ、ＲＯＭ、ハードディスクの如き記憶部２２、ＬＣＤ、プリンタの如き出力部２３、キーボード、マウスの如き入力部２４を含み、例えばパーソナルコンピュータを利用することができる。

図８に示す計算装置２０の一連の動作は、記憶部２２に格納された臨界流速算定プログラム２２ａを演算部２１が読み出して実行することにより実現される。この臨界流速算定プログラム２２ａは、前述のような手順によって得られた流速及び圧力の一般化座標の時間に関する常微分方程式と振動外乱の一般化座標の時間に関する常微分方程式の表現を含んでいる。

最初のステップＳ１においては、モデルを設定する。ここでは、前述したような管群１０のモデルを設定するものとする。入力部１４は、このモデルについて、管１１の径ａ、管１１を配置した間隔ｐ、管１１の単位長さ当たりの質量、管１１の固有振動数、流体の密度、流体の粘性係数の少なくとも１つの値を入力値として受け取る。演算部２１は、入力部２４が受け取った入力値を記憶部２２に格納する。

ステップＳ２においては、演算部２１は、記憶部２２に格納された入力値を読み出し、その数値計算部２１ａにおいて、この入力値に基づいて前記常微分方程式の表現を用いてこの管群における臨界流速について数値計算する。演算部２１は、得られた臨界流速の値を記憶部２２に格納する。

ステップＳ３においては、演算部２１は、記憶部２２に格納された臨界流速の値と、同じく記憶部２２に格納された所定の閾値２２ｂとを読み出す。演算部２１は、その判定部２１ｂにおいて、臨界流速値が閾値２２ｂを超えた場合にはＯＫと判定して一連のステップを終了する。一方、臨界流速が閾値を越えない場合にはＮＧとして判定して前のステップＳ１のモデル設定に手順を戻す。なお、閾値２２ｂは、入力部２４を介して設定することができる。

このような一連の工程において、臨界流速の値が所定の閾値を超えるまでモデル設定、数値計算、判定のループを繰り返すことにより、臨界流速が閾値を超えるようなモデル設定を可能としている。また、前述のような常微分方程式の表現を利用することにより、精度の高いモデル設定を可能としている。なお、これに限られず、臨界流速が閾値内に収まるようにモデル設定することもできる。

なお、このような臨界流速の算定は、記憶部２２に格納した臨界流速算定プログラム２２ａのような、前述の常微分方程式の表現を含み、モデル設定、数値計算、判定のステップを有するプログラムによっても提供することができる。

〔６．付録：式（１８）における粘性力の導出〕
図６に示すように液体とｚ＝０の平面２１で接する物体２２のｘ方向の速度が次式

のようにある周波数ωで振動する場合の境界層流れを考える。ｘ方向の流速をｖ_ｘ、圧力をｐとし、これらの境界層のすぐ外側の主流での値を、上添字（ｍ）を付けて表わす。

主流に関する運動方程式は

となる。

この運動方程式の解を

とおき、式（Ａ２）に代入すると、次式が得られる。

境界層内流れを解析するため、ｚ方向の流速は非常に小さいと仮定すると、ｚ方向のナビエ・ストークス方程式は

に帰着する。

従って、条件

によって次式が要求される。

境界層は非常に薄いので、近似

をｘ方向のナビエ・ストークス方程式に適用できる。このようにして次式が得られる。

式（Ａ７）を下記の境界条件の下で解く。

解を

と表わし、式（Ａ１０）と式（Ａ３）の第２式を式（Ａ７）〜（Ａ９）に代入すると、次のようになる。

式（Ａ４）を式（Ａ１１）に代入して次式を得る。

式（Ａ１４）に関する同次方程式の特性根αは、

によって定まり、解が有限である条件より、

となる。

従って、式（Ａ１４）の一般解は次のようになる。

Ｃ_ｘ（ｘ）は任意関数で、式（Ａ１３）の条件を課すことにより、次のように物体速度と主流速度で表わせる。

式（Ａ１６）を式（Ａ１５）に代入すると、

となる。粘性によってｚ＝０で接線方向に作用する力は、

である。

本文の円筒管の問題では、ｘを周方向とみて、周方向に下記の力が作用する。

虚数単位を陽に含むので

と変形し、第２項目ではｉωを時間微分作用素とみなす。また、ωには近似的に流体付加質量を考慮しない管の固有振動数を与える。

なお、上述の実施の形態は、本発明の一具体例を示すものであり、本発明を限定するものではない。例えば、管の配置は４回回転対象に限らず、例えば６回回転対象とすることがであってもよく、正多角形は正方形に限らず、例えば正六角形であってもよい。

また、この管群振動の予測が適用されるのは加圧水型原子炉の蒸気発生器に限らず、例えば一般の熱交換器やボイラに適用することもできる。また、管は中空なものに限らず、一般の構造物に適用することができ、例えば建物を風が通過する際の振動の予測にも適用することができる。

１０管群
１１管
１２セル
２０計算装置

Claims

管群振動予測方法を適用した計算装置であって、
所定径の管が所定間隔にて周期的に配置されてなる管群において、この管群の軸方向に直交して流れる流体により当該管群に発生する振動を予測するように、
前記管群の軸方向に直交する面内において、前記管群を含む領域について各管をそれぞれ含む正多角形のセルに分割するステップと、
前記正多角形のセルにおいて、この正多角形のセルの外周と、この正多角形のセルに含まれる管との境界をなす内周とについて、未知量をモード展開するステップと、
前記正多角形のセルにおいて、管の中心を始点として前記外周と前記内周を結ぶ動径上で未知量をモード展開するステップと
を有し、
前記未知量は前記流体の流速及び圧力であり、これら流速及び圧力についてナビエ・ストークスの式及び連続条件に基づいて、前記流速及び圧力の一般化座標の時間に関する常微分方程式を導き、この常微分方程式に基づいて前記流体の定常流による流れ場を決定し、
前記未知量として前記流速及び圧力にそれぞれ振動外乱を重畳し、これら流速、圧力及び振動外乱についてナビエ・ストークスの式、連続条件及び前記定常流による流れ場に基づいて、前記振動外乱の一般化座標の時間に関する常微分方程式を導き、この常微分方程式に基づいて前記管群の振動を予測し、
前記振動外乱の一般化座標の時間に関する常微分方程式に基づいて、前記管群の不安定振動が発生する前記流体の臨界流速を決定する管群予測方法を適用し、
計算装置は、
前記流速及び圧力の一般化座標の時間に関する常微分方程式と前記振動外乱の一般化座標の時間に関する常微分方程式の表現を含む所定のプログラムを格納した記憶手段と、
前記管の径、前記管を配置した間隔、前記管の単位長さ当たりの質量、前記管の固有振動数、前記流体の密度、前記流体の粘性係数の少なくとも１つの値を入力する入力手段と、
前記記憶手段に格納された前記プログラムを読み込んで実行し、前記入力手段に入力された前記少なくとも１つの値に基づいて前記臨界流速を計算し、この値が閾値を越えたかどうかを判定する演算手段と、
前記演算手段で計算した前記臨界流速と前記判定の結果を出力する出力手段と
を含むことを特徴とする計算装置。
管群振動予測方法を適用した計算装置であって、
所定径の管が所定間隔にて周期的に配置されてなる管群において、この管群の軸方向に直交して流れる流体により当該管群に発生する振動を予測するように、
前記管群の軸方向に直交する面内において、前記管群を含む領域について各管をそれぞれ含む正多角形のセルに分割するステップと、
前記正多角形のセルにおいて、この正多角形のセルの外周と、この正多角形のセルに含まれる管との境界をなす内周とについて、未知量をモード展開するステップと、
前記正多角形のセルにおいて、管の中心を始点として前記外周と前記内周を結ぶ動径上で未知量をモード展開するステップと、
前記モード展開した未知量についてガレルキン法を適用するステップと
を有し、
前記未知量は前記流体の流速及び圧力であり、これら流速及び圧力についてナビエ・ストークスの式及び連続条件に基づいて、前記流速及び圧力の一般化座標の時間に関する常微分方程式を導き、この常微分方程式に基づいて前記流体の定常流による流れ場を決定し、
前記未知量として前記流速及び圧力にそれぞれ振動外乱を重畳し、これら流速、圧力及び振動外乱についてナビエ・ストークスの式、連続条件及び前記定常流による流れ場に基づいて、前記振動外乱の一般化座標の時間に関する常微分方程式を導き、この常微分方程式に基づいて前記管群の振動を予測し、
前記振動外乱の一般化座標の時間に関する常微分方程式に基づいて、前記管群の不安定振動が発生する前記流体の臨界流速を決定する管群予測方法を適用し、
計算装置は、
前記流速及び圧力の一般化座標の時間に関する常微分方程式と前記振動外乱の一般化座標の時間に関する常微分方程式の表現を含む所定のプログラムを格納した記憶手段と、
前記管の径、前記管を配置した間隔、前記管の単位長さ当たりの質量、前記管の固有振動数、前記流体の密度、前記流体の粘性係数の少なくとも１つの値を入力する入力手段と、
前記記憶手段に格納された前記プログラムを読み込んで実行し、前記入力手段に入力された前記少なくとも１つの値に基づいて前記臨界流速を計算し、この値が閾値を越えたかどうかを判定する演算手段と、
前記演算手段で計算した前記臨界流速と前記判定の結果を出力する出力手段と
を含むことを特徴とする計算装置。