JP7698196B2 - Information processing device, information processing method, and program - Google Patents
Information processing device, information processing method, and program Download PDFInfo
- Publication number
- JP7698196B2 JP7698196B2 JP2021130708A JP2021130708A JP7698196B2 JP 7698196 B2 JP7698196 B2 JP 7698196B2 JP 2021130708 A JP2021130708 A JP 2021130708A JP 2021130708 A JP2021130708 A JP 2021130708A JP 7698196 B2 JP7698196 B2 JP 7698196B2
- Authority
- JP
- Japan
- Prior art keywords
- time
- flow velocity
- temperature
- derived
- basis vectors
- 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
Links
Images
Landscapes
- Measuring Volume Flow (AREA)
- Continuous Casting (AREA)
Description
本発明は、情報処理装置、情報処理方法、およびプログラムに関し、特に、連続鋳造設備の鋳型内の溶鋼の状態を推定するために用いて好適なものである。 The present invention relates to an information processing device, an information processing method, and a program, and is particularly suitable for use in estimating the state of molten steel in a mold of a continuous casting facility.
図18は、連続鋳造設備の概略構成の一例を示す図である。尚、各図において、l1軸、l2軸、l3軸は、向きを示すものである。○の中に×を示しているものは、紙面の手前側から奥側に向く方向を示す。l1軸、l2軸、l3軸は、三次元直交座標系の各軸である。 Fig. 18 is a diagram showing an example of a schematic configuration of a continuous casting facility. In each figure, the l1 axis, l2 axis, and l3 axis indicate directions. A circle with a cross indicates a direction from the front side to the back side of the paper. The l1 axis, l2 axis, and l3 axis are each axis of a three-dimensional orthogonal coordinate system.
取鍋1からタンディッシュ2へ供給された溶鋼3は、鋳型4へ注入される。タンディッシュ2の底部にはスライディングノズル5が設けられている。スライディングノズル5の下部には、浸漬ノズル6が設けられている。浸漬ノズル6は、鋳型4の水平断面(l1-l3断面)の中央部に配置される。浸漬ノズル6の先端の側面には、左右一対の吐出口7(l1軸方向の両端の一対の吐出口7)が形成されている。浸漬ノズル6の先端(吐出口7)は、鋳型4内(鋳型4で囲まれた領域)に供給された溶鋼3に浸漬された状態にされる。図19は、鋳型4内の溶鋼3の流れの一例を示す図である。
タンディッシュ2に供給された溶鋼3は、スライディングノズル5を介して浸漬ノズル6内を流下し、左右一対の吐出口7から鋳型4内に注入される。左右一対の吐出口7から吐出された溶鋼3は、凝固シェル8に衝突した後、図19の矢印で表されるように上昇流と下降流とに分流される。
The
図19(a)に示すように、定常時では、左右一対の吐出口7から吐出される溶鋼の量はほぼ均等になっている。しかしながら、図19(b)に示すように、左右一対の吐出口7から吐出される溶鋼の量が不均等になることがある。このように左右一対の吐出口7から吐出される溶鋼の量が不均等になることを偏流という。偏流が生じると、浸漬ノズル6を挟んだ左右両側における溶鋼3の流速に差異が生じる。
As shown in FIG. 19(a), during steady state, the amount of molten steel discharged from the pair of left and
かかる偏流が生じる原因としては、浸漬ノズル6の内面等にアルミナ等の酸化物が付着したり、吐出口7が溶損して吐出口7の形状が歪んだりすることが挙げられる。また、スライディングノズル5の構造上、溶鋼3が浸漬ノズル6内を軸対称に流下せず、左右何れかに偏って流下することも、偏流が生じる原因として挙げられる。
The causes of such drift include the adhesion of oxides such as alumina to the inner surface of the submerged
以上のような理由により鋳型4内で溶鋼3の偏流が生じると、浸漬ノズル6を挟んだ左右両側の領域のうち、溶鋼3の量が多い側の領域では、溶鋼3が凝固シェル8の内面に沿って上方および下方に勢いよく分流することになる。勢いの強い溶鋼3の上昇流は、溶鋼3の湯面(メニスカス)の盛り上がりを生起する。これにより、溶鋼3の湯面上に散布されたフラックスが鋳型4の内壁面と凝固シェル8との間に供給されるのが阻害される。そうすると、凝固シェル8の形状が不均一となりやすく、連続鋳造設備で鋳造される鋳片に発生する皺や割れ等の原因となる。
When the
また、勢いの強い溶鋼3の下降流は、溶鋼3の深くまで達して非金属介在物の浮上を妨げる。非金属性介在物がストランド深くまで進入することは、連続鋳造設備で鋳造される鋳片に非金属介在物性欠陥をもたらす等の原因となる。溶鋼3の湯面レベル(溶鋼3の湯面(メニスカス)の高さ方向(l2軸方向)の位置)は、湯面レベル計9で測定される。湯面レベル計9としては渦流センサ等が用いられる。
In addition, the powerful downward flow of the
一方、浸漬ノズル6を挟んだ左右両側の領域のうち、溶鋼3の量が少ない側の領域では、溶鋼3の吐出される流れの勢いが弱くなる。溶鋼3の吐出される流れの勢いが弱いと、吐出口7内の溶鋼3の流れによどみが発生しやすくなり、アルミナ等の析出物が浸漬ノズル6内に付着する虞がある。アルミナ等の析出物の付着は、浸漬ノズル6内の流路を閉塞する等の原因となる。
On the other hand, in the area on either side of the
以上述べたように鋳型4内に偏流が生じると、連続鋳造の操業に支障があるばかりではなく、鋳片の品質の低下を招き、好ましくない。従って、偏流が生じた場合には、例えば、鋳造の速度を変更したり、フラックスを変更したり、浸漬ノズル6を交換したりする等の処置が必要である。これらの処置を、適切かつ迅速に行うために、鋳型4内の溶鋼3の状態を、数値シミュレーションにより推定することが求められる。
As described above, if drift occurs in the
特許文献1には、溶鋼3の流速および温度の計算格子点(計算位置)におけるベクトルである状態ベクトルのデータである数値解析データを数値シミュレーションにより導出することが記載されている。しかしながら、特許文献1に記載の技術では、導出する変数の数が、計算格子点の数と物理変数との積(例えば、数千~数十万)になる。このため、溶鋼3の流速および温度を導出するのに要する時間が長くなる。よって、溶鋼3の状態をオンライン・リアルタイムで計算するためには、極めて高い処理能力を有する計算機を用いる必要がある。
そこで、非特許文献1および特許文献2に記載されているように固有直交分解(Proper Orthogonal Decomposition)による流体解析を実行することが考えられる。固有直交分解は、導出したい物理量のデータが存在するn次元の格子点ベクトルの代わりにr次元の位相空間座標系(r≪n)で当該データを表すために、当該データをr次元のモードの基底ベクトルの線形結合で表現(近似)するものである。尚、固有直交分解は、統計学の分野等では、主成分分析(Principle Component Analysis)等とも呼ばれている。
As such, it is possible to carry out fluid analysis using proper orthogonal decomposition, as described in Non-Patent
非特許文献1では、平板翼周りの流れを、固有直交分解による流体解析で導出する際に、速度ベクトルから導出される運動エネルギーの時系列データから作成したスナップショット行列を特異値分解することにより基底ベクトルを導出し、当該基底ベクトルを用いて流体の流速の各成分を表現することが開示されている。
Non-Patent
特許文献2には、空気の流れを、固有直交分解による流体解析で導出する際に、流速の各成分を1つのスナップショット行列に纏め、当該スナップショット行列を特異値分解することにより、基底ベクトルを導出することが開示されている。
しかしながら、前述したように、連続鋳造設備の鋳型4内では、浸漬ノズル6の吐出口7から吐出された溶鋼3は、凝固シェル8に衝突した後、急激に向きを変える(図19の矢印で表す上昇流と下降流とを参照)。従って、溶鋼3の流速の分布は、l1軸方向(幅方向)、l2軸方向(鋳造方向)、l3軸方向(奥行方向)で大きく異なる。このため、非特許文献1に記載されているように運動エネルギーに基づく基底ベクトルを用いて溶鋼3の流速を表現すると、より高次のモードの基底ベクトルを含めて溶鋼3の流速を表現(近似)しなければ、溶鋼3の流速を高精度に導出することが容易ではない。一方、より高次のモードの基底ベクトルを含めて溶鋼3の流速を表現(近似)すると、溶鋼3の流速を導出するのに要する時間が長くなる。
However, as described above, in the
また、特許文献2に記載されているスナップショット行列は、その各列において各計算格子点における流速の値を座標軸成分ごとにまとめて、そのまとまりを順に並べて格納したものである。従って、このスナップショット行列から固有直交分解によって得られる基底ベクトルには物理的な意味がない。この基底ベクトルの中身は、座標軸成分毎に分けることができるので、そこから各座標軸成分に対応した基底ベクトルを抽出できる。しかしながら、ここで得た各基底ベクトルは相互に正規直交化していない。このため、固有直交分解を実行するに際し、Galerkin展開によって簡素な縮約方程式(縮約シミュレータ)を構築することができず、縮約方程式が複雑になる。よって、特許文献2に記載の技術を用いて溶鋼3の流速を導出すると、溶鋼3の流速を導出するのに要する時間が長くなる。
In addition, the snapshot matrix described in
本発明は、以上のような問題点に鑑みてなされたものであり、連続鋳造設備の鋳型内の溶融金属の流速を導出するために要する時間が長くなることと、溶融金属の流速の推定精度が低下することとの双方を抑制することを目的とする。 The present invention was made in consideration of the above problems, and aims to prevent both an increase in the time required to derive the flow velocity of molten metal in a mold of a continuous casting machine and a decrease in the accuracy of estimating the flow velocity of molten metal.
本発明の情報処理装置は、連続鋳造設備の鋳型に注入される溶融金属の状態を導出することを含む処理を実行する情報処理装置であって、各計算位置および各時刻における前記溶融金属の状態を示す物理量である第1の物理量の値を格納する行列を導出する行列導出手段と、前記行列導出手段により導出された、前記第1の物理量の値を格納する行列に対して特異値分解を実行することにより、1次~M次のモードの基底ベクトルを導出することで、前記第1の物理量の1次~M次のモードの基底ベクトルを導出する基底ベクトル導出手段と、前記溶融金属の状態を示す量である状態量を導出する状態導出手段と、を有し、前記第1の物理量は、前記溶融金属の流速と、前記溶融金属の温度と、を含み、前記行列導出手段は、前記溶融金属の流速の値を格納する行列を、前記流速の軸成分ごとに導出することと、前記溶融金属の温度の値を格納する行列を導出することと、を行い、前記基底ベクトル導出手段は、前記行列導出手段により導出された、前記溶融金属の流速の値を格納する行列に対して特異値分解を実行することにより、1次~M次のモードの基底ベクトルを導出することを、前記行列導出手段により導出された、前記溶融金属の流速の値を格納する行列のそれぞれについて実行することで、前記溶融金属の流速の1次~M次のモードの基底ベクトルを、前記流速の軸成分ごとに導出することと、前記行列導出手段により導出された、前記溶融金属の温度の値を格納する行列に対して特異値分解を実行することにより、1次~M次のモードの基底ベクトルを導出することで、前記溶融金属の温度の1次~M次のモードの基底ベクトルを導出することと、を行い、前記状態導出手段は、前記基底ベクトル導出手段により前記第1の物理量の1次~M次のモードの基底ベクトルが導出された後に、前記第1の物理量の1次~M次のモードの基底ベクトルを線形結合する際に用いる係数を導出する第1の係数導出手段を有し、前記状態量は、前記第1の係数導出手段により導出された前記係数に基づいて定められ、前記第1の係数導出手段は、前記第1の物理量の1次~M次のモードの基底ベクトルを線形結合する際に用いる係数として、少なくとも前記溶融金属の流速の1次~M次のモードの基底ベクトルを線形結合する際に用いる係数を、軸成分ごとに導出することと、前記基底ベクトル導出手段により前記溶融金属の温度の1次~M次のモードの基底ベクトルが導出された後に、前記溶融金属の温度の1次~M次のモードの基底ベクトルを線形結合する際に用いる係数を導出することと、を行い、前記Mの値は、2以上の整数であり、且つ、前記流速の軸成分ごと、前記温度ごとに定まり、前記溶融金属の流速は、前記基底ベクトル導出手段により導出された、前記溶融金属の流速の1次~M次のモードの基底ベクトルを、前記第1の係数導出手段により導出された係数を用いて線形結合することにより表され、前記溶融金属の温度は、前記基底ベクトル導出手段により導出された、前記溶融金属の温度の1次~M次のモードの基底ベクトルを、前記第1の係数導出手段により導出された係数を用いて線形結合することにより表されることを特徴とする。 The information processing device of the present invention is an information processing device that executes a process including deriving a state of molten metal poured into a mold of continuous casting equipment, and includes: a matrix derivation means for deriving a matrix storing values of a first physical quantity which is a physical quantity indicating a state of the molten metal at each calculation position and each time; a basis vector derivation means for deriving basis vectors of 1st to Mth modes of the first physical quantity by performing singular value decomposition on the matrix storing the values of the first physical quantity derived by the matrix derivation means, thereby deriving basis vectors of 1st to Mth modes; and a state derivation means for deriving a state quantity which is a quantity indicating the state of the molten metal, the first physical quantity including a flow velocity of the molten metal and a temperature of the molten metal, deriving a matrix storing a value of a flow velocity for each axial component of the flow velocity , and deriving a matrix storing a value of a temperature of the molten metal, and the basis vector derivation means derives basis vectors of 1st to Mth modes by performing singular value decomposition on the matrix storing the value of the flow velocity of the molten metal derived by the matrix derivation means, for each of the matrices storing the value of the flow velocity of the molten metal derived by the matrix derivation means, thereby deriving basis vectors of 1st to Mth modes of the flow velocity of the molten metal for each axial component of the flow velocity, and deriving basis vectors of 1st to Mth modes by performing singular value decomposition on the matrix storing the value of the temperature of the molten metal derived by the matrix derivation means. and deriving basis vectors of 1st to Mth modes of the temperature of the molten metal, the state derivation means has a first coefficient derivation means for deriving a coefficient used in linearly combining the basis vectors of the 1st to Mth modes of the first physical quantity after the basis vectors of the 1st to Mth modes of the first physical quantity are derived by the basis vector derivation means, the state quantity is determined based on the coefficient derived by the first coefficient derivation means, the first coefficient derivation means derives a coefficient used in linearly combining at least the basis vectors of the 1st to Mth modes of the flow velocity of the molten metal for each axial component as a coefficient used in linearly combining the basis vectors of the 1st to Mth modes of the first physical quantity, and deriving coefficients to be used in linearly combining the basis vectors of the first to M-th modes of the molten metal temperature after the basis vectors of the first to M-th modes of the molten metal temperature have been derived, wherein the value of M is an integer of 2 or more and is determined for each axial component of the flow velocity and for each temperature , the flow velocity of the molten metal is represented by linearly combining the basis vectors of the first to M-th modes of the molten metal flow velocity derived by the basis vector derivation means, using the coefficients derived by the first coefficient derivation means, and the temperature of the molten metal is represented by linearly combining the basis vectors of the first to M-th modes of the molten metal temperature derived by the basis vector derivation means, using the coefficients derived by the first coefficient derivation means .
本発明の情報処理方法は、連続鋳造設備の鋳型に注入される溶融金属の状態を導出することを含む処理を実行する情報処理方法であって、各計算位置および各時刻における前記溶融金属の状態を示す物理量である第1の物理量の値を格納する行列を導出する行列導出工程と、前記行列導出工程により導出された、前記第1の物理量の値を格納する行列に対して特異値分解を実行することにより、1次~M次のモードの基底ベクトルを導出することで、前記第1の物理量の1次~M次のモードの基底ベクトルを導出する基底ベクトル導出工程と、前記溶融金属の状態を示す量である状態量を導出する状態導出工程と、を有し、前記第1の物理量は、前記溶融金属の流速と、前記溶融金属の温度と、を含み、前記行列導出工程は、前記溶融金属の流速の値を格納する行列を、前記流速の軸成分ごとに導出することと、前記溶融金属の温度の値を格納する行列を導出することと、を行い、前記基底ベクトル導出工程は、前記行列導出工程により導出された、前記溶融金属の流速の値を格納する行列に対して特異値分解を実行することにより、1次~M次のモードの基底ベクトルを導出することを、前記行列導出工程により導出された、前記溶融金属の流速の値を格納する行列のそれぞれについて実行することで、前記溶融金属の流速の1次~M次のモードの基底ベクトルを、前記流速の軸成分ごとに導出することと、前記行列導出工程により導出された、前記溶融金属の温度の値を格納する行列に対して特異値分解を実行することにより、1次~M次のモードの基底ベクトルを導出することで、前記溶融金属の温度の1次~M次のモードの基底ベクトルを導出することと、を行い、前記状態導出工程は、前記基底ベクトル導出工程により前記第1の物理量の1次~M次のモードの基底ベクトルが導出された後に、前記第1の物理量の1次~M次のモードの基底ベクトルを線形結合する際に用いる係数を導出する第1の係数導出工程を有し、前記状態量は、前記第1の係数導出工程により導出された前記係数に基づいて定められ、前記第1の係数導出工程は、前記第1の物理量の1次~M次のモードの基底ベクトルを線形結合する際に用いる係数として、少なくとも前記溶融金属の流速の1次~M次のモードの基底ベクトルを線形結合する際に用いる係数を、軸成分ごとに導出することと、前記基底ベクトル導出工程により前記溶融金属の温度の1次~M次のモードの基底ベクトルが導出された後に、前記溶融金属の温度の1次~M次のモードの基底ベクトルを線形結合する際に用いる係数を導出することと、を行い、前記Mの値は、2以上の整数であり、且つ、前記流速の軸成分ごと、前記温度ごとに定まり、前記溶融金属の流速は、前記基底ベクトル導出工程により導出された、前記溶融金属の流速の1次~M次のモードの基底ベクトルを、前記第1の係数導出工程により導出された係数を用いて線形結合することにより表され、前記溶融金属の温度は、前記基底ベクトル導出工程により導出された、前記溶融金属の温度の1次~M次のモードの基底ベクトルを、前記第1の係数導出工程により導出された係数を用いて線形結合することにより表されることを特徴とする。 The information processing method of the present invention is an information processing method for executing a process including deriving a state of molten metal poured into a mold of continuous casting equipment, the information processing method including: a matrix derivation step of deriving a matrix storing values of a first physical quantity which is a physical quantity indicating a state of the molten metal at each calculation position and each time; a basis vector derivation step of deriving basis vectors of first to Mth modes of the first physical quantity by performing singular value decomposition on the matrix storing the values of the first physical quantity derived by the matrix derivation step to derive basis vectors of first to Mth modes; and a state derivation step of deriving a state quantity which is a quantity indicating the state of the molten metal, the first physical quantity including a flow velocity of the molten metal and a temperature of the molten metal, The method includes the steps of: deriving a matrix storing a value of a flow velocity for each axial component of the flow velocity; and deriving a matrix storing a value of a temperature of the molten metal; the basis vector derivation step derives basis vectors of 1st to Mth modes by performing singular value decomposition on the matrix storing the value of the flow velocity of the molten metal derived in the matrix derivation step, for each of the matrices storing the value of the flow velocity of the molten metal derived in the matrix derivation step, thereby deriving basis vectors of 1st to Mth modes of the flow velocity of the molten metal for each axial component of the flow velocity; and deriving basis vectors of 1st to Mth modes by performing singular value decomposition on the matrix storing the value of the temperature of the molten metal derived in the matrix derivation step. and deriving basis vectors of 1st to Mth modes of the temperature of the molten metal. The state derivation step includes a first coefficient derivation step of deriving a coefficient used in linearly combining the basis vectors of the 1st to Mth modes of the first physical quantity after the basis vectors of the 1st to Mth modes of the first physical quantity are derived by the basis vector derivation step, and the state quantity is determined based on the coefficient derived by the first coefficient derivation step. The first coefficient derivation step derives a coefficient used in linearly combining at least the basis vectors of the 1st to Mth modes of the flow velocity of the molten metal for each axial component as a coefficient used in linearly combining the basis vectors of the 1st to Mth modes of the first physical quantity. and deriving coefficients to be used in linearly combining the basis vectors of the first to M-th modes of the molten metal temperature after the basis vectors of the first to M-th modes of the molten metal temperature have been derived, wherein the value of M is an integer of 2 or more and is determined for each axial component of the flow velocity and for each temperature , the flow velocity of the molten metal is represented by linearly combining the basis vectors of the first to M-th modes of the molten metal flow velocity derived in the basis vector derivation step, using the coefficients derived in the first coefficient derivation step, and the temperature of the molten metal is represented by linearly combining the basis vectors of the first to M-th modes of the molten metal temperature derived in the basis vector derivation step, using the coefficients derived in the first coefficient derivation step .
本発明のプログラムは、前記情報処理装置の各手段としてコンピュータを機能させるためのものである。 The program of the present invention is intended to cause a computer to function as each of the means of the information processing device.
本発明によれば、連続鋳造設備の鋳型内の溶融金属の流速を導出するために要する時間が長くなることと、溶融金属の流速の推定精度が低下することとの双方を抑制することができる。 The present invention can prevent both an increase in the time required to derive the flow velocity of molten metal in a mold of a continuous casting machine and a decrease in the accuracy of estimating the flow velocity of molten metal.
以下、図面を参照しながら、本発明の実施形態を説明する。
〔第1の実施形態〕
まず、第1の実施形態を説明する。
図1は、情報処理装置100の機能的な構成の一例を示す図である。情報処理装置100のハードウェアは、例えば、CPU、ROM、RAM、HDD、および各種のインターフェースを有する情報処理装置、または専用のハードウェアを用いることにより実現される。情報処理装置100は、連続鋳造設備の鋳型4内における溶鋼3(溶融金属)の状態として、溶鋼3の流動および温度を導出することを含む処理を実行する。本実施形態では、情報処理装置100は、連続鋳造設備の鋳型4内における溶鋼3の流動および温度を可視化するためのデータを作成する場合を例示する。本実施形態において、連続鋳造鋳型内の全部または一部の領域における溶鋼(場合によっては凝固シェル)の物理量(流速や温度など)の分布を導出することを可視化と称し、可視化によって溶鋼の物理量の分布を導出する領域を可視化対象領域と称する。また、各実施形態において、可視化(推定)の対象となる溶鋼の状態を示す物理量を第1の物理量とも称する。また、各実施形態において、センサより測定される連続鋳造設備の各観測位置における物理量を第2の物理量とも称する。
Hereinafter, an embodiment of the present invention will be described with reference to the drawings.
First Embodiment
First, the first embodiment will be described.
FIG. 1 is a diagram showing an example of a functional configuration of an
本実施形態では、図18に示した連続鋳造設備の鋳型4内における溶鋼3の流動および温度の状態を可視化する場合を例示する。鋳型4内には、温度計(熱電対)が埋め込まれている。図2は、鋳型4内に埋め込まれる温度計L1~L12、F1~F12の一例を示す図である。
図2(a)~図2(c)に示すように、鋳型4の或る高さ位置(鋳造方向(l2軸方向)の位置)において、鋳型4の一方の長辺側の面(「F面」と称する)には、鋳型4の長辺側の幅方向(l1軸方向)の中心を対称点とする点対称となる位置に、複数の温度計F1、F3、F5、F7、F9、F11が埋め込まれている。また、鋳型4の当該高さ位置(鋳造方向(l2軸方向)の位置)において、鋳型4の他方の長辺側の面(「L面」と称する)には、鋳型4の長辺側の幅方向(l1軸方向)の中心を対称点とする点対称となる位置に、複数の温度計L1、L3、L5、L7、L9、L11が埋め込まれている。
In this embodiment, a case where the flow and temperature state of
As shown in Figures 2(a) to 2(c), at a certain height position (position in the casting direction ( l2 axis direction)) of the
また、鋳型4のF面において、複数の温度計F1、F3、F5、F7、F9、F11の下方には、複数の温度計F2、F4、F6、F8、F10、F12が埋め込まれている。複数の温度計F1~F12の鋳型4の奥行方向(l3軸方向)の位置は同じである。また、複数の温度計F1、F3、F5、F7、F9、F11のl1軸方向の間隔と、複数の温度計F2、F4、F6、F8、F10、F12のl1軸方向の間隔は、同じである。また、温度計F1、F2のl1軸方向の位置、温度計F3、F4のl1軸方向の位置、温度計F5、F6のl1軸方向の位置、温度計F7、F8のl1軸方向の位置、温度計F9、F10のl1軸方向の位置、および温度計F11、F12のl1軸方向の位置は、同じである。
In addition, on the F surface of the
同様に、鋳型4のL面において、複数の温度計L1、L3、L5、L7、L9、L11の下方には、複数の温度計L2、L4、L6、L8、L10、L12が埋め込まれている。複数の温度計L1~L12の鋳型4の奥行方向(l3軸方向)の位置は同じである。また、複数の温度計L1、L3、L5、L7、L9、L11のl1軸方向の間隔と、複数の温度計L2、L4、L6、L8、L10、L12のl1軸方向の間隔は、同じである。また、温度計L1、L2のl1軸方向の位置、温度計L3、L4のl1軸方向の位置、温度計L5、L6のl1軸方向の位置、温度計L7、L8のl1軸方向の位置、温度計L9、L10のl1軸方向の位置、および温度計L11、L12のl1軸方向の位置は、同じである。
Similarly, on the L surface of the
また、鋳型4のF面において上側にある複数の温度計F1、F3、F5、F7、F9、F11のl2軸方向の位置と、鋳型4のL面において上側にある複数の温度計L1、L3、L5、L7、L9、L11のl2軸方向の位置は、全て同じである。鋳型4のF面において下側にある複数の温度計F2、F4、F6、F8、F10、F12のl2軸方向の位置と、鋳型4のL面において下側にある複数の温度計L2、L4、L6、L8、L10、L12のl2軸方向の位置は、全て同じである。
In addition, the positions in the l2 axis direction of the multiple thermometers F1, F3, F5, F7, F9, and F11 on the upper side of the F surface of the
また、温度計F1、L1のl1軸方向の位置、温度計F3、L3のl1軸方向の位置、温度計F5、L5のl1軸方向の位置、温度計F7、L7のl1軸方向の位置、温度計F9、L9のl1軸方向の位置、温度計F11、L11のl1軸方向の位置、温度計F2、L2のl1軸方向の位置、温度計F4、L4のl1軸方向の位置、温度計F6、L6のl1軸方向の位置、温度計F8、L8のl1軸方向の位置、温度計F10、L10のl1軸方向の位置、および温度計F12、L12のl1軸方向の位置は、同じである。
また、温度計F1、F2、F3、F4、F5、F6、F7、F8、F9、F10、F11、F12のF面(内壁面)からのl3軸方向の距離と、温度計L1、L2、L3、L4、L5、L6、L7、L8、L9、L10、L11、L12のL面(内壁面)からのl3軸方向の距離とは、同じである。
以下の説明では、複数の温度計F1~F12、L1~L12を、必要に応じて温度計F、Lと総称する。
In addition, the positions of thermometers F1, L1 in the l1 axis direction, the positions of thermometers F3, L3 in the l1 axis direction, the positions of thermometers F5, L5 in the l1 axis direction, the positions of thermometers F7, L7 in the l1 axis direction, the positions of thermometers F9, L9 in the l1 axis direction, the positions of thermometers F11, L11 in the l1 axis direction, the positions of thermometers F2, L2 in the l1 axis direction, the positions of thermometers F4, L4 in the l1 axis direction, the positions of thermometers F6, L6 in the l1 axis direction, the positions of thermometers F8, L8 in the l1 axis direction, the positions of thermometers F10, L10 in the l1 axis direction, and the positions of thermometers F12, L12 in the l1 axis direction are the same.
In addition, the distance in the l3 axis direction from the F surface (inner wall surface) of thermometers F1, F2, F3, F4, F5, F6, F7, F8, F9, F10, F11, and F12 is the same as the distance in the l3 axis direction from the L surface (inner wall surface) of thermometers L1, L2, L3, L4, L5, L6, L7, L8, L9, L10, L11, and L12.
In the following description, the multiple thermometers F1 to F12 and L1 to L12 will be collectively referred to as thermometers F and L as necessary.
特許文献1に記載されている技術において、連続鋳造設備の鋳型4内における溶鋼3の流動および温度の状態をオンライン・リアルタイムで可視化する場合には、溶鋼3の流速および温度を数値シミュレーションすることと、数値シミュレーションの結果を利用したデータ同化等の処理とを、(数値シミュレーションで)想定している時間の経過と同程度の処理速度で実行する必要がある。このため、溶鋼3の流速および温度を数値シミュレーションすることを、計算精度を大きく落とすことなく高速に実行することが求められる。
In the technology described in
そこで、本発明者らは、連続鋳造設備の鋳型4内における溶鋼3の流速の分布は、航空機や車の周囲の流速の分布のように等方的に乱れておらず、l1軸方向(幅方向)、l2軸方向(鋳造方向)、l3軸方向(奥行方向)で大きく異なることに着目した。そして、本発明者らは、溶鋼3の流速をその軸成分ごとに固有直交分解することにより、溶鋼3の流速を数値シミュレーションする際の、計算負荷の低減と計算精度の低下の抑制との双方を達成することができることを見出した。以下に、本実施形態の情報処理装置100が有する機能の一例を説明する。本実施形態では、特許文献1と同様に、鋳型4の奥行方向(l3軸方向)の中心の位置において、鋳型4の幅方向(l1軸方向)と鋳造方向(l2軸方向)とにより定まる2次元平面(l1-l2平面)における溶鋼3の流動および温度を可視化する場合を例示する。尚、図2(a)に示すように、鋳型4の幅方向(l1軸方向)および奥行方向(l3軸方向)の中心の位置は、浸漬ノズル6の軸の位置に対応する。以下の説明では、この2次元平面を必要に応じて可視化対象面と称する。本実施形態では、可視化対象面が可視化対象領域の一例になる。
Therefore, the inventors have noticed that the distribution of the flow velocity of the
図1において、数値シミュレーション部101、行列導出部102、基底ベクトル導出部103、および縮約シミュレータ導出部104は、オフラインで実行される。一方、状態導出部111(第1の係数導出部111aおよび物理量導出部111b)、第1の確率密度関数導出部112、観測データ取得部113、尤度関数導出部114、データ同化部115、可視化データ作成部116、および出力部117は、オンラインで実行される。オンラインで実行される処理は、連続鋳造設備の実際の操業時に、連続鋳造設備において測定された時系列データ(測定値)を、通信回線を通じて入力しながら実行される処理である。このオンラインで実行される処理は、リアルタイムで行われることが好ましい。リアルタイムで実行される処理は、観測データが取得される時刻から当該時刻での溶鋼の状態(即ち、流速と温度との少なくとも一方)を出力する時刻までのタイムラグが、偏流などの異常な溶鋼の状態を示す出力に対して適切かつ迅速な処置を行うことが可能な程度である処理である。オフラインで実行される処理は、オンライン(好ましくは、オンライン・リアルタイム)で実行される処理の準備のための処理であり、オンラインで実行される処理よりも前に事前に、連続鋳造設備において測定された時系列データ(測定値)を入力しない状態で実行される処理である。
1, the
(数値シミュレーション部101)
数値シミュレーション部101は、可視化の対象となる連続鋳造設備における非定常の流動および伝熱の数値シミュレーションを、連続鋳造設備の操業パラメータおよび溶鋼3の物性値を用いてオフラインで実行する。連続鋳造設備の操業パラメータおよび溶鋼3の物性値は、予め数値シミュレーション部101が記憶しているものである。
(Numerical Simulation Unit 101)
The
溶鋼3内の流動および伝熱の数値シミュレーションは、ナビエストークス方程式(運動量保存式)を連続の式(質量保存式)およびエネルギー保存式と共に解いて求めることができる。例えば、有限体積法を用いて各式を離散化して解くことにより数値シミュレーションは実行される。
数値シミュレーションを実行する際には、初期条件と境界条件とを設定する必要がある。流動および伝熱の境界条件は、鋳型4の内壁面と、溶鋼3の湯面と、浸漬ノズル6の外壁面と、吐出口7と、凝固シェル8と溶鋼3との境界面と、溶鋼3の下端とにおいて設定される。溶鋼3の下端とは、可視化対象領域のうち最も下の位置(l2軸の正の方向の端の位置)を指す。流動および伝熱の初期条件は、溶鋼3において設定される。流動および伝熱の境界条件は、連続鋳造設備において想定されている操業条件に従って設定すればよい。本実施形態では、数値シミュレーション部101は、測定値を用いずに数値シミュレーションを実行する場合を例示する。
A numerical simulation of the flow and heat transfer in the
When performing a numerical simulation, it is necessary to set initial conditions and boundary conditions. The boundary conditions of flow and heat transfer are set at the inner wall surface of the
数値シミュレーション部101は、時刻tでの計算格子点における溶鋼3の流速および温度を導出する。本実施形態では、計算格子点が計算位置の一例となる。尚、以下の説明では、溶鋼3の流速、溶鋼3の温度を必要に応じて流速、温度と略称する。また、流速は、l1軸方向、l2軸方向、l3軸方向の3成分の値を有する。前述したように本実施形態では、2次元平面(l1-l2平面)における溶鋼3の流動および温度を可視化する場合を例示するため、以下の説明では、溶鋼3の流動および温度を可視化するためのモデルとして2次元のモデルを仮定して、流速のl1軸成分を必要に応じてuと表記し、流速のl2軸成分を必要に応じてvと表記し、温度を必要に応じてTと表記する。また、計算格子点(の位置)を必要に応じてlと表記する。計算格子点lのl1-l2平面の座標は(l1,l2)となる。
The
本実施形態では、数値シミュレーション部101は、ステップ幅Δt1の時間隔の各時刻tで各計算格子点lにおける流速u(l,t)、v(l,t)および温度T(l,t)を導出する。尚、連続鋳造設備に対するナビエストークス方程式、連続の式、およびエネルギー保存式や、初期条件および境界条件の内容については、例えば、非特許文献3に記載されているので、ここでは、その詳細な説明を省略する。
In this embodiment, the
(行列導出部102)
行列導出部102は、数値シミュレーション部101によりステップ幅Δt1の時間隔の各時刻tで導出された、各計算格子点lにおける流速u(l,t)、v(l,t)、温度T(l,t)に基づいて、以下の(1)式、(2)式、(3)式に示すように、流速uのスナップショット行列MU、流速vのスナップショット行列MV、温度Tのスナップショット行列MTをそれぞれ導出する。
(Matrix derivation unit 102)
The
(1)式~(3)式において、行番号1、・・・、lMは、計算格子点lを示す。(1)式~(3)式では、計算格子点lの数がlM個である場合を例示する。列番号1、・・・、tMは、時刻tを表す。このように、流速uのスナップショット行列MU、流速vのスナップショット行列MV、および温度Tのスナップショット行列MTにおいて、同一の行に含まれる成分は、計算格子点lが同一の成分であり、同一の列に含まれる成分は、時刻が同一の成分である。
In formulas (1) to (3),
行列導出部102は、数値シミュレーション部101によりステップ幅Δt1の時間隔の各時刻tで導出された、各計算格子点lにおける流速u(l,t)から、ステップ幅Δt1よりも長いステップ幅Δt2の時間隔の各時刻tでの流速u(l,t)を抽出する。そして、行列導出部102は、抽出した流速u(l,t)が成分として格納されたスナップショット行列を流速uのスナップショット行列MUとして導出する。尚、行列導出部102は、ステップ幅Δt1の時間隔の各時刻tで導出された、各計算格子点lにおける流速u(l,t)を補間してから、ステップ幅Δt2の時間隔の各時刻tでの流速u(l,t)を抽出してもよい。このように処理することによって、ステップ幅Δt2をステップ幅Δt1の整数倍とならないように設定した場合でも、ステップ幅Δt2の時間隔の各時刻tでの流速u(l,t)を抽出することができる。
The
また、行列導出部102は、数値シミュレーション部101によりステップ幅Δt1の時間隔の各時刻tで導出された、各計算格子点lにおける流速v(l,t)から、ステップ幅Δt1よりも長いステップ幅Δt2の時間隔の各時刻tでの流速v(l,t)を抽出する。そして、行列導出部102は、抽出した流速v(l,t)が成分として格納されたスナップショット行列を流速vのスナップショット行列MVとして導出する。尚、行列導出部102は、ステップ幅Δt1の時間隔の各時刻tで導出された、各計算格子点lにおける流速v(l,t)を補間してから、ステップ幅Δt2の時間隔の各時刻tでの流速v(l,t)を抽出してもよい。このように処理することによって、ステップ幅Δt2をステップ幅Δt1の整数倍とならないように設定した場合でも、ステップ幅Δt2の時間隔の各時刻tでの流速v(l,t)を抽出することができる。
The
また、行列導出部102は、数値シミュレーション部101によりステップ幅Δt1の時間隔の各時刻tで導出された、各計算格子点lにおける温度T(l,t)から、ステップ幅Δt1よりも長いステップ幅Δt2の時間隔の各時刻tでの温度T(l,t)を抽出する。そして、行列導出部102は、抽出した温度T(l,t)が成分として格納されたスナップショット行列を温度Tのスナップショット行列MTとして導出する。尚、行列導出部102は、ステップ幅Δt1の時間隔の各時刻tで導出された、各計算格子点lにおける温度T(l,t)を補間してから、ステップ幅Δt2の時間隔の各時刻tでの温度T(l,t)を抽出してもよい。このように処理することによって、ステップ幅Δt2をステップ幅Δt1の整数倍とならないように設定した場合でも、ステップ幅Δt2の時間隔の各時刻tでの流速T(l,t)を抽出することができる。
The
このように、本実施形態では、数値シミュレーション部101によりステップ幅Δt1の時間隔の各時刻tで導出された、各計算格子点lにおける流速u(l,t)、v(l,t)、温度T(l,t)から、ステップ幅Δt1よりも長いステップ幅Δt2の時間隔の各時刻tでの流速u(l,t)、v(l,t)、温度T(l,t)を抽出して、流速u、v、温度Tのスナップショット行列MU、MV、MTを導出する。従って、ステップ幅Δt1の時間隔の各時刻tでの流速u(l,t)、v(l,t)、温度T(l,t)が成分として格納されたスナップショット行列を扱う場合よりも、後述する基底ベクトル導出部103における計算負荷を低減することができる。
このような観点から、スナップショット行列には、想定する流動や熱輸送の変化の特徴を捉えるのに十分なデータが含まれている必要があるが、一般的にステップ幅Δt1は、数値シミュレーションを安定に維持するための時間隔であり、流動や熱輸送の変化の特徴を捉える目的としては過小である。また、ステップ幅Δt2が大きすぎると、流動や熱輸送の変化の特徴を十分に捉えることができない。このような観点から、ステップ幅Δt2はステップ幅Δt1の10倍以上であり、対象とする連続鋳造設備における計測データ(後述する観測データ取得部113で取得される観測データ)の取得間隔よりも小さいことが好ましい。
In this manner, in this embodiment, the flow velocities u(l,t), v(l,t) and temperature T(l,t) at each time t of a time interval of a step width Δt2 longer than the step width Δt1 are extracted from the flow velocities u(l,t), v(l,t) and temperature T(l,t) at each computational grid point l derived at each time t of a time interval of a step width Δt1 by the
From this viewpoint, the snapshot matrix needs to contain sufficient data to capture the characteristics of the expected changes in flow and heat transport, but the step width Δt1 is generally a time interval for maintaining the stability of the numerical simulation, and is too small for the purpose of capturing the characteristics of changes in flow and heat transport. Moreover, if the step width Δt2 is too large, the characteristics of changes in flow and heat transport cannot be captured sufficiently. From this viewpoint, it is preferable that the step width Δt2 is 10 times or more the step width Δt1 and is smaller than the acquisition interval of the measurement data (observation data acquired by the observation data acquisition unit 113 described later) in the target continuous casting equipment.
(基底ベクトル導出部103)
前述したように本実施形態では、流速u(l,t)、v(l,t)、温度T(l,t)のそれぞれについて個別に固有直交分解を実行する。この場合、流速u(l,t)、v(l,t)、温度T(l,t)は、それぞれ、以下の(4)式、(5)式、(6)式で表される。
(Basis Vector Derivation Unit 103)
As described above, in this embodiment, proper orthogonal decomposition is performed for each of the flow velocities u(l,t), v(l,t), and temperature T(l,t). In this case, the flow velocities u(l,t), v(l,t), and temperature T(l,t) are expressed by the following formulas (4), (5), and (6), respectively.
(4)式において、φu j(l)は、流速uの基底ベクトルである。流速uの基底ベクトルφu j(l)は、流速uのスナップショット行列MUを特異値分解することにより得られる左特異ベクトルである。φu j(l)のuは、流速uの基底ベクトルであることを示し、jは、固有直交分解のモード番号を示す。特異値が大きいものから順に、モード番号jとして1から昇順に正の整数が与えられる。例えば、最大の特異値と共に導出される左特異ベクトルである流速uの基底ベクトルは、φu 1(l)となる。流速uの基底ベクトルφu j(l)は、相互に正規直交系をなす。 In equation (4), φ u j (l) is a basis vector of the flow velocity u. The basis vector φ u j (l) of the flow velocity u is a left singular vector obtained by singular value decomposition of the snapshot matrix MU of the flow velocity u. The u in φ u j (l) indicates that it is a basis vector of the flow velocity u, and j indicates the mode number of the proper orthogonal decomposition. Positive integers are given as the mode number j in ascending order from 1 in descending order of the singular value. For example, the basis vector of the flow velocity u, which is a left singular vector derived together with the maximum singular value, is φ u 1 (l). The basis vectors φ u j (l) of the flow velocity u form a mutually orthogonal system.
(4)式に示すように、固有直交分解においては、流速u(l,t)は、流速uの基底ベクトルφu j(l)の線形結合で表される。au j(t)は、流速uの基底ベクトルφu j(l)を線形結合する際に用いる係数である。具体的にau j(t)は、流速uの基底ベクトルφu j(l)に乗算される係数である。(4)式における流速uがスナップショット行列MUで表現される場合には、係数au j(t)は、流速uのスナップショット行列MUを特異値分解することにより得られるモード番号jの特異値とモード番号jの特異値に対する右特異ベクトルとの積として導出される。尚、係数au j(t)のu、jの意味は、それぞれ、φu j(l)のu、jの意味と同じである。 As shown in formula (4), in the proper orthogonal decomposition, the flow velocity u(l, t) is expressed as a linear combination of the basis vectors φ u j (l) of the flow velocity u. a u j (t) is a coefficient used when linearly combining the basis vectors φ u j (l) of the flow velocity u. Specifically, a u j (t) is a coefficient multiplied by the basis vector φ u j (l) of the flow velocity u. When the flow velocity u in formula (4) is expressed by a snapshot matrix MU, the coefficient a u j (t) is derived as the product of the singular value of the mode number j obtained by singular value decomposition of the snapshot matrix MU of the flow velocity u and the right singular vector for the singular value of the mode number j. Note that the meanings of u and j in the coefficient a u j (t) are the same as those of u and j in φ u j (l), respectively.
また、(5)式において、φv j(l)は、流速vの基底ベクトルである。流速vの基底ベクトルφv j(l)は、流速vのスナップショット行列MVを特異値分解することにより得られる左特異ベクトルである。φv j(l)のvは、流速vの基底ベクトルであることを示し、jは、固有直交分解のモード番号を示す。特異値が大きいものから順に、モード番号jとして1から昇順に正の整数が与えられる。流速vの基底ベクトルφv j(l)は、相互に正規直交系をなす。(5)式に示すように、固有直交分解においては、流速v(l,t)は、流速vの基底ベクトルφv j(l)の線形結合で表される。av j(t)は、流速vの基底ベクトルφv j(l)を線形結合する際に用いる係数である。具体的にav j(t)は、流速uの基底ベクトルφv j(l)に乗算される係数である。(5)式における流速vがスナップショット行列MVで表現される場合には、係数av j(t)は、流速vのスナップショット行列MVを特異値分解することにより得られるモード番号jの特異値とモード番号jの特異値に対する右特異ベクトルとの積として導出される。尚、係数av j(t)のv、jの意味は、それぞれ、φv j(l)のv、jの意味と同じである。 In addition, in formula (5), φ v j (l) is a basis vector of the flow velocity v. The basis vector φ v j (l) of the flow velocity v is a left singular vector obtained by singular value decomposition of the snapshot matrix MV of the flow velocity v. The v in φ v j (l) indicates that it is a basis vector of the flow velocity v, and j indicates the mode number of the proper orthogonal decomposition. Positive integers are given as the mode number j in ascending order from 1 in descending order of the singular value. The basis vectors φ v j (l) of the flow velocity v form a mutually orthogonal system. As shown in formula (5), in proper orthogonal decomposition, the flow velocity v (l, t) is expressed as a linear combination of the basis vectors φ v j (l) of the flow velocity v. a v j (t) is a coefficient used when linearly combining the basis vectors φ v j (l) of the flow velocity v. Specifically, a v j (t) is a coefficient by which the basis vector φ v j (l) of the flow velocity u is multiplied. When the flow velocity v in equation (5) is expressed by the snapshot matrix MV, the coefficient a v j (t) is derived as the product of the singular value of mode number j obtained by singular value decomposition of the snapshot matrix MV of the flow velocity v and the right singular vector for the singular value of mode number j. Note that the meanings of v and j in the coefficient a v j (t) are the same as those of v and j in φ v j (l), respectively.
また、(6)式において、φT j(l)は、温度Tの基底ベクトルである。温度Tの基底ベクトルφT j(l)は、温度Tのスナップショット行列MTを特異値分解することにより得られる左特異ベクトルである。φT j(l)のTは、温度Tの基底ベクトルであることを示し、jは、固有直交分解のモード番号を示す。特異値が大きいものから順に、モード番号jとして1から昇順に正の整数が与えられる。温度Tの基底ベクトルφT j(l)は、相互に正規直交系をなす。(6)式に示すように、固有直交分解においては、温度T(l,t)は、温度Tの基底ベクトルφT j(l)の線形結合で表される。aT j(t)は、温度Tの基底ベクトルφT j(l)を線形結合する際に用いる係数である。具体的にaT j(t)は、温度Tの基底ベクトルφT j(l)に乗算される係数である。(6)式における温度Tがスナップショット行列MTで表現される場合には、係数aT j(t)は、温度Tのスナップショット行列MTを特異値分解することにより得られるモード番号jの特異値とモード番号jの特異値に対する右特異ベクトルとの積として導出される。尚、係数aT j(t)のT、jの意味は、それぞれ、φT j(l)のT、jの意味と同じである。 In addition, in formula (6), φ T j (l) is a basis vector of temperature T. The basis vector φ T j (l) of temperature T is a left singular vector obtained by singular value decomposition of the snapshot matrix MT of temperature T. T in φ T j (l) indicates that it is a basis vector of temperature T, and j indicates the mode number of proper orthogonal decomposition. Positive integers are given as the mode number j in ascending order from 1 in descending order of the singular value. The basis vectors φ T j (l) of temperature T form a mutually orthogonal system. As shown in formula (6), in proper orthogonal decomposition, temperature T (l, t) is expressed as a linear combination of the basis vectors φ T j (l) of temperature T. a T j (t) is a coefficient used when linearly combining the basis vectors φ T j (l) of temperature T. Specifically, a T j (t) is a coefficient multiplied by the basis vector φ T j (l) of temperature T. When the temperature T in equation (6) is expressed by the snapshot matrix MT, the coefficient a T j (t) is derived as the product of the singular value of mode number j obtained by singular value decomposition of the snapshot matrix MT of the temperature T and the right singular vector for the singular value of mode number j. Note that the meanings of T and j in the coefficient a T j (t) are the same as those of T and j in φ T j (l), respectively.
ここで、流速u、v、温度Tの基底ベクトルφu j(l)、φv j(l)、φT j(l)は、計算格子点lの位置(空間座標)のみの関数である。一方、係数au j(t)、av j(t)、aT j(t)は、時刻tのみの関数である。即ち、(4)式~(6)式に示すように、固有直交分解では、第1の物理量(流速u、v、温度T)を時刻と空間について変数分離する。係数au j(t)、av j(t)、aT j(t)は、基底ベクトルφu j(l)、φv j(l)、φT(l)の寄与度に対応し、特異値は、係数au j(t)、av j(t)、aT j(t)の時間変動の振幅に対応する。 Here, the basis vectors φ u j (l), φ v j (l), and φ T j (l) of the flow velocity u, v, and temperature T are functions of only the position (spatial coordinates) of the computational grid point l. On the other hand, the coefficients a u j (t), a v j (t), and a T j (t) are functions of only the time t. That is, as shown in equations (4) to (6), in proper orthogonal decomposition, the first physical quantity (flow velocity u, v, temperature T) is separated into variables for time and space. The coefficients a u j (t), a v j (t), and a T j (t) correspond to the contributions of the basis vectors φ u j (l), φ v j (l), and φ T (l), and the singular values correspond to the amplitude of the time fluctuations of the coefficients a u j (t), a v j (t), and a T j (t).
以上のように本実施形態では、基底ベクトル導出部103は、行列導出部102により導出された、流速uのスナップショット行列MUを特異値分解することにより得られる左特異ベクトルを、流速uの基底ベクトルφu
j(l)として導出する。
また、基底ベクトル導出部103は、行列導出部102により導出された、流速vのスナップショット行列MVを特異値分解することにより得られる左特異ベクトルを、流速vの基底ベクトルφv
j(l)として導出する。
また、基底ベクトル導出部103は、行列導出部102により導出された、温度Tのスナップショット行列MTを特異値分解することにより得られる左特異ベクトルを、温度Tの基底ベクトルφT
j(l)として導出する。
As described above, in this embodiment, the basis
Furthermore, the basis
Furthermore, the basis
(4)式において、Muは、(4)式の計算において用いられる、流速uの基底ベクトルφu j(l)の数である。即ち、(4)式では、1次~Mu次のモードの流速uの基底ベクトルφu 1(l)~φu Mu(l)を線形結合することにより流速u(l,t)が導出されることを表す。(5)式において、Mvは、(5)式の計算において用いられる、流速vの基底ベクトルφv j(l)の数である。即ち、(5)式では、1次~Mv次のモードの流速vの基底ベクトルφv 1(l)~φv Mv(l)を線形結合することにより流速v(l,t)が導出されることを表す。(6)式において、MTは、(6)式の計算において用いられる、温度Tの基底ベクトルφT j(l)の数である。即ち、(6)式では、1次~MT次のモードの温度Tの基底ベクトルφT 1(l)~φT MT(l)を線形結合することにより温度T(l,t)が導出されることを表す。 In formula (4), M u is the number of basis vectors φ u j (l) of flow velocity u used in the calculation of formula (4). That is, formula (4) indicates that flow velocity u(l, t ) is derived by linearly combining basis vectors φ u 1 (l) to φ u Mu (l) of flow velocity u in the first to M u modes. In formula (5), M v is the number of basis vectors φ v j (l) of flow velocity v used in the calculation of formula (5). That is, formula (5) indicates that flow velocity v(l, t) is derived by linearly combining basis vectors φ v 1 (l) to φ v Mv (l) of flow velocity v in the first to M v modes. In formula (6), M T is the number of basis vectors φ T j (l) of temperature T used in the calculation of formula (6). That is, equation (6) indicates that the temperature T(l, t) is derived by linearly combining the basis vectors φ T 1 (l) to φ T MT (l) of the temperature T of the 1st to M Tth modes.
このように本実施形態では、基底ベクトル導出部103は、流速u、v、温度Tの基底ベクトルφu
j(l)、φv
j(l)、φT
j(l)のうち、それぞれMu個、Mv個、MT個、の基底ベクトルφu
1(l)~φu
Mu(l)、φv
1(l)~φv
Mv(l)、φT
1(l)~φT
MT(l)を、(4)式、(5)式、(6)式の計算に用いる。(4)式、(5)式、(6)式の計算に用いる流速u、v、温度Tの基底ベクトルφu
j(l)、φv
j(l)、φT
j(l)の数Mu、Mv、MTが取り得る値の最大値(即ち、最大次数)は、それぞれ、流速uのスナップショット行列MUの列数、流速vのスナップショット行列MVの列数、温度Tのスナップショット行列MTの列数となる。
Thus, in this embodiment, the basis
ここで、(4)式、(5)式、(6)式の計算に用いる流速u、v、温度Tの基底ベクトルφu
j(l)、φv
j(l)、φT
j(l)の数Mu、Mv、MTを導出する方法の一例を説明する。
流速u、v、温度Tの基底ベクトルφu
j(l)、φv
j(l)、φT
j(l)を、特異値が大きいものから順に並べると、それらの大きさは指数関数的に減少する。従って、低次のモードの基底ベクトルφu
j(l)、φv
j(l)、φT
j(l)と比較すると、高次のモードの基底ベクトルφu
j(l)、φv
j(l)、φT
j(l)の、流速u、v、温度Tに対する寄与は低くなる。このため、流速u、v、温度Tは、低次のモードの基底ベクトルφu
j(l)、φv
j(l)、φT
j(l)だけの線形結合で近似しても流速u、v、温度Tの精度は大きく低下しない。
Here, an example of a method for deriving the numbers M u , M v , and M T of the basis vectors φ u j (l), φ v j (l), and φ T j (l) of the flow velocities u, v , and temperature T used in the calculations of equations ( 4) , (5), and (6) will be described.
When the basis vectors φuj ( l ), φvj (l), and φTj (l) of the flow velocities u , v, and temperature T are arranged in order from largest singular value to smallest, their magnitudes decrease exponentially. Therefore, compared with the basis vectors φuj (l), φvj (l), and φTj (l ) of the low-order modes, the contribution of the basis vectors φuj (l), φvj ( l), and φTj (l ) of the high - order modes to the flow velocities u, v, and temperature T is low . Therefore, even if the flow velocities u, v, and temperature T are approximated by a linear combination of only the basis vectors φuj ( l), φvj (l), and φTj (l) of the low-order modes, the accuracy of the flow velocities u, v, and temperature T does not decrease significantly.
そこで本実施形態では、基底ベクトル導出部103は、流速u、v、温度Tのスナップショット行列MU、MV、MTを特異値分解することにより得られる特異値を、モード番号jの小さい特異値(即ち大きな値の特異値)から順に積算する。そして、基底ベクトル導出部103は、積算値が所定の条件を最初に満足するときのモード番号jの最大値を、(4)式、(5)式、(6)式の計算に用いる流速u、v、温度Tの基底ベクトルφu
j(l)、φv
j(l)、φT
j(l)の数Mu、Mv、MTとする。流速u、v、温度Tのスナップショット行列MU、MV、MTを特異値分解することにより得られる特異値はそれぞれ異なる。従って、(4)式、(5)式、(6)式の計算に用いる流速u、v、温度Tの基底ベクトルφu
j(l)、φv
j(l)、φT
j(l)の数Mu、Mv、MTも異なる値になり得る。
Therefore, in this embodiment, the basis
所定の条件としては、例えば、全ての特異値に対する割合が所定値以上であることを採用することができる。所定値としては、例えば、99.999%を採用することができる。このようにする場合、流速uのスナップショット行列MUを特異値分解することにより得られるモード番号1~jの特異値の積算値を、流速uのスナップショット行列MUを特異値分解することにより得られる全ての特異値の積算値で割った値が0.99999以上になるか否かを判定することを、1を初期値としてモード番号jを1つずつ増加させながら実行する。そして、流速uのスナップショット行列MUを特異値分解することにより得られるモード番号1~jの特異値の積算値を、流速uのスナップショット行列MUを特異値分解することにより得られる全ての特異値の積算値で割った値が0.99999以上になったと最初に判定されたときのモード番号jをMuとする。Mv、MTについても、以上の説明において、流速uのスナップショット行列MUを、流速vのスナップショット行列MV、温度Tのスナップショット行列MTに置き換えることにより導出される。
As the predetermined condition, for example, a ratio to all singular values may be equal to or greater than a predetermined value. As the predetermined value, for example, 99.999% may be used. In this case, a determination is made as to whether or not a value obtained by dividing the integrated value of the singular values of
(縮約シミュレータ導出部104)
縮約シミュレータ導出部104は、溶鋼3の流動を記述する運動方程式を、軸成分ごとに、流速u、vの基底ベクトルφu
j(l)、φv
j(l)に射影することで、係数au
j(t)、av
j(t)についての時間発展方程式(微分方程式)を導出する。また、縮約シミュレータ導出部104は、溶鋼3の熱伝導を記述する熱伝導方程式を、温度Tの基底ベクトルφT
j(l)に射影することで、係数au
j(t)、av
j(t)、aT
j(t)についての時間発展方程式(微分方程式)を導出する。
(Contraction simulator derivation unit 104)
The reduced
例えば、溶鋼3のl1軸方向(幅方向)成分の流動を記述する運動方程式、溶鋼3のl2軸方向(鋳造方向)成分の流動を記述する運動方程式、溶鋼3の熱伝導を記述する熱伝導方程式は、それぞれ、以下の(7)式、(8)式、(9)式で表される。
For example, an equation of motion describing the flow of the l1 axial direction (width direction) component of the
ここで、pは圧力であり、ρは密度であり、νは粘性係数であり、kは熱伝導率である。尚、運動方程式は、数値シミュレーション部101で用いるナビエストークス方程式に対応し、熱伝導方程式は、数値シミュレーション部101で用いるエネルギー保存式に対応する。
Here, p is pressure, ρ is density, ν is viscosity coefficient, and k is thermal conductivity. The equation of motion corresponds to the Navier-Stokes equation used in the
(7)式に、(4)式および(5)式を代入した式の両辺に対し、i次のモードの流速uの基底ベクトルφu i(l)との内積をとり、時間微分の項が左辺に来るように整理すると、以下の(10)式が得られる。また、(8)式に、(4)式および(5)式を代入した式の両辺に対し、i次のモードの流速vの基底ベクトルφv i(l)との内積をとり、時間微分の項が左辺に来るように整理すると、以下の(11)式が得られる。また、(9)式に、(4)式~(6)式を代入した式の両辺に対し、i次のモードの温度Tの基底ベクトルφT i(l)との内積をとり、時間微分の項が左辺に来るように整理すると、以下の(12)式が得られる。 Substituting equations (4) and (5) into equation (7), we take the inner product of both sides of the equation with the basis vector φ u i (l) of the flow velocity u in the i-th mode, and rearrange the equation so that the time differential term is on the left side to obtain the following equation (10). Substituting equations (4) and (5) into equation (8), we take the inner product of both sides of the equation with the basis vector φ v i (l) of the flow velocity v in the i-th mode, and rearrange the equation so that the time differential term is on the left side to obtain the following equation (11). Substituting equations (4) to (6) into equation (9), we take the inner product of both sides of the equation with the basis vector φ T i (l) of the temperature T in the i-th mode, and rearrange the equation so that the time differential term is on the left side to obtain the following equation (12).
尚、流速uの基底ベクトルφu j(l)、流速vの基底ベクトルφv j(l)、温度Tの基底ベクトルφT j(l)は、それぞれ正規直交系をなす。従って、<φu m(l),φu n(l)>=δmn、<φv m(l),φv n(l)>=δmn、<φT m(l),φT n(l)>=δmnが成り立つ。ここで、<a,b>は、aとbとの内積を示す。また、δmnは、クロネッカーのデルタ(m=nのときδmn=1、m≠nのときδmn=0)である。 In addition, the basis vector φuj (l) of flow velocity u, the basis vector φvj (l ) of flow velocity v, and the basis vector φTj (l) of temperature T each form an orthonormal system. Therefore, < φu (l), φu (l)> = δmn , < φvm (l), φvn (l)> = δmn , and < φTm ( l), φTn ( l)> = δmn hold. Here, <a, b> indicates the inner product of a and b. Also, δmn is the Kronecker delta ( δmn = 1 when m = n, and δmn = 0 when m ≠ n).
尚、(10)式において、Au iのuとiとは、Au iが流速uのi次のモードに対する係数au i(t)の時間微分を与える式に用いられる定数であることを示す。Buu ijのuuの1つ目のuとijのiとは、Buu ijが流速uのi次のモードに対する係数au i(t)の時間微分を与える式に用いられる定数であることを示し、uuの2つ目のuとijのjとは、Buu ijが流速uのj次のモードに対する係数au j(t)に乗算される定数であることを示す。Buv ijのuvのuとijのiとは、Buv ijが流速uのi次のモードに対する係数au i(t)の時間微分を与える式に用いられる定数であることを示し、uvのvとijのjとは、Buv ijが流速vのj次のモードに対する係数av j(t)に乗算される定数であることを示す。Cuu ijkのuuの1つ目のuとijkのiとは、Cuu ijkが流速uのi次のモードに対する係数au i(t)の時間微分を与える式に用いられる定数であることを示し、ijkのjは、Cuv ijkが流速uのj次のモードに対する係数au j(t)に乗算される定数であることを示し、uuの2つ目のuとijkのkとは、Cuu ijkが流速uのk次のモードに対する係数au k(t)に乗算される定数であることを示す。Cuv ijkのuvのuとijkのiとは、Cuv ijkが流速uのi次のモードに対する係数au i(t)の時間微分を与える式に用いられる定数であることを示し、ijkのjは、Cuv ijkが流速uのj次のモードに対する係数au j(t)に乗算される定数であることを示し、uvのvとijkのkとは、Cuv ijkが流速vのk次のモードに対する係数av k(t)に乗算される定数であることを示す。 In addition, in formula (10), u and i of A u i indicate that A u i is a constant used in the formula that gives the time derivative of the coefficient a u i (t) for the i-th mode of the flow velocity u. The first u of uu and i of ij of B uu ij indicate that B uu ij is a constant used in the formula that gives the time derivative of the coefficient a u i (t) for the i-th mode of the flow velocity u, and the second u of uu and j of ij indicate that B uu ij is a constant by which the coefficient a u j (t) for the j-th mode of the flow velocity u is multiplied. The u of uv and i of ij of B uv ij indicate that B uv ij is a constant used in the formula that gives the time derivative of the coefficient a u i (t) for the i-th mode of the flow velocity u, and the v of uv and j of ij indicate that B uv ij is a constant by which the coefficient a v j (t) for the j-th mode of the flow velocity v is multiplied. The first u in uu of C uu ijk and the i in ijk indicate that C uu ijk is a constant used in the equation giving the time derivative of the coefficient a u i (t) for the i-th mode of flow velocity u, the j in ijk indicates that C uv ijk is a constant by which the coefficient a u j (t) for the j-th mode of flow velocity u is multiplied, and the second u in uu and the k in ijk indicate that C uu ijk is a constant by which the coefficient a u k (t) for the k-th mode of flow velocity u is multiplied. The u in uv and the i in ijk of C uv ijk indicate that C uv ijk is a constant used in an equation giving the time derivative of the coefficient a u i (t) for the i-th mode of flow velocity u, the j in ijk indicates that C uv ijk is a constant by which the coefficient a u j (t) for the j-th mode of flow velocity u is multiplied, and the v in uv and the k in ijk indicate that C uv ijk is a constant by which the coefficient a v k (t) for the k-th mode of flow velocity v is multiplied.
また、(11)式において、Av iのvとiとは、Av iが流速vのi次のモードに対する係数av i(t)の時間微分を与える式に用いられる定数であることを示す。Bvu ijのvuのvとijのiとは、Bvu ijが流速vのi次のモードに対する係数av i(t)の時間微分を与える式に用いられる定数であることを示し、vuのuとijのjとは、Bvu ijが流速uのj次のモードに対する係数au j(t)に乗算される定数であることを示す。Bvv ijのvvの1つ目のvとijのiとは、Bvv ijが流速vのi次のモードに対する係数av i(t)の時間微分を与える式に用いられる定数であることを示し、vvの2つ目のvとijのjとは、Bvv ijが流速vのj次のモードに対する係数av j(t)に乗算される定数であることを示す。Cvu ijkのvuのvとijkのiとは、Cvu ijkが流速vのi次のモードに対する係数av i(t)の時間微分を与える式に用いられる定数であることを示し、ijkのjは、Cvu ijkが流速vのj次のモードに対する係数av j(t)に乗算される定数であることを示し、vuのuとijkのkとは、Cvu ijkが流速uのk次のモードに対する係数au k(t)に乗算される定数であることを示す。Cvv ijkのvvの1つ目のvとijkのiとは、Cvv ijkが流速vのi次のモードに対する係数av i(t)の時間微分を与える式に用いられる定数であることを示し、ijkのjは、Cvu ijkが流速vのj次のモードに対する係数av j(t)に乗算される定数であることを示し、vvの2つ目のvとijkのkとは、Cvv ijkが流速vのk次のモードに対する係数av k(t)に乗算される定数であることを示す。 In addition, in formula (11), the v and i of A v i indicate that A v i is a constant used in a formula that gives the time derivative of the coefficient a v i (t) for the i-th mode of the flow velocity v. The v of vu and the i of ij of B vu ij indicate that B vu ij is a constant used in a formula that gives the time derivative of the coefficient a v i (t) for the i-th mode of the flow velocity v, and the u of vu and the j of ij indicate that B vu ij is a constant by which the coefficient a u j (t) for the j-th mode of the flow velocity u is multiplied. The first v of vv and the i of ij of B vv ij indicate that B vv ij is a constant used in a formula that gives the time derivative of the coefficient a v i (t) for the i-th mode of the flow velocity v, and the second v of vv and the j of ij indicate that B vv ij is a constant by which the coefficient a v j (t) for the j-th mode of the flow velocity v is multiplied. The v in vu and the i in ijk of Cvu ijk indicate that Cvu ijk is a constant used in an equation giving the time derivative of the coefficient avi (t) for the i-th mode of flow velocity v , the j in ijk indicates that Cvu ijk is a constant by which the coefficient auvj (t) for the j-th mode of flow velocity v is multiplied, and the u in vu and the k in ijk indicate that Cvu ijk is a constant by which the coefficient auk ( t ) for the k-th mode of flow velocity u is multiplied. The first v in vv and the i in ijk of C vv ijk indicate that C vv ijk is a constant used in an equation giving the time derivative of the coefficient a v i (t) for the i-th mode of flow velocity v, the j in ijk indicates that C vu ijk is a constant by which the coefficient a v j (t) for the j-th mode of flow velocity v is multiplied, and the second v in vv and the k in ijk indicate that C vv ijk is a constant by which the coefficient a v k (t) for the k-th mode of flow velocity v is multiplied.
また、(12)式において、AT iのTとiとは、AT iが温度Tのi次のモードに対する係数aT i(t)の時間微分を与える式に用いられる定数であることを示す。BTT ijのTTの1つ目のTとijのiとは、BTT ijが温度Tのi次のモードに対する係数aT i(t)の時間微分を与える式に用いられる定数であることを示し、TTの2つ目のTとijのjとは、BTT ijが温度Tのj次のモードに対する係数aT j(t)に乗算される定数であることを示す。BTu ijのTuのTとijのiとは、BTu ijが温度Tのi次のモードに対する係数aT i(t)の時間微分を与える式に用いられる定数であることを示し、Tuのuとijのjとは、BTu ijが流速uのj次のモードに対する係数au j(t)に乗算される定数であることを示す。BTv ijのTvのTとijのiとは、BTv ijが温度Tのi次のモードに対する係数aT i(t)の時間微分を与える式に用いられる定数であることを示し、Tvのvとijのjとは、BTv ijが流速vのj次のモードに対する係数av j(t)に乗算される定数であることを示す。CTu ijkのTuのTとijkのiとは、CTu ijkが温度Tのi次のモードに対する係数aT i(t)の時間微分を与える式に用いられる定数であることを示し、Tuのuとijkのjとは、CTu ijkが流速uのj次のモードに対する係数au j(t)に乗算される定数であることを示し、ijkのkは、CTu ijkが温度Tのk次のモードに対する係数aT k(t)に乗算される定数であることを示す。CTv ijkのTvのTとijkのiとは、CTv ijkが温度Tのi次のモードに対する係数aT i(t)の時間微分を与える式に用いられる定数であることを示し、Tvのvとijkのjとは、CTv ijkが流速vのj次のモードに対する係数av j(t)に乗算される定数であることを示し、ijkのkは、CTv ijkが温度Tのk次のモードに対する係数aT k(t)に乗算される定数であることを示す。 In addition, in equation (12), the T and i in A T i indicate that A T i is a constant used in an equation that gives the time derivative of the coefficient a T i (t) for the i-th mode of temperature T. The first T in TT and the i in ij in B TT ij indicate that B TT ij is a constant used in an equation that gives the time derivative of the coefficient a T i (t) for the i-th mode of temperature T, and the second T in TT and the j in ij indicate that B TT ij is a constant by which the coefficient a T j (t) for the j-th mode of temperature T is multiplied. The T and i in Tu of B Tu ij indicate that B Tu ij is a constant used in an equation that gives the time derivative of the coefficient a T i (t) for the i-th mode of temperature T, and the u in Tu and the j in ij indicate that B Tu ij is a constant by which the coefficient a u j (t) for the j-th mode of flow velocity u is multiplied. The T in Tv and the i in ij of B Tv ij indicate that B Tv ij is a constant used in an equation that gives the time derivative of the coefficient a T i (t) for the i-th mode of temperature T, and the v in Tv and the j in ij indicate that B Tv ij is a constant by which the coefficient a v j (t) for the j-th mode of flow velocity v is multiplied. The T in Tu and the i in ijk of C Tu ijk indicate that C Tu ijk is a constant used in an equation that gives the time derivative of the coefficient a T i (t) for the i-th mode of temperature T, the u in Tu and the j in ijk indicate that C Tu ijk is a constant by which the coefficient a u j (t) for the j-th mode of flow velocity u is multiplied, and the k in ijk indicates that C Tu ijk is a constant by which the coefficient a T k (t) for the k-th mode of temperature T is multiplied. The T in Tv and the i in ijk of C Tv ijk indicate that C Tv ijk is a constant used in the equation that gives the time derivative of the coefficient a T i (t) for the i-th mode of temperature T, the v in Tv and the j in ijk indicate that C Tv ijk is a constant by which the coefficient a v j (t) for the j-th mode of flow velocity v is multiplied, and the k in ijk indicates that C Tv ijk is a constant by which the coefficient a T k (t) for the k-th mode of temperature T is multiplied.
以上のように本実施形態では、縮約シミュレータ導出部104は、基底ベクトル導出部103により導出された、流速u、v、温度Tの基底ベクトルφu
1(l)~φu
Mu(l)、φv
1(l)~φv
Mv(l)、φT
1(l)~φT
MT(l)および基底ベクトルの数Mu、Mv、MTと、(4)式~(9)式を用いて、(10)式の定数Ai
u、Buu
ij、Buv
ij、Cuu
ijk、Cuv
ijkと、(11)式の定数Ai
v、Bvu
ij、Bvv
ij、Cvu
ijk、Cvv
ijkと、(12)式の定数Ai
T、BTT
ij、BTu
ij、BTv
ij、CTu
ijk、CTv
ijkと、を導出する。本実施形態では、このようにして導出された定数が設定された(10)式~(12)式が縮約シミュレータとなる。縮約シミュレータは、非特許文献1に記載されている縮約モデルに相当するものである。支配方程式(基礎方程式)を基底に射影して縮約モデルを導出すること自体は、非特許部文献1、2に記載されているので、ここでは、その詳細な説明を省略する。
As described above, in this embodiment, the reduced
(状態導出部111)
状態導出部111は、溶鋼3の状態を示す量である状態量を成分として含む状態ベクトルxt導出する。本実施形態では、状態導出部111は、後述するデータ同化部115によりデータ同化を行う際の状態ベクトルxtに含める状態量を導出する。また、本実施形態では、状態量が、第1の物理量の一例である溶鋼3の流速u、vおよび温度Tである場合を例に挙げて説明する。状態導出部111は、第1の係数導出部111aと、物理量導出部111bと、を有する。
(State Derivation Unit 111)
The
((第1の係数導出部111a))
第1の係数導出部111aは、時刻tにおいて、時刻t-Δtでの計算格子点lにおける溶鋼3の流速と、時刻t-Δtでの計算格子点lにおける溶鋼3の温度と、時刻t-Δtから時刻tまでの間での境界条件のパラメータ(例えば、流速、圧力、熱流束)の少なくとも何れか1つとのうち、少なくとも何れか1つが異なる複数のケースのそれぞれについて、(10)式~(12)式の微分方程式の解(au
1(t)~au
Mu(t)、av
1(t)~av
Mv(t)、aT
1(t)~aT
MT(t))を導出する。
((First
The first
第1の係数導出部111aは、最初の時刻t(=t0)での溶鋼3の流速と、最初の時刻t(=t0)での溶鋼3の温度と、時刻t0から時刻t0+Δtまでの間での境界条件のパラメータの少なくとも1つとのうち、少なくとも何れか1つが異なる複数のケースに対して、初期値として、予め設定された値を用いる。2番目の時刻t以降に関しては、後述する観測データ取得部113で取得された観測データ(情報)を基に後述するデータ同化部115で導出された結果に基づいて、時刻tでの溶鋼3の流速および温度、並びに時刻tから時刻t+Δtまでの間での境界条件のパラメータを導出する。この点の詳細については後述する。
The first
((物理量導出部111b))
物理量導出部111bは、第1の係数導出部111aにより導出された、流速uに対する係数au
1(t)~au
Mu(t)と、基底ベクトル導出部103により導出された、流速uの基底ベクトルφu
1(l)~φu
Mu(l)とを(4)式に代入することにより、時刻tでの各計算格子点lにおける溶鋼3の流速uを導出する。また、物理量導出部111bは、第1の係数導出部111aにより導出された、流速vに対する係数av
1(t)~av
Mv(t)と、基底ベクトル導出部103により導出された、流速vの基底ベクトルφv
1(l)~φv
Mv(l)とを(5)式に代入することにより、時刻tでの各計算格子点lにおける溶鋼3の流速vを導出する。また、物理量導出部111bは、第1の係数導出部111aにより導出された温度Tに対する係数aT
1(t)~aT
MT(t)と、基底ベクトル導出部103により導出された温度Tの基底ベクトルφT
1(l)~φT
MT(l)とを(6)式に代入することにより、時刻tでの各計算格子点lにおける溶鋼3の温度Tを導出する。物理量導出部111bは、前述した複数のケースの係数au
1(t)~au
Mu(t)、av
1(t)~av
Mv(t)、aT
1(t)~aT
MT(t)のそれぞれについて、時刻tでの各計算格子点lにおける溶鋼3の流速u、v、温度Tを導出する。
((Physical
The physical
本実施形態の以下の説明では、計算格子点lにおける溶鋼3の流速u、vおよび温度Tの全てで構成されるベクトルを必要に応じて状態ベクトルと称し、物理量導出部111bにより導出される溶鋼3の流速u、vおよび温度Tを必要に応じて数値解析データと称する。物理量導出部111bにより、前述した複数のケースのそれぞれについて数値解析データが導出される。
In the following description of this embodiment, the vector composed of all of the flow velocities u, v and temperature T of the
(第1の確率密度関数導出部112)
第1の確率密度関数導出部112は、後述する観測データ取得部113で取得された時刻t-Δtまでの観測データ(情報)yt0:t-Δtを基に推定される時刻tでの状態ベクトルxtの条件付き確率密度関数p(xt|yt0:t-Δt)を導出する。確率密度関数p(xt|yt0:t-Δt)は、物理量導出部111bにより導出された複数の数値解析データ{xt
(k)}kが示す確率密度関数である。以下の説明では、この確率密度関数p(xt|yt0:t-Δt)を必要に応じて第1の確率密度関数とも称する。本実施形態では、第1の確率密度関数導出部112は、ステップ幅Δtの時間隔で、各時刻tでの状態ベクトルxtの第1の確率密度関数p(xt|yt0:t-Δt)を導出する。
(First probability density function derivation unit 112)
The first probability density
第1の確率密度関数導出部112は、物理量導出部111bにより導出された時刻tでの複数の数値解析データ{xt|t-Δt
(k)}kの相対度数(度数を全データ数で除した値)を用いた度数分布を導出する。この度数分布が、時刻tでの状態ベクトルxtの第1の確率密度関数p(xt|yt0:t-Δt)となる。また、複数の数値解析データ{xt|t-Δt
(k)}kの各成分(各計算格子点lにおける各第1の物理量)の度数分布を導出して、これらの積を取ったものを複数の数値解析データ{xt|t-Δt
(k)}kの度数分布としてもよい。また、これらの度数分布を所定の関数(例えば、スプライン関数)で近似したものを第1の確率密度関数としてもよい。
The first probability density
また、第1の確率密度関数導出部112は、物理量導出部111bにより導出された時刻tでの複数の数値解析データ{xt|t-Δt
(k)}kを用いてアンサンブル近似(モンテカルロ近似)を行うことで、時刻tでの状態ベクトルxtの確率密度関数p(xt|yt0:t-Δt)を導出することもできる。
また、第1の確率密度関数導出部112は、時刻tでの状態ベクトルxtの第1の確率密度関数p(xt|yt0:t-Δt)がガウス分布であると仮定して、物理量導出部111bにより導出された時刻tでの複数の数値解析データ{xt
(k)}kから平均値と分散値を導出してもよい。
尚、第1の確率密度関数導出部112は、以上のようにして導出される時刻tでの状態ベクトルxtの第1の確率密度関数p(xt|yt0:t-Δt)を、ガウス分布のような確率密度関数をもつノイズを用いて修正してもよい。
In addition, the first probability density
In addition, the first probability density
Furthermore, the first probability density
(観測データ取得部113)
観測データ取得部113は、連続鋳造設備において測定された時系列データ(測定値)に基づいて観測データを所定の周期で取得する。
本実施形態では、観測データ取得部113は、温度計F、Lで測定される温度(鋳型4内の温度)と、湯面レベル計9で測定される湯面レベルとを含む測定値を入力する。
(Observation data acquisition unit 113)
The observation data acquiring unit 113 acquires observation data at a predetermined cycle based on time-series data (measured values) measured in the continuous casting equipment.
In this embodiment, the observation data acquisition unit 113 inputs measurement values including the temperatures measured by the thermometers F and L (the temperatures inside the mold 4 ) and the molten metal level measured by the molten metal level gauge 9 .
温度計F、Lで測定される温度は、鋳型4の温度であり、溶鋼3の温度ではない。そこで、本実施形態では、観測データ取得部113は、鋳型4の奥行方向(l3軸方向)において相互に対向する2つの温度計F、Lの合計12組のそれぞれについて、当該2つの温度計F、Lで同時刻に測定された温度の算術平均値を導出する。尚、2つの温度計F、Lの組は、温度計F1、L1、温度計F2、L2などである。これにより12個の温度が得られる。これらの温度は、当該温度の導出元となる2つの温度計F、Lが存在しているl1-l2平面上の位置(座標(l1,l2))の温度となる。
The temperatures measured by the thermometers F and L are the temperatures of the
観測データ取得部113は、鋳型4の長辺における伝熱状態が奥行方向(l3軸方向)の1次元定常熱伝導であるとの仮定をおいて、鋳型4の冷却水の温度と、鋳型4の熱伝導度とに基づいて、これら12個の温度を外挿し、可視化対象面の各位置での溶鋼3の温度を導出する。
また、観測データ取得部113は、湯面レベル計9で測定された複数の湯面レベルの単位時間当たりの変化から、可視化対象面の溶鋼3の湯面に対応する各位置での鋳造方向(l2軸方向)の流速を導出する。
The observation data acquisition unit 113 assumes that the heat transfer state on the long side of the
In addition, the observation data acquisition unit 113 derives the flow velocity in the casting direction (l2 axial direction) at each position corresponding to the
以上のように本実施形態では、観測データ取得部113は、可視化の対象となる2次元平面の各位置での溶鋼3の温度と、可視化の対象となる2次元平面の溶鋼3の湯面に対応する各位置での鋳造方向の流速とを観測データとして取得する。本実施形態における以下の説明では、当該位置を必要に応じて観測位置と称し、観測位置において第2の物理量の一例として取得される溶鋼3の温度および鋳造方向の流速の全てで構成されるベクトルを観測ベクトルと称する。尚、本実施形態では、観測データ取得部113は、連続鋳造設備において測定された時系列データから、物理量導出部111bで数値シミュレーションを行う時刻に対応する時刻のデータを抽出し、抽出したデータを用いて、可視化の対象となる2次元平面の各観測位置での溶鋼3の温度と、可視化の対象となる2次元平面の溶鋼3の湯面に対応する各観測位置での鋳造方向の流速とを観測データとして取得する。
As described above, in this embodiment, the observation data acquisition unit 113 acquires, as observation data, the temperature of the
(尤度関数導出部114)
尤度関数導出部114は、観測データ取得部113で取得された時刻tでの観測データyt
(0)(yt
(0)は観測位置において取得される溶鋼3の温度および鋳造方向の流速のデータの全てで構成されるベクトル)の観測ノイズを、過去の当該観測データから評価することにより、時刻tでの状態ベクトルxtが得られたときの観測ベクトルytの条件付き確率密度関数p(yt|xt)に対する状態ベクトルxtの尤度関数L(xt|yt)を導出する。尤度関数L(xt|yt)とは、条件付き確率密度関数p(yt|xt)を状態ベクトルxtの関数と見做したもので、観測ベクトルytに観測データyt
(0)を代入することによって状態ベクトルxtの関数として確定する。本実施形態では、尤度関数導出部114は、ステップ幅Δtの時間隔で、各時刻tでの状態ベクトルxtの尤度関数L(xt|yt)を導出する。
(Likelihood function derivation unit 114)
The likelihood
観測ノイズは、例えば、平均ベクトルが0(ゼロベクトル)のガウス分布に従うものとし、観測データの分散共分散行列を観測ノイズの共分散行列として用いることができる。この場合、過去の同種の観測データの分散共分散行列値を、全ての種類の観測データについて予め計算しておき、尤度関数導出部114が記憶しておく。即ち、尤度関数導出部114は、観測データの種類毎の分散共分散行列値を記憶する。また、観測データの分散共分散行列値は、観測データの各成分(各観測位置での値)の分散値とし、共分散値は0(ゼロ)としたものでもよい。
The observation noise may, for example, follow a Gaussian distribution with a mean vector of 0 (zero vector), and the variance-covariance matrix of the observation data may be used as the covariance matrix of the observation noise. In this case, the variance-covariance matrix values of past observation data of the same type are calculated in advance for all types of observation data, and stored by the likelihood
尤度関数導出部114は、時刻tでの観測データyt
(0)と、観測ノイズの分散共分散行列値とに基づいて、時刻tでの状態ベクトルxtの尤度関数L(xt|yt)を導出する。
The likelihood
観測ノイズは、ガウス分布に従うとすることに限定されない。例えば、尤度関数導出部114は、時刻tでの観測データyt
(0)と、当該時刻tよりも前の所定の期間の当該観測データyt
(0)とに基づいて、これらの複数の観測データ{yt-nΔt
(0),・・・,yt
(0)}の平均ベクトルm(yt
(0))を導出し,複数のデータ{yt-nΔt
(0)-m(yt
(0)),・・・,yt
(0)-m(yt
(0))}の相対度数(度数を全データ数で除した値)を用いた度数分布を導出する。そして、尤度関数導出部114は、当該度数分布を、時刻tでの観測ノイズの確率密度関数として導出する。また、複数のデータ{yt-nΔt
(0)-m(yt
(0)),・・・,yt
(0)-m(yt
(0))}の各成分(各観測位置αにおける各データ)の度数分布を導出して、これらの積を取ったものを複数のデータ{yt-nΔt
(0)-m(yt
(0)),・・・,yt
(0)-m(yt
(0))}の度数分布としてもよい。また、尤度関数導出部114は、この度数分布を所定の関数(例えば、スプライン関数)で近似したものを観測ノイズの確率密度関数としてもよい。
The observation noise is not limited to being Gaussian distribution. For example, the likelihood
(データ同化部115)
データ同化部115は、第1の確率密度関数導出部112により導出された時刻tでの状態ベクトルxtの第1の確率密度関数p(xt|y0:t-Δt)と、尤度関数導出部114により導出された時刻tでの観測データが得られたときの当該時刻tでの状態ベクトルxtの尤度関数L(xt|yt)とを、ベイズの定理を基礎としたベイズ統計のモデリングによるデータ同化を行うフィルタのアルゴリズムに与えることにより、時刻tまでの観測データが得られたときの時刻tでの状態ベクトルxtの条件付き確率密度関数p(xt|yt0:t)を導出する。以下の説明では、この確率密度関数p(xt|yt0:t)を必要に応じて第2の確率密度関数とも称する。本実施形態では、データ同化部115は、ステップ幅Δtの時間隔で、各時刻tでの状態ベクトルxtの第2の確率密度関数p(xt|y0:t)を導出する。
(Data Assimilation Unit 115)
The
データ同化を行う際には、システム方程式(状態方程式)と観測方程式とを定める必要がある。システム方程式は、ステップ幅Δtの時間隔で、前後の時刻の状態ベクトルxt+Δt、xtの関係を定義する式である。本実施形態では、物理量導出部111bによる複数の数値シミュレーションでこの関係は決定される。また、観測方程式は、観測ベクトルyt(観測量)と、状態ベクトルxt(状態量)との関係式を示すものであり、(13)式で与えられる。
yt=Htxt+wt ・・・(13)
ここで、Htは、観測行列であり、wtは、観測ノイズである。例えば、観測ベクトルの第α成分が状態ベクトルの第β成分の物理量と同一であれば、観測行列Htは、その(α,β)成分が1で、第α行の他の成分が0(ゼロ)である行列となる。尚、ここでの物理量が同一とは、軸成分も含めて同一であることを指す。従って、流速u、vは異なる物理量となる。また、観測位置αがどの計算格子点lとも一致しなければ、状態ベクトルで補間するように観測行列Htを定めることもできる。
When performing data assimilation, it is necessary to determine a system equation (state equation) and an observation equation. The system equation is an equation that defines the relationship between state vectors xt + Δt and xt at previous and next times at a time interval of a step width Δt. In this embodiment, this relationship is determined by multiple numerical simulations by the physical
yt = Htxt + wt ... (13)
Here, Ht is the observation matrix, and wt is the observation noise. For example, if the α-th component of the observation vector is the same as the physical quantity of the β-th component of the state vector, the observation matrix Ht will be a matrix whose (α, β) component is 1 and the other components of the α-th row are 0 (zero). Note that the physical quantities being the same here refer to the same quantities including the axial components. Therefore, the flow velocities u and v are different physical quantities. Also, if the observation position α does not coincide with any of the computational grid points l, the observation matrix Ht can be determined so as to be interpolated with the state vector.
ベイズ統計のモデリングによるデータ同化を行うフィルタとしては、例えば、アンサンブルカルマンフィルタを用いることができる。
この場合、データ同化部115は、第1の確率密度関数導出部112により導出された時刻tでの状態ベクトルxtの第1の確率密度関数p(xt|y0:t-Δt)と、尤度関数導出部114により導出された当該時刻tでの状態ベクトルxtの尤度関数L(xt|yt)とを、アンサンブル近似し、アンサンブルメンバー(粒子)を導出する。時刻t-Δtの各アンサンブルメンバーxt-Δt
(k)を、システム方程式に基づいて更新し、一期先(時刻t)の予測分布のアンサンブルメンバーxt|t-Δt
(k)を第1の確率密度関数導出部112により導出する。そして、データ同化部115は、一期先の予測分布のアンサンブルメンバーxt|t-Δt
(k)から、フィルタ分布のアンサンブルメンバーxt
(k)を導出する。このようなフィルタ分布のアンサンブルメンバーの導出を行うことにより、時刻tでの状態ベクトルxtの第2の確率密度関数が導出される。
As a filter that performs data assimilation through Bayesian statistical modeling, for example, an ensemble Kalman filter can be used.
In this case, the
連続鋳造設備においては、全ての位置における溶鋼3の流速を測定することはできない。しかしながら、システム方程式において、溶鋼3の温度と流速はカップリングされるので、状態ベクトルの確率密度関数は、溶鋼3の流速の確率密度関数を包含している。本実施形態では、溶鋼3の湯面での鋳造方向の流速については測定しているので、溶鋼3の湯面での鋳造方向の流速についてもデータ同化を行うことにより、溶鋼3の流速の確率密度関数の精度を向上させることができる。
In continuous casting equipment, it is not possible to measure the flow velocity of
ベイズ統計のモデリングによるデータ同化を行うフィルタは、アンサンブルカルマンフィルタに限定されない。例えば、粒子フィルタを用いてもよい。アンサンブルカルマンフィルタや粒子フィルタについては、非特許文献4、5に記載されている。また、アンサンブルカルマンフィルタや粒子フィルタの具体的な計算アルゴリズムについては、非特許文献6に記載されている。従って、これらについての詳細な説明を省略する。
The filter that performs data assimilation using Bayesian statistical modeling is not limited to the ensemble Kalman filter. For example, a particle filter may be used. The ensemble Kalman filter and particle filter are described in
データ同化部115は、時刻tでの状態ベクトルxtの第2の確率密度関数p(xt|yt)の最頻値xmを、時刻tでの状態ベクトルxtの推定値として導出する。また、状態ベクトルxtの各成分(各計算格子点lにおける各第1の物理量)について、第2の確率密度関数p(xt|yt)の周辺確率密度関数を導出し、周辺確率密度関数の各成分の最頻値で構成されるベクトルを第2の確率密度関数p(xt|yt)の最頻値xmとしてもよい。
以上のように本実施形態では、データ同化部115は、第1の確率密度関数導出部112により導出された時刻tでの状態ベクトルxtの第1の確率密度関数p(xt|y0:t-Δt)と、尤度関数導出部114により導出された当該時刻tでの状態ベクトルxtの尤度関数L(xt|yt)とを融合させて、時刻tでの最も合理的な状態ベクトルxtの第2の確率密度関数p(xt|yt)を導出し、その最頻値xmを、時刻tでの状態ベクトルxtの推定値として導出する。
The
As described above, in this embodiment, the
また、データ同化部115は、状態ベクトルxtの第2の確率密度関数p(xt|yt)を第1の係数導出部111aに出力する。第1の係数導出部111aは、状態ベクトルxtの第2の確率密度関数p(xt|yt)のアンサンブルメンバーに基づいて、時刻tでの各計算格子点lにおける溶鋼3の流速u(l,t)、v(l,t)および温度T(l,t)と、時刻tから時刻t+Δtまでの間での境界条件のパラメータとの組を複数導出する。第1の係数導出部111aは、時刻tでの各計算格子点lにおける溶鋼3の流速u(l,t)、v(l,t)および温度T(l,t)と、流速u、vの基底ベクトルφu
1(l)~φu
Mu(l)、φv
1(l)~φv
Mv(l)および温度Tの基底ベクトルφT
1(l)~φT
MT(l)と、に基づいて、以下の(14)式~(16)式により、時刻tでの係数au
1(t)~au
Mu(t)、av
1(t)~av
Mv(t)、aT
1(t)~aT
MT(t)を導出する。尚、流速u、vの基底ベクトルφu
1(l)~φu
Mu(l)、φv
1(l)~φv
Mv(l)および温度Tの基底ベクトルφT
1(l)~φT
MT(l)は、基底ベクトル導出部103により導出されたものである。また、第1の係数導出部111aは、このようにして導出した時刻tでの係数au
1(t)~au
Mu(t)、av
1(t)~av
Mv(t)、aT
1(t)~aT
MT(t)と、時刻tから時刻t+Δtまでの間での境界条件のパラメータと、に基づいて、(10)式~(12)式により、係数au
1(t+Δt)~au
Mu(t+Δt)、av
1(t+Δt)~av
Mv(t+Δt)、aT
1(t+Δt)~aT
MT(t+Δt)を導出する。係数au
1(t+Δt)、av
1(t+Δt)、aT
1(t+Δt)の導出は、前述した複数のケースのそれぞれについて実行される。従って、前述した複数のケースについての係数au
1(t+Δt)~au
Mu(t+Δt)、av
1(t+Δt)~av
Mv(t+Δt)、aT
1(t+Δt)~aT
MT(t+Δt)のアンサンブルメンバーが導出される。
Furthermore, the
(可視化データ作成部116)
可視化データ作成部116は、データ同化部115により導出された時刻tでの状態ベクトルの推定値を構成する各位置(各計算格子点l)における溶鋼3の温度Tと流速u、vの少なくとも何れか1つの表示データを作成する。
可視化データ作成部116は、各時刻での各位置における溶鋼3の温度と流速を、情報処理装置100のユーザが認識できるようにしていれば、どのような表示データを作成してもよい。ただし、可視化データ作成部116は、情報処理装置100のユーザが直感的に分かり易い表示データを作成するのが好ましい。
(Visualization Data Creation Unit 116)
The visualization data creation unit 116 creates display data for at least one of the temperature T and flow velocities u, v of the
The visualization data creation unit 116 may create any type of display data as long as it enables the user of the
例えば、可視化データ作成部116は、溶鋼3の温度として想定される温度範囲を複数の領域に分けた場合の、それぞれの領域に対し異なる表示態様を予め記憶する。表示態様としては、例えば、色、模様、および濃度の少なくとも1つを採用することができる。可視化データ作成部116は、データ同化部115により導出された時刻tでの各位置における溶鋼3の温度から、各位置における温度を特定し、可視化対象面の各位置の画像が、特定した温度に対応する表示態様で表示されるように表示データを作成する。
For example, the visualization data creation unit 116 prestores a different display mode for each region when the temperature range assumed for the temperature of the
また、可視化データ作成部116は、溶鋼3の流速として想定される流速範囲を複数の領域に分けた場合の、それぞれの領域に対し異なる表示態様を予め記憶する。表示態様としては、例えば、色、模様、および濃度の少なくとも1つを採用することができる。可視化データ作成部116は、データ同化部115により導出された時刻tでの各位置における溶鋼3の流速から、各位置における流速を特定し、可視化対象面の各位置の画像が、特定した流速に対応する表示態様で表示されるように表示データを作成する。また、可視化データ作成部116は、可視化対象面に対して予め設定された領域ごとに溶鋼3の代表的な向きを導出し、当該向きの方向を向く矢印線が、可視化対象面の当該領域の画像に重ねて表示されるように表示データを作成する。
The visualization data creation unit 116 also prestores a different display mode for each region when the flow velocity range assumed as the flow velocity of the
(出力部117)
出力部117は、可視化データ作成部116で作成された表示データをコンピュータディスプレイに表示する。また、出力部117は、このような表示データに代えてまたは加えて、データ同化部115により導出された各時刻での状態ベクトルの推定値を構成する各位置における溶鋼3の流速と温度のデータを出力することができる。出力の態様は、情報処理装置100の内部または外部の記憶媒体への送信、外部装置への送信が挙げられる。このようにすれば、各時刻での各位置における溶鋼3の流速と温度を表示するための表示データを、後で別途作成することができる。
(Output unit 117)
The
尚、第1の確率密度関数導出部112、観測データ取得部113、尤度関数導出部114、データ同化部115、可視化データ作成部116、および出力部117は、特許文献1に記載されたものと同様にして実現することができる。
The first probability density
(動作フローチャート)
次に、図3のフローチャートを参照しながら、情報処理装置100のオフラインにおける処理の一例を説明する。
まず、ステップS301において、数値シミュレーション部101は、ステップ幅Δt1の時間隔の各時刻tでの各計算格子点lにおける流速u(l,t)、v(l,t)および温度T(l,t)を導出する。
(Operation Flowchart)
Next, an example of offline processing by the
First, in step S301, the
次に、ステップS302において、行列導出部102は、ステップS301でステップ幅Δt1の時間隔の各時刻tで導出された、各計算格子点lにおける流速u(l,t)、v(l,t)、温度T(l,t)から、ステップ幅Δt1よりも長いステップ幅Δt2の時間隔の各時刻tでの流速u(l,t)、v(l,t)、温度T(l,t)を抽出する。そして、行列導出部102は、抽出した流速u(l,t)、v(l,t)、温度T(l,t)が成分として格納された流速uのスナップショット行列MU、流速vのスナップショット行列MV、温度Tのスナップショット行列MTをそれぞれ導出する。
Next, in step S302, the
次に、ステップS303において、基底ベクトル導出部103は、ステップS302で導出された、流速u、v、温度Tのスナップショット行列MU、MV、MTを特異値分解することにより得られる左特異ベクトルを、流速uの基底ベクトルφu
j(l)、流速vの基底ベクトルφv
j(l)、温度Tの基底ベクトルφT
j(l)として導出する。また、基底ベクトル導出部103は、(4)式、(5)式、(6)式の計算に用いる流速u、v、温度Tの基底ベクトルφu
j(l)、φv
j(l)、φT
j(l)の数Mu、Mv、MTを導出する。
Next, in step S303, the basis
次に、ステップS304において、縮約シミュレータ導出部104は、ステップS303で導出された、流速u、v、温度Tの基底ベクトルφu
1(l)~φu
Mu(l)、φv
1(l)~φv
Mv(l)、φT
1(l)~φT
MT(l)および基底ベクトルの数Mu、Mv、MTに基づいて、縮約シミュレータを導出して記憶する。縮約シミュレータは、(10)式~(12)式において、定数Ai
u、Buu
ij、Buv
ij、Cuu
ijk、Cuv
ijk、Ai
v、Bvu
ij、Bvv
ij、Cvu
ijk、Cvv
ijk、Ai
T、BTT
ij、BTu
ij、BTv
ij、CTu
ijk、CTv
ijk、Mu、Mv、MTが設定されたものである。ステップS304の処理が終了すると、図3のフローチャートによる処理は終了する。
Next, in step S304, the reduced
次に、図4のフローチャートを参照しながら、情報処理装置100のオンラインにおける処理の一例を説明する。
まず、ステップS401において、第1の係数導出部111aは、時刻tを初期値(t0)に設定する。
次に、ステップS402において、第1の係数導出部111aは、時刻tでの溶鋼3の流速および温度と、時刻tから時刻t+Δtまでの間での境界条件のパラメータとの組の初期値として、複数の組の初期値を設定する。
Next, an example of online processing by the
First, in step S401, the first
Next, in step S402, the first
次に、ステップS403において、第1の係数導出部111aは、時刻tを時刻t+Δtに更新して、処理は、ステップS404に進む。
次に、ステップS404において、第1の係数導出部111aは、時刻t-Δtでの各計算格子点lにおける溶鋼3の流速u(l,t-Δt)、v(l,t-Δt)、温度T(l,t-Δt)と、図3のステップS303で導出された、流速u、v、温度Tの基底ベクトルφu
1(l)~φu
Mu(l)、φv
1(l)~φv
Mv(l)、φT
1(l)~φT
MT(l)とを、(14)式、(15)式、(16)式にそれぞれ代入することにより、流速u、v、温度Tに対する係数au
1(t-Δt)~au
Mu(t-Δt)、av
1(t-Δt)~av
Mv(t-Δt)、aT
1(t-Δt)~aT
MT(t-Δt)を導出する。そして、第1の係数導出部111aは、時刻tでの(10)式~(12)式の微分方程式の解(au
1(t)~au
Mu(t)、av
1(t)~av
Mv(t)、aT
1(t)~aT
MT(t))を導出する。係数au
1(t)~au
Mu(t)、av
1(t)~av
Mv(t)、aT
1(t)~aT
MT(t)は、時刻t-Δtでの各計算格子点lにおける溶鋼3の流速u(l,t-Δt)、v(l,t-Δt)および温度T(l,t-Δt)と、時刻t-Δtから時刻tまでの間での境界条件のパラメータとの組のそれぞれについて導出される。
Next, in step S403, the first
Next, in step S404, the first
次に、ステップS405において、物理量導出部111bは、流速u、v、温度Tに対する係数au
1(t)~au
Mu(t)、av
1(t)~av
Mv(t)、aT
1(t)~aT
MT(t)と、図3のステップS303で導出された、流速u、v、温度Tの基底ベクトルφu
1(l)~φu
Mu(l)、φv
1(l)~φv
Mv(l)、φT
1(l)~φT
MT(l)とを、(4)式、(5)式、(6)式にそれぞれ代入することにより、時刻tでの各計算格子点lにおける溶鋼3の流速u(l,t)、v(l,t)、温度T(l,t)を導出する。時刻tでの各計算格子点lにおける溶鋼3の流速u(l,t)、v(l,t)、温度T(l,t)は、時刻t-Δtでの各計算格子点lにおける溶鋼3の流速u(l,t-Δt)、v(l,t-Δt)および温度T(l,t-Δt)と、時刻t-Δtから時刻tまでの間での境界条件のパラメータとの組のそれぞれについて導出される。
Next, in step S405, the physical
次に、ステップS406において、第1の確率密度関数導出部112は、時刻tでの状態ベクトルxt(各計算格子点lにおける溶鋼3の流速および温度)の第1の確率密度関数p(xt|yt-Δt)を導出する。
次に、ステップS407において、観測データ取得部113は、時刻tでの観測データyt
(0)(各観測位置αにおける溶鋼3の温度と鋳造方向(l2軸方向)の流速)を取得する。
Next, in step S406, the first probability density
Next, in step S407, the observation data acquisition unit 113 acquires observation data y t (0) (the temperature of the
次に、ステップS408において、尤度関数導出部114は、観測データyt
(0)に基づいて、時刻tでの状態ベクトルxtの尤度関数L(xt|yt)を導出する。
次に、ステップS409において、データ同化部115は、時刻tでの状態ベクトルxtの第1の確率密度関数p(xt|y0:t-Δt)と、当該時刻tでの状態ベクトルxtの尤度関数L(xt|yt)とを、ベイズ統計のモデリングによるデータ同化を行うフィルタのアルゴリズムに与えることにより、時刻tでの状態ベクトルxtの第2の確率密度関数p(xt|yt0:t)を導出する。
Next, in step S408, the likelihood
Next, in step S409, the
次に、ステップS410において、データ同化部115は、時刻tでの状態ベクトルxtの第2の確率密度関数p(xt|yt)の最頻値を、時刻tでの状態ベクトルの推定値として導出する。
次に、ステップS411において、可視化データ作成部116は、ステップS410で導出された時刻tでの状態ベクトルの推定値を構成する各位置における溶鋼3の流速と温度の表示データを作成する。
Next, in step S410, the
Next, in step S411, the visualization data creation unit 116 creates display data of the flow velocity and temperature of the
次に、ステップS412において、出力部117は、時刻tが、予め設定された時刻teになったか否かを判定する。この判定の結果、時刻tが、予め設定された時刻teになっていない場合、処理は、ステップS413に進む。
ステップS413において、第1の係数導出部111aは、(10)式~(12)式の微分方程式の解(au
1(t+Δt)~au
Mu(t+Δt)、av
1(t+Δt)~av
Mv(t+Δt)、aT
1(t+Δt)~aT
MT(t+Δt))を導出する際に使用するデータである時刻tでの溶鋼3の流速および温度と、時刻tから時刻t+Δtまでの間での境界条件のパラメータとの組を、ステップS409でデータ同化部115により導出された、時刻tでの状態ベクトルxtの第2の確率密度関数p(xt|yt0:t)に基づいて複数導出する。
Next, in step S412, the
In step S413, the first
そして、処理は、ステップS403に進み、時刻tが、予め設定された時刻teになるまで、ステップS403~S413の処理が繰り返し実行される。尚、ステップS404においては、ステップS404の直前のステップS413で導出された溶鋼3の流速および温度と、境界条件のパラメータとの組が用いられる。
以上のようにして、ステップS412において、時刻tが、予め設定された時刻teになると、処理は、ステップS414に進む。ステップS414において、出力部117は、可視化データ(ステップS411で作成された表示データと、ステップS410で導出された各時刻での状態ベクトルの推定値を構成する各位置における溶鋼3の流速と温度のデータ)を出力する。尚、表示データは、ステップS411で表示データが作成される度にコンピュータディスプレイに表示されるようにしてもよい。ステップS414の処理が終了すると、図4のフローチャートによる処理は終了する。
Then, the process proceeds to step S403, and the processes of steps S403 to S413 are repeatedly executed until the time t reaches a preset time t e . In step S404, the set of the flow velocity and temperature of the
In this manner, when the time t reaches a preset time t e in step S412, the process proceeds to step S414. In step S414, the
(実施例1)
次に、実施例1を説明するが、本発明は、以下の実施例1に限定されるものではない。
((発明例1))
まず、発明例1について説明する。
発明例1では、図18に示す構成の鉄鋼の連続鋳造機について、情報処理装置100によって可視化を行った。鋳型4の長辺の長さ(内法)は1500mm、短辺の長さ(内法)は250mmである。
Example 1
Next, a first embodiment will be described, but the present invention is not limited to the following first embodiment.
((Example 1))
First, Example 1 of the invention will be described.
In Example 1, visualization was performed by the
発明例1では、数値シミュレーション部101は、凝固シェル8の生成を考慮せずに(即ち、凝固シェル8は生成されないものとして)、有限体積法により、可視化対象面の各時刻tでの各計算格子点lにおける溶鋼3の流速および温度を計算した。可視化対象面の幅(l1軸方向の長さ)は約1500mm、深さ(l2軸方向の長さ)は、約8000mmとした。この可視化対象面の領域内に2600点の計算格子点lを設定した。また、溶鋼3の熱伝導度、比熱、密度、および粘性係数などの必要な物性値を設定した。
In Example 1, the
更に、伝熱の境界条件は、溶鋼3の外周境界において熱流束を定数で与えるものとした。また、流動の境界条件は、溶鋼3の湯面において流速の鉛直成分(l2軸方向の成分)を与えるものとし、鋳型4の内壁面と、浸漬ノズル6の外壁面と、溶鋼3の下端とにおいて流速の全成分を与えるものとし、吐出口7においては圧力を与えるものとした。そして、湯面における流速の鉛直成分(l2軸方向の成分)の初期値として、0(ゼロ)を設定(スリップ条件)した。また、鋳型4の内壁面と、浸漬ノズル6の外壁面とにおける流速の全成分の初期値として、0(ゼロ)を設定した。また、湯面から8000mm下側(l2軸の正の方向)に離れた位置(溶鋼3の下端)における流速の初期値として、l1軸方向を0(ゼロ)、l2軸方向を鋳造速度(鋳片の引き抜き速度)と同じに設定した。また、吐出口7における圧力の初期値として、浸漬ノズル6の外壁面との境界では0(ゼロ)以外の値をもつが、それらの平均値が0(ゼロ)になるように設定した。
Furthermore, the boundary condition of heat transfer was set to a constant heat flux at the outer boundary of the
発明例1では、数値シミュレーション部101は、各計算格子点lにおける溶鋼3の流速u、vおよび温度Tの初期値、並びに、溶鋼3の外周境界における熱流束分布、吐出口7における圧力、および鋳造速度を入力する。そうすると、数値シミュレーション部101は、0.025秒の時間隔の各時刻tでの各計算格子点lにおける溶鋼3の流速u(l,t)、v(l,t)および温度T(l,t)を数値シミュレーションにより算出する。
In Example 1, the
行列導出部102は、数値シミュレーション部101により導出された、各時刻tでの各計算格子点lにおける溶鋼3の流速u(l,t)、v(l,t)および温度T(l,t)から、0.25秒の時間隔で400個の時系列データを抽出し、流速u、v、温度Tのスナップショット行列MU、MV、MTを導出する。流速u、v、温度Tのスナップショット行列MU、MV、MTは、それぞれ2600行400列の行列である。
The
基底ベクトル導出部103は、スナップショット行列MU、MV、MTに対して特異値分解を実行することにより、流速u、v、温度Tの基底ベクトルφu
j(l)、φv
j(l)、φT
j(l)を導出する。また、基底ベクトル導出部103は、本実施形態で説明したようにして、(4)式、(5)式、(6)式の計算に用いる流速u、v、温度Tの基底ベクトルφu
j(l)、φv
j(l)、φT
j(l)の数Mu、Mv、MTを導出する。
The basis
図5は、溶鋼3の流速u、v、温度Tのスナップショット行列MU、MV、MTを特異値分解することにより得られるモード番号jの特異値(図5(a))と、モード番号1~jの特異値の積算値の割合(図5(b))の一例を示す図である。モード番号1~jの特異値の積算値の割合は、流速u、v、温度Tのスナップショット行列MU、MV、MTを特異値分解することにより得られるモード番号1~jの特異値の積算値を、流速u、v、温度Tのスナップショット行列MU、MV、MTを特異値分解することにより得られる全ての特異値の積算値で割った値を百分率で表したものである。尚、流速u、v、温度Tのスナップショット行列MU、MV、MTの列数は400であるので、モード番号が1~400のそれぞれについて特異値が導出される。しかしながら、表記の都合上、図5(a)および図5(b)では、高次のモードの特異値および特異値の積算値の割合の図示を省略する。モード番号1~jの特異値の積算値の割合が99.999%以上になるモード番号jの最小値を、(4)式、(5)式、(6)式の計算に用いる流速u、v、温度Tの基底ベクトルφu
j(l)、φv
j(l)、φT
j(l)の数Mu、Mv、MTとする。
5A and 5B are diagrams showing an example of singular values of mode number j obtained by singular value decomposition of snapshot matrices MU, MV, and MT of the flow velocity u, v, and temperature T of the molten steel 3 ( FIG. 5A ), and a ratio of the integrated value of singular values of
図5(b)に示す結果から、発明例1では、(4)式、(5)式、(6)式の計算に用いる流速u、v、温度Tの基底ベクトルφu j(l)、φv j(l)、φT j(l)の数Mu、Mv、MTを、それぞれ、15個、16個、20個とした。 From the results shown in FIG. 5(b), in Example 1, the numbers M u , M v , and M T of basis vectors φ u j (l), φ v j (l), and φ T j (l) of flow velocities u , v, and temperature T used in the calculations of equations (4), (5), and (6) were set to 15, 16, and 20, respectively.
縮約シミュレータ導出部104は、流速u、v、温度Tの基底ベクトルφu
1(l)~φu
15(l)、φv
1(l)~φv
16(l)、φT
1(l)~φT
20(l)および基底ベクトルの数Mu(=15)、Mv(=16)、MT(=20)に基づいて、(10)式~(12)式の定数Ai
u、Buu
ij、Buv
ij、Cuu
ijk、Cuv
ijk、Ai
v、Bvu
ij、Bvv
ij、Cvu
ijk、Cvv
ijk、Ai
T、BTT
ij、BTu
ij、BTv
ij、CTu
ijk、CTv
ijkを導出する。
The reduced
流速u、v、温度Tの基底ベクトルφu
j(l)、φv
j(l)、φT
j(l)は、各計算格子点lにおいて導出される。図6は、発明例1における流速u、v、温度Tの基底ベクトルφu
2(l)、φv
2(l)、φT
2(l)の分布を示す図である。図6では、流速u、v、温度Tの基底ベクトルφu
2(l)、φv
2(l)、φT
2(l)の大きさを等値線で示す。図6では、モード番号jが2である場合の流速u、v、温度Tの基底ベクトルφu
2(l)、φv
2(l)、φT
2(l)の分布を示す。図6に示すように、流速u、vの基底ベクトルφu
2(l)、φv
2(l)の分布は大きく異なることが分かる。
以上でオフラインにおける処理が終了する。
The basis vectors φ u j (l), φ v j (l), and φ T j (l) of the flow velocities u, v, and temperature T are derived at each computational grid point l. FIG. 6 is a diagram showing the distribution of the basis vectors φ u 2 (l), φ v 2 (l), and φ T 2 (l) of the flow velocities u, v, and temperature T in the invention example 1. In FIG. 6, the magnitudes of the basis vectors φ u 2 (l), φ v 2 (l), and φ T 2 (l) of the flow velocities u, v, and temperature T are shown by isopleths. In FIG. 6, the distribution of the basis vectors φ u 2 (l), φ v 2 (l), and φ T 2 (l) of the flow velocities u, v , and temperature T when the mode number j is 2 is shown. As shown in FIG. 6, it can be seen that the distributions of the basis vectors φ u 2 (l) and φ v 2 (l) of the flow velocities u and v are significantly different.
This completes the offline processing.
その後、オンラインにおいて、第1の係数導出部111aは、各時刻tにおいて、時刻t-Δtでの溶鋼3の流速と、時刻t-Δtでの溶鋼3の温度と、時刻t-Δtから時刻tまでの間での境界条件のパラメータとの少なくとも1つとのうち、少なくとも1つが異なる100ケースのそれぞれについて、(10)式~(12)式の微分方程式の解(au
1(t)~au
15(t)、av
1(t)~av
16(t)、aT
1(t)~aT
20(t))を導出する。
Thereafter, online, the first
物理量導出部111bは、各時刻tにおいて第1の係数導出部111aにより導出された、流速u、v、温度Tに対する係数au
1(t)~au
15(t)、av
1(t)~av
16(t)、aT
1(t)~aT
20(t)と、基底ベクトル導出部103により導出された、流速u、v、温度Tの基底ベクトルφu
1(l)~φu
15、φv
1(l)~φv
16、φT
1(l)~φT
20(l)とを、(4)式、(5)式、(6)式にそれぞれ代入することにより、時刻tでの各計算格子点lにおける溶鋼3の流速u(l,t)、v(l,t)、温度T(l,t)を導出する。これにより、各時刻tにおいて、多数の計算結果(数値解析データ(各計算格子点lにおける溶鋼3の流速および温度)のアンサンブル{xt
(k)}k)が得られる。
The physical
第1の確率密度関数導出部112は、多数の計算結果(数値解析データ(各計算格子点lにおける溶鋼3の流速および温度)のアンサンブル{xt
(k)}k)から、各時刻tでの状態ベクトルxtの第1の確率密度関数p(xt|yt-Δt)を導出する。
The first probability density
発明例1では、鋳型4に埋め込まれた温度計F1~F12、L1~L12で測定された温度を利用して観測データを得た。
温度計F1~F12、L1~L12は、図2に示すようにして配置されている。従って、観測データ取得部113は、情報処理装置100の測定対象面の数値シミュレーションと対応させるために、鋳型4の奥行方向(l3軸方向)において相互に対向する2つの温度計F、Lの合計12組のそれぞれについて、当該2つの温度計F、Lで同時刻に測定された温度の算術平均値を導出する。そして、観測データ取得部113は、この算術平均値を、当該温度の導出元となる2つの温度計F、Lが存在しているl1-l3平面上の位置(座標(l1,l3))の温度とし、その時系列データを求める。観測データ取得部113は、このようにして求めた12個の温度の時系列データを、鋳型4の冷却水の温度と、鋳型4の熱伝導度とに基づいて外挿し、可視化対象面の各観測位置αでの溶鋼3の温度を観測データyt
(0)として導出する。
In the invention example 1, the observation data was obtained using the temperatures measured by the thermometers F1 to F12 and L1 to L12 embedded in the
The thermometers F1 to F12 and L1 to L12 are arranged as shown in Fig. 2. Therefore, in order to correspond to the numerical simulation of the measurement target surface of the
尤度関数導出部114は、観測ベクトルytの観測ノイズとして、測定直前の一定期間(例えば、6時間または24時間)に導出された、複数の観測データ{yt-nΔt
(0),・・・,yt
(0)}を統計処理する。本実施例では、尤度関数導出部114は、時刻tでの観測ノイズは、平均ベクトルが0(ゼロベクトル)のガウス分布に従うものとし、その分散共分散行列値を算出する。これにより、各時刻tでの状態ベクトルxtの尤度関数L(xt|yt)が得られる。
The likelihood
データ同化部115は、各時刻tでの状態ベクトルxtの第1の確率密度関数p(xt|y0:t-Δt)と、各時刻tでの状態ベクトルxtの尤度関数L(xt|yt)とを、アンサンブルカルマンフィルタのアルゴリズムに与えることにより、フィルタ演算を行う。これにより、時刻tでの状態ベクトルxtの第2の確率密度関数p(xt|yt0:t)が得られる。(1)式において、観測行列Htは、各観測位置αにおける溶鋼3の温度を状態ベクトルxtから導出するように設定し、観測ノイズwtは、平均ベクトルが0(ゼロベクトル)で、分散共分散行列の分散値は、各温度計F1~F12、L1~L12の位置で評価した観測データytの変動の過去24時間分の分散値、共分散値は0(ゼロ)とした。
The
データ同化部115は、時刻tでの状態ベクトルxtの第2の確率密度関数p(xt|yt0:t)の最頻値を、時刻tでの状態ベクトルxtの推定値として導出する。
また、第1の係数導出部111aは、データ同化部115により導出された、時刻tでの状態ベクトルxtの第2の確率密度関数p(xt|yt0:t)から、時刻tでの各計算格子点lにおける溶鋼3の流速u(l,t)、v(l,t)および温度T(l,t)と、時刻tから時刻t+Δtまでの間での境界条件のパラメータとの組をランダムに100組導出する。そして、第1の係数導出部111aは、時刻tでの各計算格子点lにおける溶鋼3の温度T(l,t)および流速u(l,t)、v(l,t)の100組のそれぞれと、流速u、vの基底ベクトルφu
1(l)~φu
15、φv
1(l)~φv
16および温度Tの基底ベクトルφT
1(l)~φT
20(l)とに基づいて、(14)式~(16)式により、係数au
1(t)~au
15(t)、av
1(t)~av
16(t)、aT
1(t)~aT
20(t)のアンサンブルメンバーを100組導出する。そして、第1の係数導出部111aは、係数au
1(t)~au
15(t)、av
1(t)~av
16(t)、aT
1(t)~aT
20(t)のアンサンブルメンバーと、時刻tから時刻t+Δtまでの間での境界条件のパラメータとに基づいて、(10)式~(12)式により、係数au
1(t+Δt)~au
15(t+Δt)、av
1(t+Δt)~av
16(t+Δt)、aT
1(t+Δt)~aT
20(t+Δt)のアンサンブルメンバーを100組導出する。新たな境界条件は、浸漬ノズル6の吐出口7の位置での圧力の分布のみを、時刻tでの状態ベクトルxtの第2の確率密度関数p(xt|yt0:t)に基づいて修正し、その他の境界条件は、変更初期のままで固定した。時間隔Δtは1秒とした。
The
In addition, the first
可視化データ作成部116は、データ同化部115により導出された、時刻tでの状態ベクトルの推定値を構成する各位置における溶鋼3の流速と温度に基づいて、可視化データを作成する。出力部117は、可視化データ作成部116により可視化データが作成される度に、可視化データをコンピュータディスプレイに表示する。
The visualization data creation unit 116 creates visualization data based on the flow velocity and temperature of the
((発明例2))
次に、発明例2について説明する。
発明例2では、発明例1に対し、可視化対象面の領域内に設定する計算格子点lの数を2600から3600に増やした。発明例2では、発明例1と同様にして流速u、v、温度Tの基底ベクトルφu
j(l)、φv
j(l)、φT
j(l)の数Mu、Mv、MTを導出した結果、それぞれ、14個、14個、18個となった。発明例2におけるその他の処理および構成は、発明例1と同じである。
((Example 2))
Next, Example 2 of the invention will be described.
In Example 2, the number of computational grid points l set within the region of the visualization target surface is increased from 2600 to 3600 compared to Example 1. In Example 2, the numbers M u , M v , and M T of basis vectors φ u j (l), φ v j (l), and φ T j (l) of flow velocities u , v , and temperature T are derived in the same manner as in Example 1, and the numbers are 14, 14, and 18, respectively. Other processing and configurations in Example 2 are the same as those in Example 1.
((比較例1))
次に、比較例1について説明する。
比較例1では、発明例1において、溶鋼3の流速u、vに代えて、溶鋼3の運動エネルギーE(=(u2+v2)÷2)を用いる。また、比較例1では、流速u(l,t)、v(l,t)、温度T(l,t)は、それぞれ(4)式、(5)式、(6)式に代えて、以下の(17)式、(18)式、(19)式で表される。
(Comparative Example 1)
Next, comparative example 1 will be described.
In Comparative Example 1, the kinetic energy E (=( u2 + v2 )÷2) of the
ここで、φE
j(l)は、運動エネルギーEの基底ベクトルである。
比較例1では、数値シミュレーションにより、各時刻tでの各計算格子点lにおける溶鋼3の流速u、vから、各時刻tでの各計算格子点lにおける溶鋼3の運動エネルギーEを導出する。温度Tについては発明例1と同様にして、数値シミュレーションにより、各時刻tでの各計算格子点lにおける溶鋼3の温度Tを導出する。
Here, φ E j (l) is a basis vector of the kinetic energy E.
In Comparative Example 1, the kinetic energy E of the
各時刻tでの各計算格子点lにおける溶鋼3の運動エネルギーEおよび温度Tから、0.25秒の時間隔で400個の時系列データを抽出し、運動エネルギーE、温度Tのスナップショット行列ME、MTを導出する。運動エネルギーE、温度Tのスナップショット行列ME、MTは、それぞれ2600行400列の行列である。
400 time series data are extracted at 0.25 second intervals from the kinetic energy E and temperature T of the
スナップショット行列ME、MTに対して特異値分解を実行することにより、運動エネルギーE、温度Tの基底ベクトルφE j(l)、φT j(l)を導出する。また、(17)式、(18)式、(19)式の計算に用いる運動エネルギーE、温度Tの基底ベクトルφE j(l)、φT j(l)の数ME、MTを導出する。 By performing singular value decomposition on the snapshot matrices ME and MT, basis vectors φ E j (l) and φ T j (l) of kinetic energy E and temperature T are derived. In addition, the numbers M E and M T of basis vectors φ E j (l) and φ T j (l) of kinetic energy E and temperature T used in the calculations of equations (17), (18), and ( 19 ) are derived.
図7は、運動エネルギーEのスナップショット行列MEを特異値分解することにより得られるモード番号jの特異値(図7(a))と、モード番号1~jの特異値の積算値の割合(図7(b))の一例を示す図である。モード番号1~jの特異値の積算値の割合は、図5で説明したのと同様のものであり、運動エネルギーEのスナップショット行列MEを特異値分解することにより得られるモード番号1~jの特異値の積算値を、運動エネルギーEのスナップショット行列MEを特異値分解することにより得られる全ての特異値の積算値で割った値を百分率で表したものである。尚、図7(a)および図7(b)でも図5(a)および図5(b)と同様に、高次のモードの特異値および特異値の積算値の割合の図示を省略する。発明例1と同様に、モード番号1~jの特異値の積算値の割合が99.999%以上になるモード番号jの最小値を、(17)式、(18)式、(19)式の計算に用いる運動エネルギーE、温度Tの基底ベクトルφE
j(l)、φT
j(l)の数ME、MTとする。
7A and 7B are diagrams showing an example of the singular values of mode number j obtained by singular value decomposition of the snapshot matrix ME of the kinetic energy E ( FIG. 7A ), and the ratio of the integrated value of the singular values of
図7(b)に示す結果から、比較例1では、(17)式、(18)式の計算に用いる運動エネルギーEの基底ベクトルφe j(l)の数MEを38個とした。比較例1における温度Tのスナップショット行列MTは発明例1と同じであるので、温度Tの基底ベクトルφT j(l)の数MTは発明例1と同様に20個になる。 7B, in Comparative Example 1, the number M E of basis vectors φ e j (l) of kinetic energy E used in the calculations of Equations (17) and (18) is set to 38. Since the snapshot matrix MT of temperature T in Comparative Example 1 is the same as that in Invention Example 1, the number M T of basis vectors φ T j (l) of temperature T is 20, similar to Invention Example 1.
(7)式および(8)式に、(17)式および(18)式を代入した式の両辺に対して、i次のモードの運動エネルギーEの基底ベクトルφE i(l)との内積をとり、時間微分の項が左辺に来るように整理すると、それぞれ(10)式および(11)式と同じ式が得られる。また、(9)式に(17)式~(19)式を代入した式の両辺に対して、i次のモードの温度Tの基底ベクトルφT i(l)との内積をとり、時間微分の項が左辺に来るように整理すると、(12)式と同じ式が得られる。 Substituting equations (17) and (18) into equations (7) and (8), taking the inner product of both sides with the basis vector φ E i (l) of the kinetic energy E of the i-th mode and rearranging the equation so that the time derivative term is on the left side, gives the same equations as equations (10) and (11), respectively. Also, substituting equations (17) to (19) into equation (9), taking the inner product of both sides with the basis vector φ T i (l) of the temperature T of the i-th mode and rearranging the equation so that the time derivative term is on the left side, gives the same equation as equation (12).
比較例1では、運動エネルギーE、温度Tの基底ベクトルφE 1(l)~φE 38(l)、φT 1(l)~φT 20(l)および基底ベクトルの数Mu(=ME=38)、Mv(=ME=38)、MT(=20)に基づいて、(10)式~(11)式の定数Ai u、Buu ij、Buv ij、Cuu ijk、Cuv ijk、Ai v、Bvu ij、Bvv ij、Cvu ijk、Cvv ijk、Ai T、BTT ij、BTE ij、CTE ijkを導出する。 In comparative example 1, the constants A i u , B uu ij, B uv ij, C uu ijk, C uv ijk, A i v , B vu ij , B vv ij , C vu ijk, C vv ijk , A i T , B TT ij, B TE ij , and C TE ijk in equations (10) and (11) are derived based on the basis vectors φ E 1 (l) to φ E 38 (l), φ T 1 (l) to φ T 20 (l) of kinetic energy E and temperature T , and the numbers of basis vectors M u (= M E = 38 ) , M v ( = M E = 38 ) , and M T (= 20) .
図8は、比較例1における運動エネルギーEの基底ベクトルφE
2(l)の分布を示す図である。図8では、運動エネルギーEの基底ベクトルφE
2(l)の大きさを等値線で示す。図8では、モード番号jが2である場合の運動エネルギーEの基底ベクトルφE
2(l)の分布を示す。尚、温度Tの基底ベクトルφT
2(l)の分布は、図6に示したものと同じになる。図8に示すように、運動エネルギーEの基底ベクトルφE
2(l)の分布は、図6に示した流速u、vの基底ベクトルφu
2(l)、φv
2(l)の分布が合わさることにより、流速u、vの基底ベクトルφu
2(l)、φv
2(l)の分布とは大きく異なり、流速u、vの基底ベクトルφu
2(l)、φv
2(l)の分布の特徴が反映されなくなることが分かる。
以上でオフラインにおける処理が終了する。
FIG. 8 is a diagram showing the distribution of the basis vector φ E 2 (l) of the kinetic energy E in Comparative Example 1. In FIG. 8, the magnitude of the basis vector φ E 2 (l) of the kinetic energy E is shown by isopleths. In FIG. 8, the distribution of the basis vector φ E 2 (l) of the kinetic energy E when the mode number j is 2 is shown. The distribution of the basis vector φ T 2 (l) of the temperature T is the same as that shown in FIG. 6. As shown in FIG. 8, the distribution of the basis vector φ E 2 (l) of the kinetic energy E is significantly different from the distribution of the basis vectors φ u 2 (l) and φ v 2 (l) of the flow velocities u and v shown in FIG. 6 due to the combination of the distributions of the basis vectors φ u 2 (l) and φ v 2 (l) of the flow velocities u and v, and it can be seen that the characteristics of the distribution of the basis vectors φ u 2 (l) and φ v 2 (l) of the flow velocities u and v are no longer reflected.
This completes the offline processing.
その後、オンラインにおいて、発明例1と同様に、各時刻tにおいて、100ケースのそれぞれについて、(10)~(12)式の微分方程式の解(au
1(t)~au
38(t)、av
1(t)~av
38(t)、aT
1(t)~aT
20(t))を導出する。そして、各時刻tにおいて、流速u、v、温度Tに対する係数au
1(t)~au
38(t)、av
1(t)~av
38(t)、aT
1(t)~aT
20(t)と、オフラインで導出した運動エネルギーE、温度Tの基底ベクトルφE
1(l)~φE
38、φT
1(l)~φT
20を、(17)式、(18)式、(19)式にそれぞれ代入することにより、時刻tでの各計算格子点lにおける溶鋼3の流速u(l,t)、v(l,t)、温度T(l,t)を導出する。これにより、各時刻tにおいて、多数の計算結果(数値解析データ(各計算格子点lにおける溶鋼3の流速および温度)のアンサンブル{xt
(k)}k)が得られる。以降のオンラインにおける処理は、発明例1と同様であるので、当該処理の詳細な説明を省略する。
Then, online, in the same manner as in Example 1, solutions (a u 1 (t) to a u 38 (t), a v 1 (t) to a v 38 (t), a T 1 (t) to a T 20 ( t )) of the differential equations ( 10 ) to ( 12 ) are derived for each of the 100 cases at each time t. Then, at each time t, the flow velocity u ( l ,t), v( l ,t) and temperature T(l,t) of the
((比較例2))
次に、比較例2について説明する。
比較例2では、特許文献1と同様にして溶鋼3の流速および温度を導出して可視化データを作成した。即ち、数値シミュレーション部101、行列導出部102、基底ベクトル導出部103、および縮約シミュレータ導出部104を有さず、第1の係数導出部111aおよび物理量導出部111bに代えて特許文献1に記載されている数値シミュレーション部を用いることにより、溶鋼3の流速および温度を導出して可視化データを作成した。計算格子点lの数は、発明例1および比較例1と同様に2600個とした。
(Comparative Example 2)
Next, comparative example 2 will be described.
In Comparative Example 2, the flow velocity and temperature of the
((比較例3))
次に、比較例3について説明する。
比較例3は、比較例2に対し、可視化対象面の領域内に設定する計算格子点lの数を2600から3600に増やしたものである。即ち、比較例3は、特許文献1と同様にして溶鋼3の温度および流速を導出して可視化データを作成する際の計算格子点lの数を比較例2よりも増やしたものである。
(Comparative Example 3)
Next, comparative example 3 will be described.
In Comparative Example 3, the number of computational grid points l set within the region of the visualization target surface is increased from 2600 to 3600 compared to Comparative Example 2. That is, in Comparative Example 3, the number of computational grid points l when deriving the temperature and flow velocity of
((発明例と比較例との比較))
比較例2、3では、計算時間は長くなるが計算精度は高い。そこで、比較例2を基準として、発明例1および比較例1の計算精度を比較する。
図9は、発明例1における溶鋼3の温度と時間との関係を示す図である。図10は、比較例1における溶鋼3の温度と時間との関係を示す図である。尚、溶鋼3の温度は、可視化データの作成に用いた溶鋼3の温度である。
((Comparison between the invention example and the comparative example))
In Comparative Examples 2 and 3, the calculation time is long but the calculation accuracy is high. Therefore, the calculation accuracy of Example 1 and Comparative Example 1 will be compared with Comparative Example 2 as the standard.
Fig. 9 is a diagram showing the relationship between the temperature of the
図9において、グラフ901は、発明例1における溶鋼3の温度と時間との関係を示す。グラフ902は、比較例2における溶鋼3の温度と時間との関係を示す。図10において、グラフ1001は、比較例1における溶鋼3の温度と時間との関係を示す。グラフ1002は、比較例2における溶鋼3の温度と時間との関係を示し、グラフ902と同じである。何れのグラフ901、902、1001、1002も同一の計算格子点lにおける溶鋼3の温度を示す。
In FIG. 9,
図9に示すように、発明例1における溶鋼3の温度と時間との関係を示すグラフ901では、比較例2における溶鋼3の温度と時間との関係を示すグラフ902と同様に変化し、略同じ結果が得られた。これに対し、図10に示すように、比較例1における溶鋼3の温度と時間との関係を示すグラフ1001は、比較例2における溶鋼3の温度と時間との関係を示すグラフ1002と大きく異なる。特に、比較例1における溶鋼3の温度と時間との関係を示すグラフ1001では、20秒から30秒の時間帯における高温の状態や80秒以降の低温の状態を検出することができていない。
As shown in FIG. 9,
以上のような発明例1および比較例1における計算精度の差は、使用する基底ベクトルφu j(l)、φv j(l)、φE j(l)の違いに起因すると考えられる。そこで、流速u、vの基底ベクトルφu j(l)、φv j(l)を用いて比較例2における流速および温度を再現した場合と、運動エネルギーEの基底ベクトルφE j(l)を用いて比較例2における流速および温度を再現した場合とで、どのくらいの精度に差が生じるのかを調査した。 The difference in calculation accuracy between Example 1 and Comparative Example 1 is believed to be due to the difference in the basis vectors φuj ( l), φvj (l), and φEj ( l ) used. Therefore, an investigation was conducted to determine the difference in accuracy between reproducing the flow velocity and temperature in Comparative Example 2 using the basis vectors φuj ( l ) , φvj (l) of the flow velocities u and v, and reproducing the flow velocity and temperature in Comparative Example 2 using the basis vector φEj(l) of the kinetic energy E.
まず、比較例2において導出した流速u(l,t)、v(l,t)および温度T(l,t)と、比較例1において導出した運動エネルギーE、温度Tの基底ベクトルφE j(l)、φT j(l)とを用いて、以下の(20)式、(21)式と、(16)式とにより、流速u、v、温度Tに対する係数au j(t)、av j(t)、aT j(t)を導出した。また、比較例2において導出した流速u(l,t)、v(l,t)および温度T(l,t)と、発明例1において導出した流速u、v、温度Tの基底ベクトルφu j(l)、φv j(l)、φT j(l)とを用いて、(14)式、(15)式、(16)式により、係数au j(t)、av j(t)、aT j(t)を導出した。 First, coefficients a u j (t), a v j (t), and a T j (t) for the flow velocities u, v, and temperature T were derived from the following equations (20), (21), and (16) using the flow velocities u(l,t), v(l,t) and temperature T( l ,t) derived in Comparative Example 2, and the basis vectors φ E j (l), φ T j ( l ) of the kinetic energy E and temperature T derived in Comparative Example 1. Furthermore, coefficients a u j (t), a v j (t), and a T j (t) were derived from the following equations (14), (15), and (16) using the flow velocities u(l,t), v ( l ,t) and temperature T (l, t ) derived in Comparative Example 2, and the basis vectors φ u j (l), φ v j (l), and φ T j (l) of the flow velocities u, v, and temperature T derived in Invention Example 1.
そして、(20)式、(21)式、(16)式より導出した係数au j(t)、av j(t)、aT j(t)を、それぞれ、(17)式、(18)式、(19)式に代入する。このようにして得られる(17)式、(18)式、(19)式により、比較例2における流速u(l,t)、v(l,t)、温度T(l,t)を再現することができる。このような流速u(l,t)、v(l,t)、温度T(l,t)の再現を、運動エネルギーEの基底ベクトルφE j(l)の数MEおよび温度Tの基底ベクトルφT j(l)の数MTを異ならせて実行した。以下の説明では、このようにして再現した流速u(l,t)、v(l,t)、温度T(l,t)を、必要に応じて、比較例1の結果から再現した流速u^(l,t)、v^(l,t)、温度T^(l,t)と称する。 Then, the coefficients a u j (t), a v j (t), and a T j (t) derived from equations (20), (21), and (16) are substituted into equations (17), (18), and (19), respectively. The flow velocities u(l, t), v(l, t), and temperature T(l, t) in Comparative Example 2 can be reproduced by equations (17), (18), and (19) thus obtained. The reproduction of such flow velocities u(l, t), v(l, t), and temperature T(l, t) was performed by varying the number M E of basis vectors φ E j (l) of kinetic energy E and the number M T of basis vectors φ T j (l) of temperature T. In the following explanation, the flow velocities u(l,t), v(l,t), and temperature T(l,t) reproduced in this manner will be referred to as the flow velocities u^(l,t), v^(l,t), and temperature T^(l,t) reproduced from the results of Comparative Example 1, as necessary.
また、(14)式、(15)式、(16)式により導出した係数au j(t)、av j(t)、aT j(t)を、(4)式、(5)式、(6)式に代入する。このようにして得られる(4)式、(5)式、(6)式により、比較例2における流速u(l,t)、v(l,t)、温度T(l,t)を再現することができる。このような流速u(l,t)、v(l,t)、温度T(l,t)の再現を、流速uの基底ベクトルφu j(l)の数Mu、流速vの基底ベクトルφv j(l)の数Mv、および温度Tの基底ベクトルφT j(l)の数MTを異ならせて実行した。以下の説明では、このようにして再現した流速u(l,t)、v(l,t)、温度T(l,t)を、必要に応じて、発明例1の結果から再現した流速u^(l,t)、v^(l,t)、温度T^(l,t)と称する。 Furthermore, the coefficients a u j (t), a v j (t), and a T j (t) derived from equations (14), (15), and (16) are substituted into equations (4), (5), and (6). The flow velocities u(l, t), v(l, t), and temperature T(l, t) in Comparative Example 2 can be reproduced by equations (4), (5), and (6) thus obtained. The reproduction of such flow velocities u(l, t), v(l, t), and temperature T(l, t) was performed by varying the number M u of basis vectors φ u j (l) of flow velocity u, the number M v of basis vectors φ v j (l) of flow velocity v, and the number M T of basis vectors φ T j (l) of temperature T. In the following explanation, the flow velocities u(l,t), v(l,t), and temperature T(l,t) reproduced in this manner will be referred to as the flow velocities u^(l,t), v^(l,t), and temperature T^(l,t) reproduced from the results of Invention Example 1, as necessary.
そして、比較例1の結果から再現した流速u^(l,t)、v^(l,t)、温度T^(l,t)の、比較例2における流速u(l,t)、v(l,t)、温度T(l,t)に対する二乗平均平方根誤差(RMSE(Root Mean Squared Error))を導出した。同様に、発明例1の結果から再現した流速u^(l,t)、v^(l,t)、温度T^(l,t)の、比較例2における流速u(l,t)、v(l,t)、温度T(l,t)に対する二乗平均平方根誤差を導出した。尚、二乗平均平方根誤差は、使用した基底ベクトルの数ごと(ME、Mu、Mv、MTごと)に導出した。また、各物理変数(流速u、v、温度T)の二乗平均平方根誤差は、比較例2における当該物理変数(流速u、v、温度T)の平均値からの変動の二乗平均平方根値で割って規格化した。 Then, the root mean squared error (RMSE) of the flow velocities u(l,t), v(l,t), and temperature T(l,t) reproduced from the results of Comparative Example 1 with respect to the flow velocities u(l,t), v(l,t), and temperature T(l,t) in Comparative Example 2 was derived. Similarly, the root mean squared error of the flow velocities u(l,t), v(l,t), and temperature T(l,t) reproduced from the results of Invention Example 1 with respect to the flow velocities u(l,t), v(l,t), and temperature T(l,t) in Comparative Example 2 was derived. The root mean squared error was derived for each number of basis vectors used (for each M E , M u , M v , and M T ). In addition, the root mean square error of each physical variable (flow velocities u, v, temperature T) was normalized by dividing it by the root mean square value of the fluctuation of the physical variable (flow velocities u, v, temperature T) from the average value in Comparative Example 2.
図11は、流速uの二乗平均平方根誤差と、基底ベクトルφu j(l)、φE j(l)の数Mu、MEとの関係の一例を示す図である。図11(a)は、発明例1の結果から再現した流速u^(l,t)の、比較例2における流速u(l,t)に対する二乗平均平方根誤差(RMSE)を示す。図11(b)は、比較例1の結果から再現した流速u^(l,t)の、比較例2における流速u(l,t)に対する二乗平均平方根誤差(RMSE)を示す。 Fig. 11 is a diagram showing an example of the relationship between the root mean square error of the flow velocity u and the numbers M u and M E of basis vectors φ u j (l) and φ E j (l). Fig. 11(a) shows the root mean square error (RMSE) of the flow velocity u^(l,t) reproduced from the results of invention example 1 with respect to the flow velocity u(l,t) in comparative example 2. Fig. 11(b) shows the root mean square error (RMSE) of the flow velocity u^(l,t) reproduced from the results of comparative example 1 with respect to the flow velocity u(l,t) in comparative example 2.
図12は、流速vの二乗平均平方根誤差と、基底ベクトルφv j(l)、φE j(l)の数Mv、MEとの関係の一例を示す図である。図12(a)は、発明例1の結果から再現した流速v^(l,t)の、比較例2における流速v(l,t)に対する二乗平均平方根誤差(RMSE)を示す。図12(b)は、比較例1の結果から再現した流速v^(l,t)の、比較例2における流速v(l,t)に対する二乗平均平方根誤差(RMSE)を示す。 Fig. 12 is a diagram showing an example of the relationship between the root mean square error of the flow velocity v and the numbers Mv and M E of basis vectors φvj (l) and φEj (l) . Fig. 12(a) shows the root mean square error (RMSE) of the flow velocity v^(l,t) reproduced from the results of invention example 1 with respect to the flow velocity v(l,t) in comparative example 2. Fig. 12(b) shows the root mean square error (RMSE) of the flow velocity v^(l,t) reproduced from the results of comparative example 1 with respect to the flow velocity v(l,t) in comparative example 2.
図13は、温度Tの二乗平均平方根誤差と、基底ベクトルφT j(l)の数MTとの関係の一例を示す図である。尚、温度Tの基底ベクトルφT j(l)は、発明例1および比較例1で同じである。従って、温度Tの二乗平均平方根誤差と、基底ベクトルφT j(l)の数MTとの関係は、発明例1の結果から再現した温度T^(l,t)を用いても、比較例1の結果から再現した温度T^(l,t)を用いても同じ関係になる。 13 is a diagram showing an example of the relationship between the root mean square error of temperature T and the number M T of basis vectors φ T j (l). The basis vectors φ T j (l) of temperature T are the same in invention example 1 and comparative example 1. Therefore, the relationship between the root mean square error of temperature T and the number M T of basis vectors φ T j (l) is the same whether the temperature T^(l, t) reproduced from the results of invention example 1 is used or the temperature T^(l, t) reproduced from the results of comparative example 1 is used.
図11(a)において、発明例1で用いた流速uの基底ベクトルφu j(l)の数Mu(=15)であるときに、発明例1の結果から再現した流速u^(l,t)の、比較例2における流速u(l,t)に対する二乗平均平方根誤差は0.038(3.8%)程度になる。これに対し、図11(b)において、比較例1で用いた運動エネルギーEの基底ベクトルφE j(l)の数ME(=38)であるときに、比較例1の結果から再現した流速u^(l,t)の、比較例2における流速u(l,t)に対する二乗平均平方根誤差は0.65(65%)程度になる。 11(a), when the number M u (=15) of basis vectors φ u j (l) of flow velocity u used in invention example 1, the root mean square error of the flow velocity u^(l,t) reproduced from the results of invention example 1 with respect to the flow velocity u(l,t) in comparative example 2 is about 0.038 (3.8%). In contrast to this, in FIG 11(b), when the number M E (=38) of basis vectors φ E j (l) of kinetic energy E used in comparative example 1, the root mean square error of the flow velocity u^(l,t) reproduced from the results of comparative example 1 with respect to the flow velocity u(l,t) in comparative example 2 is about 0.65 (65%).
また、図12(a)において、発明例1で用いた流速vの基底ベクトルφv j(l)の数Mv(=16)であるときに、発明例1の結果から再現した流速u^(l,t)の、比較例2における流速u(l,t)に対する二乗平均平方根誤差は0.011(1.1%)程度になる。これに対し、図12(b)において、比較例1で用いた運動エネルギーEの基底ベクトルφE j(l)の数ME(=38)であるときに、比較例1の結果から再現した流速v^(l,t)の、比較例2における流速v(l,t)に対する二乗平均平方根誤差は0.50(50%)程度になる。 12(a), when the number M v (=16) of basis vectors φ v j (l) of flow velocity v used in invention example 1, the root mean square error of flow velocity u^(l,t) reproduced from the results of invention example 1 with respect to flow velocity u(l,t) in comparative example 2 is about 0.011 (1.1%). In contrast to this, when the number M E (=38) of basis vectors φ E j (l) of kinetic energy E used in comparative example 1 in Fig. 12(b), the root mean square error of flow velocity v^(l,t) reproduced from the results of comparative example 1 with respect to flow velocity v(l,t) in comparative example 2 is about 0.50 (50%).
このように、運動エネルギーEの基底ベクトルφE
j(l)を用いる比較例1では、比較例2における流速u、vの再現誤差が、本来の流速u、vの変動幅に匹敵するほど大きいことが分かる。比較例1では、本来の溶鋼3の流動現象を正しく反映できていないために、データ同化の結果も良好にならないためであると考えられる。これに対し、流速u、vの基底ベクトルφu
j(l)、φv
j(l)を個別に用いる発明例1では、比較例2における流速u、vを良好に再現することができることが分かる。
As described above, it is seen that in Comparative Example 1, which uses the basis vector φ E j (l) of the kinetic energy E, the reproduction error of the flow velocities u and v in Comparative Example 2 is large enough to rival the fluctuation range of the original flow velocities u and v. This is thought to be because Comparative Example 1 does not correctly reflect the original flow phenomenon of the
尚、図13に示すように、発明例1および比較例1で用いた温度Tの基底ベクトルφT j(l)の数MT(=20)であるときに、発明例1および比較例1の結果から再現した温度T^(l,t)の、比較例2における温度T(l,t)に対する二乗平均平方根誤差は0.090(9.0%)程度になる。従って、発明例1では、比較例2における温度Tを良好に再現することができることが分かる。 13, when the number M T (=20) of basis vectors φ T j (l) of temperature T used in invention example 1 and comparative example 1 is, the root mean square error of temperature T^(l, t) reproduced from the results of invention example 1 and comparative example 1 with respect to temperature T(l, t) in comparative example 2 is about 0.090 (9.0%). Therefore, it can be seen that invention example 1 can reproduce temperature T in comparative example 2 well.
次に、各発明例と各比較例の計算時間を比較する。表1は、発明例1、発明例2、比較例1、比較例2、および比較例3の計算時間を示す。発明例1および発明例2の計算時間は、オンラインにおいて、係数au
1(t)~au
Mu(t)、av
1(t)~av
Mv(t)、aT
1(t)~aT
MT(t)を導出するのに要する時間(第1の係数導出部111aにおける処理時間)である。比較例1の計算時間は、係数aE
1(t)~aE
ME(t)、aT
1(t)~aT
MT(t)を導出するのに要する時間である。比較例2および比較例3の計算時間は、特許文献1においてオンラインで実行される数値シミュレーションに要する時間(特許文献1に記載されている数値シミュレーション部における処理時間)である。尚、表1に示す数値の単位は秒である。
Next, the calculation time of each invention example and each comparative example is compared. Table 1 shows the calculation time of invention example 1, invention example 2, comparative example 1, comparative example 2, and comparative example 3. The calculation time of invention example 1 and invention example 2 is the time required to derive coefficients a u 1 (t) to a u Mu (t), a v 1 (t) to a v Mv (t), and a T 1 (t) to a T MT (t) online (processing time in the first
表1において、想定時間は、シミュレーションで想定している時間である。実操業において想定時間の欄に示される時間分のシミュレーションを実行するために必要な時間が、発明例1、発明例2、比較例1、比較例2、および比較例3の各欄に記載された計算時間である。表1に示すように、発明例1、発明例2、および比較例1では、3600秒分のシミュレーションを実行するのに必要な時間(想定時間が3600秒であるときの計算時間)は、その半分の1800秒を大きく下回る。従って、データ同化を実行しても、3600秒以内に可視化データを作成することができる。よって、リアルタイムで可視化データを作成することができる。一方、比較例2では、3600秒分のシミュレーションを実行するのに必要な時間は、3059秒であり、データ同化を実行し、可視化データを作成するまでの時間が3600秒を超える可能性がある。また、比較例3では、3600秒分のシミュレーションを実行するのに必要な時間は、6154秒であり、データ同化を実行する前に、計算時間が3600秒を大きく上回る。計算時間が3600秒を超える場合は、観測データが取得される時刻から可視化データを作成する時刻までのタイムラグが時間の経過とともに蓄積され、リアルタイムで可視化データを作成することができない。 In Table 1, the estimated time is the time assumed in the simulation. The time required to execute the simulation for the time shown in the estimated time column in actual operation is the calculation time listed in each column of invention example 1, invention example 2, comparative example 1, comparative example 2, and comparative example 3. As shown in Table 1, in invention example 1, invention example 2, and comparative example 1, the time required to execute a 3600-second simulation (the calculation time when the estimated time is 3600 seconds) is much less than half that, 1800 seconds. Therefore, even if data assimilation is performed, visualization data can be created within 3600 seconds. Therefore, visualization data can be created in real time. On the other hand, in comparative example 2, the time required to execute a 3600-second simulation is 3059 seconds, and the time required to execute data assimilation and create visualization data may exceed 3600 seconds. In comparative example 3, the time required to execute a 3600-second simulation is 6154 seconds, and the calculation time before executing data assimilation is much more than 3600 seconds. If the calculation time exceeds 3,600 seconds, the time lag between the time the observation data is acquired and the time the visualization data is created will accumulate over time, making it impossible to create the visualization data in real time.
発明例1、発明例2、比較例1、比較例2、および比較例3の結果を纏めると、表2のようになる。 The results of Example 1, Example 2, Comparative Example 1, Comparative Example 2, and Comparative Example 3 are summarized in Table 2.
POD適用は、固有直交分解の適用の有無を示す(〇は、固有直交分解の適用があることを示し、×は、固有直交分解の適用がないことを示す)。POD基底の数は、オンライン処理で使用する基底ベクトルφu
j(l)、φv
j(l)、φT
1(l)、φE
1(l)の総数(Mu、Mv、MT、MEの総数)を示す。解像度は、計算格子点lの総数を示す。精度は、溶鋼3の流速u、v、温度Tの計算精度の良否を示す(〇が良を示し、×が否を示す)。計算速度は、溶鋼3の流速u、v、温度Tのオンラインにおける計算速度の良否を示す(〇は、リアルタイムでの処理が十分に可能であることを示し、△は、リアルタイムで処理をするには時間の余裕がないことを示し、×は、リアルタイムでの処理ができないことを示す)。
POD application indicates whether proper orthogonal decomposition is applied (◯ indicates that proper orthogonal decomposition is applied, and × indicates that proper orthogonal decomposition is not applied). The number of POD bases indicates the total number of basis vectors φ u j (l), φ v j (l), φ T 1 (l), and φ E 1 (l) used in online processing (total number of M u , M v , M T , and M E ). Resolution indicates the total number of computational grid points l. Accuracy indicates whether the computation accuracy of the flow velocity u, v, and temperature T of the
表2に示すように、発明例1および発明例2では、精度および計算速度が共に良好である。これに対し、比較例1~3では、精度および計算速度の一方が発明例1および発明例2に対して劣ることが分かる。 As shown in Table 2, in invention examples 1 and 2, both the accuracy and calculation speed are good. In contrast, in comparison examples 1 to 3, it can be seen that either the accuracy or the calculation speed is inferior to invention examples 1 and 2.
(まとめ)
以上のように本実施形態では、情報処理装置100は、各計算格子点lおよび各時刻tにおける溶鋼3の流速u、vの値を格納するスナップショット行列を、スナップショット行列MU、MVという具合に流速の軸成分ごとにオフラインで導出する。そして、情報処理装置100は、スナップショット行列MU、MVに対して特異値分解を実行することにより、流速u、vの基底ベクトルφu
1(l)~φu
Mu(l)、φv
1(l)~φv
Mv(l)を導出することを、スナップショット行列MU、MVのそれぞれについてオフラインで実行する。その後、情報処理装置100は、流速u、vの基底ベクトルφu
1(l)~φu
Mu(l)、φv
1(l)~φv
Mv(l)を線形結合する際に用いる係数au
1(t)~au
Mu(t)、av
1(t)~av
Mv(t)を、軸成分ごとにオンラインで導出する。そして、情報処理装置100は、流速u、vの基底ベクトルφu
1(l)~φu
Mu(l)、φv
1(l)~φv
Mv(l)を係数au
1(t)~au
Mu(t)、av
1(t)~av
Mv(t)を用いて線形結合することにより、溶鋼3の流速u、vを軸成分ごとにオンラインで導出する。従って、非特許文献1に記載されているように運動エネルギーを用いて固有直交分解を実行する場合に比べて溶鋼3の流速u、vの推定精度を向上させることができる。また、非特許文献2に記載されているように複雑な計算式を用いなくても溶鋼の流速u、vを導出することができる。よって、連続鋳造設備4の鋳型内の溶鋼3の流速u、vを導出するために要する時間が長くなることと、溶鋼3の流速u、vの推定精度が低下することとの双方を抑制することができる。
(summary)
As described above, in this embodiment, the
また、本実施形態では、情報処理装置100は、それぞれが流速u、vの基底ベクトルφu
1(l)~φu
Mu(l)、φv
1(l)~φv
Mv(l)に乗算される係数au
1(t)~au
Mu(t)、av
1(t)~av
Mv(t)の時間変化を表す微分方程式であって、軸成分ごとに構成される微分方程式をオンラインで連立して解くことにより、係数au
1(t)~au
Mu(t)、av
1(t)~av
Mv(t)を、軸成分ごとにオンラインで導出する。前記微分方程式のそれぞれは、溶鋼3の流動を記述する運動方程式であって、溶鋼3の流速u、vを、流速u、vの基底ベクトルφu
1(l)~φu
Mu(l)、φu
1(l)~φv
Mv(l)の線形結合で表した運動方程式の両辺に対し、流速u、vの基底ベクトルφu
1(l)~φu
Mu(l)、φu
1(l)~φv
Mv(l)との内積をとることにより構成される微分方程式((10)式、(11)式のiを1,・・・,Mu、1,・・・,Mvとした微分方程式)である。従って、流速u、vの基底ベクトルφu
j(l)、φv
j(l)がそれぞれ正規直交系をなすことを利用して、縮約シミュレータを容易に導出することができる。
Furthermore, in this embodiment, the
また、本実施形態では、情報処理装置100は、各時刻tでの各計算格子点lにおける溶鋼3の温度Tの値を格納するスナップショット行列MTをオフラインで導出する。そして、情報処理装置100は、スナップショット行列MTに対して特異値分解を実行することにより、温度Tの基底ベクトルφT
1(l)~φT
MT(l)を導出することをオフラインで実行する。その後、情報処理装置100は、温度Tの基底ベクトルφu
1(l)~φT
MT(l)を線形結合する際に用いる係数aT
1(t)~aT
MT(t)をオンラインで導出する。そして、情報処理装置100は、温度Tの基底ベクトルφT
1(l)~φT
MT(l)を係数aT
1(t)~aT
MT(t)を用いて線形結合することにより、溶鋼3の温度Tをオンラインで導出する。従って、溶鋼3の流速u、vと共に温度Tも導出することができる。よって、連続鋳造設備4の鋳型内の溶鋼3の流速u、vおよび温度Tを導出するために要する時間が長くなることと、溶鋼3の流速u、vおよび温度Tの推定精度が低下することとの双方を抑制することができる。
Moreover, in this embodiment, the
また、本実施形態では、情報処理装置100は、温度Tの基底ベクトルφT
1(l)~φT
MT(l)に乗算される係数aT
1(t)~aT
MT(t)の時間変化を表す微分方程式をオンラインで連立して解くことにより、係数aT
1(t)~aT
MT(t)をオンラインで導出する。前記微分方程式のそれぞれは、溶鋼3の熱伝導を記述する熱伝導方程式であって、溶鋼3の温度Tを、温度Tの基底ベクトルφT
1(l)~φT
MT(l)の線形結合で表した運動方程式の両辺に対し、温度Tの基底ベクトルφT
1(l)~φT
MTとの内積をとることにより構成される微分方程式((12)式においてiを1,・・・,MTとした微分方程式)である。従って、温度Tの基底ベクトルφT
j(l)が正規直交系をなすことを利用して、縮約シミュレータを容易に導出することができる。
Moreover, in this embodiment, the
また、本実施形態では、情報処理装置100は、溶鋼3の流速u(l,t)、v(l,t)、温度T(l,t)を数値シミュレーションによりステップ幅Δt1の時間隔で、オフラインで導出する。そして、情報処理装置100は、ステップ幅Δt1の時間隔の各時刻tで導出された、各計算格子点lにおける流速u(l,t)、v(l,t)、温度T(l,t)から、ステップ幅Δt1よりも長いステップ幅Δt2の時間隔の各時刻tでの流速u(l,t)、v(l,t)、温度T(l,t)を抽出し、抽出した流速u(l,t)、v(l,t)、温度T(l,t)が成分として格納されたスナップショット行列MU、MV、MTを導出する。従って、特異値分解の際の計算負荷を軽減することができる。
In this embodiment, the
また、本実施形態では、情報処理装置100は、スナップショット行列MU、MV、MTに対して特異値分解を実行することにより得られる特異値に基づいて、スナップショット行列MU、MV、MTから導出できるモードの最大次数よりも小さい次数をMu、Mv、Mnの値として決定する。従って、過度に高次の基底ベクトルを用いなくても、実用上要求される精度で、溶鋼3の流速u、v、温度Tを、基底ベクトルφu
j(l)、φv
j(l)、φT
j(l)の線形結合で近似することができる。
Furthermore, in this embodiment, the
また、本実施形態では、情報処理装置100は、時刻tでの状態ベクトルxtの第1の確率密度関数p(xt|y0:t-Δt)と、当該時刻tでの状態ベクトルxtの尤度関数L(xt|yt)とを、データ同化を行うフィルタのアルゴリズムに与えることにより、時刻tでの状態ベクトルxtの第2の確率密度関数p(xt|y0:t)を導出する。情報処理装置100は、この第2の確率密度関数p(xt|y0:t)の最頻値を、時刻tでの状態ベクトルxtの推定値として導出する。また、情報処理装置100は、この第2の確率密度関数p(xt|y0:t)に基づいて、時刻t+Δtでの数値解析データxtを求めるための溶鋼3の流速および温度、並びに境界条件のパラメータを複数導出する。ここで、状態ベクトルxtの成分として含まれる溶鋼3の流速u、v、温度Tは、溶鋼3の流速u、v、温度Tの基底ベクトルφu
j(l)、φv
j(l)、φT
j(l)を、係数au
j(t)、av
j(t)、aT
j(t)を用いて線形結合することにより導出される。従って、連続鋳造機の内部における溶融金属の流動および温度の分布を、観測データによる事実と、物理法則に基づいた数値シミュレーションによる合理的な数値データとして可視化することを、大きな計算負荷をかけることなく実現することができる。よって、溶鋼3の偏流の様な異常現象の発生の様子を、時々刻々と知ることができる。これにより、鋳造の操業条件を臨機応変に制御したり、問題のある部材の変更、交換をタイミングよく行ったりすることが可能になり、高品質の鋼材を歩留まり良くかつ生産性良く製造することができる。
In this embodiment, the
(変形例)
本実施形態では、鋳型4の幅方向(l1軸方向)と鋳造方向(l2軸方向)とにより定まる2次元平面(l1-l2平面)における溶鋼3の流動および温度を可視化する場合を例示した。しかしながら、3次元空間における溶鋼3の流動および温度を可視化してもよい。このようにする場合、溶鋼3の流速のl1軸成分およびl2軸成分(流速u、v)に加えて、溶鋼3の流速のl3軸成分も溶鋼3の軸成分に含める。溶鋼3の流速のl3軸成分に対する処理は、本実施形態の溶鋼3の流速のl1軸成分およびl2軸成分(流速u、v)の説明において、l1軸成分およびl2軸成分をl3軸成分に置き換えることにより実現することができる。
(Modification)
In the present embodiment, the case where the flow and temperature of the
また、本実施形態において、情報処理装置100は、第2の確率密度関数p(xt|y0:t)の最頻値を、時刻tでの状態ベクトルxtの推定値として導出したが、第2の確率密度関数p(xt|y0:t)の平均値を、時刻tでの状態ベクトルxtの推定値として導出してもよいし、第2の確率密度関数p(xt|y0:t)の中央値を、時刻tでの状態ベクトルxtの推定値として導出してもよい。
Furthermore, in this embodiment, the
また、本実施形態において、可視化対象領域は、溶鋼3が存在する領域に限定したが、溶鋼3と凝固シェル8とが存在する領域にしてもよい。この場合、溶鋼3の物性値と、凝固シェル8の物性値と、凝固潜熱とを用いて、凝固を伴う非定常の流動および伝熱の数値シミュレーションを行えばよい。
In addition, in this embodiment, the visualization target area is limited to the area where the
また、本実施形態のようにデータ同化を実行し、溶鋼3の流速および温度を可視化すれば、溶鋼3の偏流の様な異常現象の発生の様子を、時々刻々と知ることができるので好ましい。しかしながら、必ずしもこのようにする必要はない。このようにする場合、例えば、尤度関数導出部114、データ同化部115、および可視化データ作成部116は不要になる。この場合、(10)式~(12)式の境界条件は、例えば、観測データ取得部113で取得された観測データに基づいて定めればよい。また、出力部117は、物理量導出部111bにより導出された各時刻tでの各計算格子点lにおける溶鋼3の流速u、v、温度Tを示す情報を出力する。また、溶鋼3の温度Tについては導出しなくてもよい。
In addition, by performing data assimilation as in this embodiment and visualizing the flow velocity and temperature of the
また、本実施形態では、数値シミュレーション部101、行列導出部102、基底ベクトル導出部103、および縮約シミュレータ導出部104を、オフラインで実行する場合を例示した。しかしながら、オフラインで実行される処理は、オンラインで実行される処理(第1の係数導出部111aにおける処理)の開始時よりも前の操業時にオンラインで実行してもよい。また、オンラインで実行される処理は、オフラインで実行してもよい。例えば、データ同化(第1の確率密度関数導出部112、観測データ取得部113、尤度関数導出部114、データ同化部115による処理)を実行せずに、数値シミュレーション部101、行列導出部102、基底ベクトル導出部103、縮約シミュレータ導出部104、第1の係数導出部111a、および物理量導出部111bにより、溶鋼3の流速uおよび温度Tを導出してもよい。このようにする場合、第1の係数導出部111aおよび物理量導出部111bは、オフラインで実行してもよい。
In addition, in this embodiment, the
〔第2の実施形態〕
次に、第2の実施形態を説明する。第1の実施形態では、各計算格子点lにおける溶鋼3の流速u(l,t)、v(l,t)、温度T(l,t)を状態量とし、各観測位置αにおける溶鋼3の温度Tと流速vを観測量とする。従って、第1の実施形態では、時刻t-Δtでの各計算格子点lにおける溶鋼3の流速u(l,t-Δt)、v(l,t-Δt)、温度T(l,t-Δt)と、流速u、v、温度Tの基底ベクトルφu
1(l)~φu
Mu(l)、φv
1(l)~φv
Mv(l)、φT
1(l)~φT
MT(l)とに基づいて、流速u、v、温度Tに対する係数au
1(t-Δt)~au
Mu(t-Δt)、av
1(t-Δt)~av
Mv(t-Δt)、aT
1(t-Δt)~aT
MT(t-Δt)を導出すること(ステップS404)と、流速u、v、温度Tに対する係数au
1(t)~au
Mu(t)、av
1(t)~av
Mv(t)、aT
1(t)~aT
MT(t)と、流速u、v、温度Tの基底ベクトルφu
1(l)~φu
Mu(l)、φv
1(l)~φv
Mv(l)、φT
1(l)~φT
MT(l)とに基づいて、時刻tでの各計算格子点lにおける溶鋼3の流速u(l,t)、v(l,t)、温度T(l,t)を導出すること(ステップS405)を、全てのアンサンブルメンバーに対してデータ同化を行う度に繰り返し実行する必要がある。このステップS405の計算負荷は、計算格子点lの数が多くなるほど高くなる。また、アンサンブルメンバーの数は、多く(例えば、100組)設定するほど、溶鋼3の流速u、vおよび温度Tの推定精度が向上することが期待される。このため、計算格子点lの数が多い場合には、ステップS405の処理を繰り返し実行すると、リアルタイムでの処理ができない虞がある。
Second Embodiment
Next, a second embodiment will be described. In the first embodiment, the flow velocity u(l, t), v(l, t), and temperature T(l, t) of the
そこで、本実施形態では、流速u、v、温度Tに対する係数au
1(t)~au
Mu(t)、av
1(t)~av
Mv(t)、aT
1(t)~aT
MT(t)を状態量として状態ベクトルxtを構成し、各観測位置αにおける溶鋼3の温度Tと流速vとから導出される温度T、流速vに対する後述の係数a^v
1(t)~a^v
Mv(t)、a^T
1(t)~a^T
MT(t)を観測量として観測ベクトルytを構成してデータ同化を実行することにより状態ベクトルの推定値を導出し、状態ベクトルの推定値を構成する係数を逆変換して、溶鋼3の流速と温度を導出する。従って、データ同化を行う度に、流速u、v、温度Tを係数au
1(t-Δt)~au
Mu(t-Δt)、av
1(t-Δt)~av
Mv(t-Δt)、aT
1(t-Δt)~aT
MT(t-Δt)に変換する上述の計算と、係数au
1(t)~au
Mu(t)、av
1(t)~av
Mv(t)、aT
1(t)~aT
MT(t)を流速u、v、温度Tに逆変換する上述の計算と、を全てのアンサンブルメンバーに対して実行する代わりに、流速v、温度Tを係数a^v
1(t)~a^v
Mv(t)、a^T
1(t)~a^T
MT(t)に変換する計算(後述するステップS1706)を流速v、温度Tの測定値のみに対して行えばよく、係数au
1(t)~au
Mu(t)、av
1(t)~av
Mv(t)、aT
1(t)~aT
MT(t)を流速u、v、温度Tに逆変換する計算(後述するステップS1712)を状態ベクトルの推定値のみに対して行えばよいので、計算負荷を軽減することができる。
Therefore, in this embodiment, a state vector xt is constructed using the coefficients a u 1 (t) to a u Mu (t), a v 1 (t) to a v Mv (t), and a T 1 (t) to a T MT (t) for the flow velocity u, v, and temperature T as state quantities, and an observation vector yt is constructed using the coefficients a ^ v 1 (t) to a ^ v Mv (t), a ^ T 1 (t) to a ^ T MT (t) for the temperature T and flow velocity v derived from the temperature T and flow velocity v of the
以上のように第1の実施形態では、第1の物理量を状態量とし、第2の物理量を観測量とするのに対し、本実施形態では、第1の物理量を基底ベクトルの線形結合で表す場合に基底ベクトルに乗算される係数を状態量とし、第2の物理量を基底ベクトルの線形結合で表す場合に基底ベクトルに乗算される係数を観測量とする。本実施形態と第1の実施形態とは、このことによる構成および処理が主として異なる。従って、本実施形態の説明において、第1の実施形態と同一の部分については図1~図13、図18~図19に付した符号と同一の符号を付す等して詳細な説明を省略する。 As described above, in the first embodiment, the first physical quantity is the state quantity and the second physical quantity is the observation quantity, whereas in this embodiment, when the first physical quantity is expressed as a linear combination of the basis vectors, the coefficient by which the basis vector is multiplied is the state quantity, and when the second physical quantity is expressed as a linear combination of the basis vectors, the coefficient by which the basis vector is multiplied is the observation quantity. This is the main difference between this embodiment and the first embodiment in terms of configuration and processing. Therefore, in the description of this embodiment, the same parts as those in the first embodiment are denoted by the same reference numerals as those in Figures 1 to 13 and Figures 18 to 19, and detailed description will be omitted.
図14は、情報処理装置1400の機能的な構成の一例を示す図である。情報処理装置1400のハードウェアは、例えば、図1に示した第1の実施形態の情報処理装置100と同じハードウェアで実現される。
Figure 14 is a diagram showing an example of the functional configuration of an
図1に示した第1の実施形態の情報処理装置100では、数値シミュレーション部101、行列導出部102、基底ベクトル導出部103、および縮約シミュレータ導出部104が、オフラインで実行される。これに対し本実施形態の情報処理装置1400は、数値シミュレーション部101、行列導出部102、基底ベクトル導出部103、および縮約シミュレータ導出部104に加えて、第2の係数導出部1405および観測行列導出部1406を有し、数値シミュレーション部101、行列導出部102、基底ベクトル導出部103、縮約シミュレータ導出部104、第2の係数導出部1405、および観測行列導出部1406が、オフラインで実行される。
In the
また、第1の実施形態の情報処理装置100では、状態導出部111、第1の確率密度関数導出部112、観測データ取得部113、尤度関数導出部114、データ同化部115、可視化データ作成部116、および出力部117が、オンラインで実行される。これに対し本実施形態の情報処理装置1400は、状態導出部1411、第1の確率密度関数導出部1412、観測データ取得部1413、尤度関数導出部1414、データ同化部1415、可視化データ作成部1416、および出力部1417に加えて、推定物理量導出部1418を有し、状態導出部1411、第1の確率密度関数導出部1412、観測データ取得部1413、尤度関数導出部1414、データ同化部1415、推定物理量導出部1418、可視化データ作成部1416、および出力部1417が、オンラインで実行される。また、状態導出部1411、第1の確率密度関数導出部1412、観測データ取得部1413、尤度関数導出部1414、データ同化部1415、可視化データ作成部1416、および出力部1417の機能の一部が、第1の実施形態の状態導出部111、第1の確率密度関数導出部112、観測データ取得部113、尤度関数導出部114、データ同化部115、可視化データ作成部116、および出力部117と異なる。
In addition, in the
(第2の係数導出部1405)
第1の実施形態で説明したように、(4)式に示す係数au
j(t)は、数値シミュレーション部101によりステップ幅Δt1の時間隔の各時刻tで導出された、各計算格子点lにおける流速u(l,t)、v(l,t)、温度T(l,t)に対して、流速uのスナップショット行列MUを特異値分解することにより得られるモード番号jの特異値とモード番号jの特異値に対する右特異ベクトルとの積として導出される。(5)式に示す係数av
j(t)は、数値シミュレーション部101によりステップ幅Δt1の時間隔の各時刻tで導出された、各計算格子点lにおける流速u(l,t)、v(l,t)、温度T(l,t)に対して、流速vのスナップショット行列MVを特異値分解することにより得られるモード番号jの特異値とモード番号jの特異値に対する右特異ベクトルとの積として導出される。(6)式に示す係数aT
j(t)は、数値シミュレーション部101によりステップ幅Δt1の時間隔の各時刻tで導出された、各計算格子点lにおける流速u(l,t)、v(l,t)、温度T(l,t)に対して、温度Tのスナップショット行列MTを特異値分解することにより得られるモード番号jの特異値とモード番号jの特異値に対する右特異ベクトルとの積として導出される。
(Second Coefficient Derivation Unit 1405)
As described in the first embodiment, the coefficient a u j (t) shown in the formula (4) is derived as a product of a singular value of mode number j obtained by singular value decomposition of a snapshot matrix MU of the flow velocity u for the flow velocity u(l,t), v(l,t), and temperature T(l,t) at each computational grid point l derived at each time t with a time interval of the step width Δt1 by the
第2の係数導出部1405は、行列導出部102により導出された流速uのスナップショット行列MUを特異値分解することにより得られるモード番号jの特異値とモード番号jの特異値に対する右特異ベクトルとの積を、係数au
j(t)として導出する。
また、第2の係数導出部1405は、行列導出部102により導出された流速vのスナップショット行列MVを特異値分解することにより得られるモード番号jの特異値とモード番号jの特異値に対する右特異ベクトルとの積を、係数av
j(t)として導出する。
また、第2の係数導出部1405は、行列導出部102により導出された温度Tのスナップショット行列MTを特異値分解することにより得られるモード番号jの特異値とモード番号jの特異値に対する右特異ベクトルとの積を、係数aT
j(t)として導出する。
The second coefficient derivation unit 1405 derives, as a coefficient a u j (t), the product of the singular value of mode number j obtained by singular value decomposition of the snapshot matrix MU of the flow velocity u derived by the
In addition, the second coefficient derivation unit 1405 derives, as a coefficient a v j (t), the product of the singular value of mode number j obtained by singular value decomposition of the snapshot matrix MV of the flow velocity v derived by the
In addition, the second coefficient derivation unit 1405 derives, as a coefficient a T j (t), the product of the singular value of mode number j obtained by singular value decomposition of the snapshot matrix MT of temperature T derived by the
(観測行列導出部1406)
観測行列導出部1406は、観測行列Htを導出する。(13)式に示したように、観測行列Htは、観測ベクトルyt(観測量)と、状態ベクトルxt(状態量)との関係を示す行列である。
(Observation matrix derivation unit 1406)
The observation
第1の実施形態では、計算格子点lにおける溶鋼3の流速u、vおよび温度Tの全てで構成されるベクトルを必要に応じて状態ベクトルと称し、物理量導出部111bにより導出される溶鋼3の流速u、vおよび温度Tを必要に応じて数値解析データと称することとした。これに対し、本実施形態では、溶鋼3の流速u、vおよび温度Tを基底ベクトルの線形結合で表す場合の当該基底ベクトルに乗算される係数au
1(t)~au
Mu(t)、av
1(t)~av
Mv(t)、aT
1(t)~aT
MT(t)を状態量とし、これらの係数の全てで構成されるベクトルを必要に応じて状態ベクトルと称することとする。また、後述する第1の係数導出部1411aにより導出される状態ベクトルを構成する係数au
1(t)~au
Mu(t)、av
1(t)~av
Mv(t)、aT
1(t)~aT
MT(t)を必要に応じて数値解析データと称することとする。
In the first embodiment, a vector constituted by all of the flow velocities u, v and temperature T of the
このように、溶鋼3の流速u、vおよび温度Tといった第1の物理量そのものではなく、第1の物理量を基底ベクトルの線形結合で表す場合に基底ベクトルに乗算される係数au
1(t)~au
Mu(t)、av
1(t)~av
Mv(t)、aT
1(t)~aT
MT(t)を用いて状態ベクトルxtを構成する場合、観測ベクトルytは、第2の物理量の測定値(温度計F、Lで測定される温度から導出される可視化対象面の各位置での溶鋼3の温度Tおよび湯面レベル計9で測定される湯面レベルから導出される流速v)ではなく、当該第2の物理量の測定値を、基底ベクトル導出部103で導出された流速v、温度Tの基底ベクトルφv
j(l)、φT
j(l)の線形結合で表す場合の当該基底ベクトルに乗算される係数a^v
j(t)、a^T
j(t)の全てで構成されるベクトルとする。
In this manner, when the state vector xt is constructed using the coefficients a u 1 (t) to a u Mu (t), a v 1 (t) to a v Mv (t), and a T 1 (t) to a T MT (t) by which the basis vector is multiplied when the first physical quantity is expressed as a linear combination of the basis vectors, rather than the first physical quantity itself, such as the flow velocity u, v , and temperature T of the molten steel 3, the observation vector yt is not the measured value of the second physical quantity (the temperature T of the
本発明者らは、以上のように係数で構成されるベクトルを状態ベクトルおよび観測ベクトルとする場合、係数a^v j(t)、a^T j(t)のそれぞれ異なる物理量に対する係数au j(t)、av j(t)、aT j(t)への依存度は、物理量そのもので構成されるベクトルを状態ベクトルおよび観測ベクトルとする場合よりも高くなるという知見を得た。このため、本発明者らは、観測ベクトルytの第α成分の係数と、状態ベクトルxtの第β成分の係数とが異なる物理量(流速u、vおよび温度T)に対する係数となっている場合でも、それらの係数の相関に応じて、観測行列Htの(α,β)成分を定めることが好ましいと考えた。尚、ここでの物理量が同一とは、軸成分も含めて同一であることを指す。従って、流速u、vは異なる物理量となる。 The inventors have found that when vectors composed of coefficients as described above are used as state vectors and observation vectors, the degree of dependence of the coefficients a^ vj (t) and a ^ Tj (t) on the coefficients auj (t ) , avj ( t ), and aTj ( t ) for different physical quantities is higher than when vectors composed of the physical quantities themselves are used as state vectors and observation vectors. For this reason, the inventors have found that even if the coefficient of the α-th component of the observation vector yt and the coefficient of the β-th component of the state vector xt are coefficients for different physical quantities (flow velocities u, v, and temperature T), it is preferable to determine the (α, β) components of the observation matrix Ht according to the correlation between these coefficients. Note that the physical quantities here are the same, including the axial components. Therefore, the flow velocities u and v are different physical quantities.
そこで、観測行列導出部1406は、流速uのスナップショット行列MU、流速vのスナップショット行列MV、温度Tのスナップショット行列MTを特異値分解することによって第2の係数導出部1405により導出される係数au
j(t)、av
j(t)、aT
j(t)についての相関を示す値を、観測行列Htの成分として導出する。観測行列Htは、以下の(22)式により表されるものとする。
Therefore, the observation
ここで、Huvは、観測ベクトルytの成分となる係数a^v j(t)と、状態ベクトルxtの成分となる係数au 1(t)~au Mu(t)との関係を示す成分を格納する観測行列Htの部分行列であり、Mv行Mu列の行列である。Hvvは、観測ベクトルytの成分となる係数a^v j(t)と、状態ベクトルxtの成分となる係数av 1(t)~av Mv(t)との関係を示す成分を格納する観測行列Htの部分行列であり、Mv行Mv列の行列である。HTvは、観測ベクトルytの成分となる係数a^v j(t)と、状態ベクトルxtの第β成分の係数aT 1(t)~aT MT(t)との関係を示す成分を格納する部分行列であり、Mv行MT列の行列である。HuTは、観測ベクトルytの成分となる係数a^T j(t)と、状態ベクトルxtの成分となる係数au 1(t)~au Mu(t)との関係を示す成分を格納する観測行列Htの部分行列であり、MT行Mu列の行列である。HvTは、観測ベクトルytの成分となる係数a^T j(t)と、状態ベクトルxtの成分となる係数au 1(t)~au Mu(t)との関係を示す成分を格納する観測行列Htの部分行列であり、MT行Mv列の行列である。HTTは、観測ベクトルytの成分となる係数a^T j(t)と、状態ベクトルxtの成分となる係数aT 1(t)~aT MT(t)との関係を示す成分を格納する観測行列Htの部分行列であり、MT行MT列の行列である。 Here, H uv is a submatrix of the observation matrix H t that stores components indicating the relationship between the coefficient a^ vj (t) that is a component of the observation vector y t and the coefficients a u 1 (t) to a u Mu ( t ) that are components of the state vector x t , and is a matrix with M v rows and M u columns. H vv is a submatrix of the observation matrix H t that stores components indicating the relationship between the coefficient a^ vj (t) that is a component of the observation vector y t and the coefficients a v 1 (t) to a vMv ( t) that are components of the state vector x t , and is a matrix with M v rows and M v columns. H Tv is a submatrix that stores components indicating the relationship between the coefficient a^ vj (t) that is a component of the observation vector y t and the coefficients a T 1 (t) to a T MT (t) of the β-th component of the state vector x t , and is a matrix with M v rows and M T columns. H uT is a submatrix of the observation matrix H t that stores components indicating the relationship between the coefficient a^ Tj ( t ) that is a component of the observation vector yt and the coefficients a u1 (t) to a uMu (t) that are components of the state vector xt , and is a matrix with M T rows and M u columns. H vT is a submatrix of the observation matrix H t that stores components indicating the relationship between the coefficient a^ Tj (t) that is a component of the observation vector yt and the coefficients a u1 (t) to a uMu (t) that are components of the state vector xt , and is a matrix with M T rows and M v columns. H TT is a submatrix of the observation matrix H t that stores components indicating the relationship between the coefficient a^ Tj (t) that is a component of the observation vector yt and the coefficients a T1 (t) to a TMT (t) that are components of the state vector xt , and is a matrix with M T rows and M T columns.
具体的に部分行列Huv、Hvv、HTv、HuT、HvT、HTTは、それぞれ、以下の(23)式~(28)式で表される。 Specifically, the submatrices H uv , H vv , H Tv , H uT , H vT , and H TT are respectively expressed by the following equations (23) to (28).
図15Aは、観測行列Htの部分行列HuTの具体例を示す図である。図15Bは、観測行列Htの部分行列HvTの具体例を示す図である。図15Cは、観測行列Htの部分行列HTTの具体例を示す図である。ここでは、表記スペースの関係で、温度Tについては固有値が大きいもの(固有直交分解のモード番号jが小さいもの)から順に8個、流速u、vについては固有値が大きいもの(固有直交分解のモード番号jが小さいもの)から順に6個の部分行列HuT、HvT、HTTのみを示す。 Fig. 15A is a diagram showing a specific example of a submatrix H uT of the observation matrix H t . Fig. 15B is a diagram showing a specific example of a submatrix H vT of the observation matrix H t . Fig. 15C is a diagram showing a specific example of a submatrix H TT of the observation matrix H t . Here, due to space limitations, only the submatrix H uT , H vT , and H TT are shown in order of the largest eigenvalue (smallest mode number j of proper orthogonal decomposition) for the temperature T, and the submatrix H uT , H vT , and H TT are shown in order of the largest eigenvalue (smallest mode number j of proper orthogonal decomposition) for the flow velocities u and v.
縮約シミュレータを構成する(10)式~(12)式では、物理変数を無次元化していない。従って、温度Tに関する係数aT j(t)と流速u、vに関する係数au j(t)、av j(t)との数値の大きさを補正するために、プラントル数Ptを用いると、(22)式の観測行列Htは、以下の(29)式のようになる。 In the equations (10) to (12) constituting the reduced simulator, the physical variables are not made dimensionless. Therefore, if the Prandtl number Pt is used to correct the magnitude of the values of the coefficient a T j (t) related to the temperature T and the coefficients a u j (t) and a v j (t) related to the flow velocities u and v, the observation matrix H t in equation (22) becomes the following equation (29).
本実施形態では、(29)式を観測行列Htとして用いる場合を例示する。尚、縮約シミュレータを構成する(10)式~(12)式において、物理変数を無次元化している場合には、(22)式を観測行列Htとして用いればよい。 In this embodiment, the case where equation (29) is used as the observation matrix Ht is exemplified. In addition, when the physical variables are non-dimensionalized in equations (10) to (12) constituting the reduced simulator, equation (22) may be used as the observation matrix Ht .
(状態導出部1411)
状態導出部1411は、溶鋼3の状態を示す量である状態量を成分として含む状態ベクトルxtを導出する。第1の実施形態の状態導出部111は、第1の係数導出部111aと、物理量導出部111bと、を有し、計算格子点lにおける溶鋼3の流速u、vおよび温度Tの全てで構成されるベクトルを状態ベクトルxtとして導出する。これに対し、本実施形態では、状態ベクトルxtは、係数au
1(t)~au
Mu(t)、av
1(t)~av
Mv(t)、aT
1(t)~aT
MT(t)の全てで構成されるベクトルであり、第1の係数導出部1411aにより導出される。従って、本実施形態の状態導出部1411は、第1の係数導出部1411aを有し、物理量導出部111bを有さない。
(State Derivation Unit 1411)
The
第1の係数導出部1411aは、第1の実施形態の第1の係数導出部111aと同じ処理を実行する。ただし、本実施形態では、第1の係数導出部1411aは、最初の時刻t(=t0)での係数au
1(t0)~au
Mu(t0)、av
1(t0)~av
Mv(t0)、aT
1(t0)~aT
MT(t0)と、時刻t0から時刻t0+Δtまでの間での境界条件のパラメータの少なくとも1つとのうち、少なくとも何れか1つが異なる複数のケースに対して、初期値として、予め設定された値を用いる。2番目の時刻t以降に関しては、後述するデータ同化部1415で導出された結果に基づいて、時刻tでの係数au
1(t0)~au
Mu(t0)、av
1(t0)~av
Mv(t0)、aT
1(t0)~aT
MT(t0)、並びに時刻tから時刻t+Δtまでの間での境界条件のパラメータを導出する。
尚、上述のように、本実施形態では、第1の物理量の一例である流速u、v、温度Tを係数au
1(t)~au
Mu(t)、av
1(t)~av
Mv(t)、aT
1(t)~aT
MT(t)に変換する計算は、実行しない。
The first
As described above, in this embodiment, calculations are not performed to convert the flow velocities u, v, and temperature T, which are examples of the first physical quantity, into coefficients a u 1 (t) to a u Mu (t), a v 1 (t) to a v Mv (t), and a T 1 (t) to a T MT (t).
(第1の確率密度関数導出部1412)
第1の確率密度関数導出部1412は、第1の実施形態の第1の確率密度関数導出部112と同様に、後述する観測データ取得部1413で取得された時刻t-Δtまでの観測データ(情報)yt0:t-Δtを基に推定される時刻tでの状態ベクトルxtの条件付き確率密度関数p(xt|yt0:t-Δt)を導出する。第1の確率密度関数導出部1412は、第1の係数導出部1411aにより導出された時刻tでの複数の数値解析データ{xt|t-Δt
(k)}kの相対度数(度数を全データ数で除した値)を用いた度数分布を導出する。この度数分布が、時刻tでの状態ベクトルxtの第1の確率密度関数p(xt|yt0:t-Δt)となる。第1の実施形態では、物理量導出部111bにより導出される溶鋼3の流速u、vおよび温度Tが、数値解析データとなる。これに対し、本実施形態では、第1の係数導出部1411aにより導出される係数au
1(t)~au
Mu(t)、av
1(t)~av
Mv(t)、aT
1(t)~aT
MT(t)が、数値解析データとなる。このように本実施形態の第1の確率密度関数導出部1412と第1の実施形態の第1の確率密度関数導出部112とは、数値解析データの内容が異なるだけであり、これらの処理内容は同じである。
(First probability density function derivation unit 1412)
The first probability density
(観測データ取得部1413)
観測データ取得部1413は、連続鋳造設備において測定された時系列データ(測定値)に基づいて観測データを所定の周期で取得する。第1の実施形態では、観測データyt
(0)は、第2の物理量の測定値(例えば、温度計F、Lで測定される温度から導出される可視化対象面の各位置での溶鋼3の温度Tおよび湯面レベル計9で測定される湯面レベルから導出される流速v)の全てで構成されるベクトルである。これに対し本実施形態では、観測データyt
(0)は、当該第2の物理量の測定値を基底ベクトル導出部103で導出された基底ベクトルの線形結合で表す場合の当該基底ベクトルに乗算される係数(例えば、流速v、温度Tを流速v、温度Tの基底ベクトルφv
j(l)、φT
j(l)の線形結合で表す場合の当該基底ベクトルに乗算される係数a^v
j(t)、a^T
j(t))の全てで構成されるベクトルである。
(Observation data acquisition unit 1413)
The observation data acquisition unit 1413 acquires observation data at a predetermined cycle based on time-series data (measured values) measured in the continuous casting equipment. In the first embodiment, the observation data y t (0) is a vector consisting of all of the measured values of the second physical quantity (for example, the temperature T of the
そこで、観測データ取得部1413は、第1の実施形態の観測データ取得部113と同様に、温度計F、Lで測定される温度(鋳型4内の温度)と、湯面レベル計9で測定される湯面レベルとを含む測定値を入力する。そして、観測データ取得部1413は、第1の実施形態の観測データ取得部113と同様に、鋳型4の長辺における伝熱状態が奥行方向(l3軸方向)の1次元定常熱伝導であるとの仮定をおいて、鋳型4の冷却水の温度と、鋳型4の熱伝導度とに基づいて、これら12個の温度を外挿し、可視化対象面の各位置での溶鋼3の温度を導出する。また、観測データ取得部1413は、第1の実施形態の観測データ取得部113と同様に、湯面レベル計9で測定された複数の湯面レベルの単位時間当たりの変化から、可視化対象面の溶鋼3の湯面に対応する各位置での鋳造方向(l2軸方向)の流速を導出する。本実施形態でも、第1の実施形態と同様に、以上のようにして導出される溶鋼3の温度の導出位置(可視化対象面の各位置)と、流速vの導出位置(可視化対象面の溶鋼3の湯面に対応する各位置)を観測位置とする。
Thus, the observation data acquisition unit 1413, like the observation data acquisition unit 113 in the first embodiment, inputs measured values including the temperatures (temperatures in the mold 4) measured by the thermometers F and L and the molten steel level measured by the molten steel level gauge 9. Then, like the observation data acquisition unit 113 in the first embodiment, the observation data acquisition unit 1413 assumes that the heat transfer state in the long side of the
そして、観測データ取得部1413は、可視化対象面Sの範囲内で、基底ベクトル導出部103によりオフラインで導出された温度Tの基底ベクトルφT
j(l)と、観測位置での温度T(l,t)との内積をとることを、モード番号jのそれぞれについて個別に実行することにより、観測位置での温度T(l,t)を、基底ベクトル導出部103により導出された温度Tの基底ベクトルφT
j(l)を用いた線形結合で表す場合に当該基底ベクトルφT
j(l)に乗算される係数a^T
j(t)を導出する。
Then, the observation data acquisition unit 1413 takes the inner product of the basis vector φ T j (l) of the temperature T derived offline by the basis
また、観測データ取得部1413は、融液自由表面Cの範囲内で、基底ベクトル導出部103によりオフラインで導出された流速vの基底ベクトルφv
j(l)と、観測位置での流速v(l,t)との内積をとることを、モード番号jのそれぞれについて個別に実行することにより、観測位置での流速v(l,t)を、基底ベクトル導出部103により導出された流速vの基底ベクトルφv
j(l)を用いた線形結合で表す場合に当該基底ベクトルφv
j(l)に乗算される係数a^v
j(t)を導出する。融液自由表面Cは、可視化対象面Sの境界領域のうち、溶鋼3の湯面に対応する境界領域である。
Furthermore, the observation data acquisition unit 1413 individually performs, for each mode number j, an inner product of the basis vector φ v j (l) of the flow velocity v derived offline by the basis
具体的に、観測データ取得部1413は、以下の(30)式、(31)式の計算を行うことにより、係数a^T j(t)、a^v j(t)を導出する。 Specifically, the observed data acquiring unit 1413 derives the coefficients a^ T j (t) and a^ v j (t) by performing calculations according to the following equations (30) and (31).
尚、係数a^T
j(t)、a^v
j(t)はそれぞれ、(30)式の左辺の記号、(31)式の左辺の記号(aの上に^が付され、aの右上にT、vが付され、T、vの下にjが付された記号の右隣に(t)が付された記号)を示す。
本実施形態では、以上のようにして観測データ取得部1413で導出される係数a^T
j(t)、a^v
j(t)の全てで構成されるベクトルが、時刻tでの観測データyt
(0)になる。(30)式、(31)式に示すように、係数a^T
j(t)、a^v
j(t)を導出する際の積分範囲は、可視化対象面S、融液自由表面Cであるので、観測データ取得部113に対する観測データ取得部1413の計算負荷の増加は過大にならない。
The coefficients a^ Tj (t) and a ^ vj (t) respectively indicate the symbols on the left side of equation (30) and equation (31) (the symbols with ^ above a, T and v to the upper right of a, and (t) to the right of the symbol with j below T and v).
In this embodiment, a vector composed of all of the coefficients a^ Tj (t) and a^ vj ( t) derived by the observation data acquisition unit 1413 in the above manner becomes the observation data yt (0) at time t. As shown in equations (30) and (31), the integration ranges when deriving the coefficients a ^Tj ( t) and a ^ vj (t) are the visualization target surface S and the melt free surface C, so that the increase in the calculation load of the observation data acquisition unit 1413 on the observation data acquisition unit 113 does not become excessive.
尚、本実施形態では、第2の物理量の一例である流速v、温度Tは、第1の物理量の一例である流速u、v、温度Tに含まれている必要がある。第1の実施形態と異なり、本実施形態では、観測データ取得部1413において、第2の物理量の基底ベクトルを用いる必要があるからである。これに対し、第1の実施形態では、第2の物理量を観測量としているので、第1の物理量が第2の物理量を含んでいなくても、(13)式に示す観測方程式を構成することができる。従って、第1の実施形態では、必ずしも第1の物理量が第2の物理量を含んでいる必要はない。 In this embodiment, the flow velocity v and temperature T, which are examples of the second physical quantity, must be included in the flow velocity u, v, and temperature T, which are examples of the first physical quantity. This is because, unlike the first embodiment, in this embodiment, the observation data acquisition unit 1413 must use a basis vector of the second physical quantity. In contrast to this, in the first embodiment, the second physical quantity is used as the observation quantity, so that the observation equation shown in equation (13) can be constructed even if the first physical quantity does not include the second physical quantity. Therefore, in the first embodiment, it is not necessarily necessary for the first physical quantity to include the second physical quantity.
(尤度関数導出部1414)
尤度関数導出部1414は、第1の実施形態の尤度関数導出部114と同様に、観測データ取得部1413で取得された時刻tでの観測データyt
(0)の観測ノイズを、過去の当該観測データから評価することにより、時刻tでの状態ベクトルxtが得られたときの観測ベクトルytの条件付き確率密度関数p(yt|xt)に対する状態ベクトルxtの尤度関数L(xt|yt)を導出する。第1の実施形態で説明したように、尤度関数L(xt|yt)とは、条件付き確率密度関数p(yt|xt)を状態ベクトルxtの関数と見做したもので、観測ベクトルytに観測データyt
(0)を代入することによって状態ベクトルxtの関数として確定する。
(Likelihood function derivation unit 1414)
The likelihood
第1の実施形態では、観測ベクトルyt
)は、観測位置において取得される溶鋼3の温度Tおよび鋳造方向の流速u、vの全てで構成されるベクトルであり、観測ベクトルytのデータである観測データyt
(0)は、観測データ取得部113により取得される。これに対し、本実施形態では、観測ベクトルytは、係数a^T
j(t)、a^v
j(t)の全てで構成されるベクトルであり、観測ベクトルytのデータである観測データyt
(0)は、観測データ取得部1413により取得される。このように本実施形態の尤度関数導出部1414と第1の実施形態の尤度関数導出部114とは、観測ベクトルytおよび観測データyt
(0)の内容が異なるだけであり、これらの処理内容は同じである。
In the first embodiment, the observation vector yt ) is a vector composed of all of the temperature T of the
(データ同化部1415)
データ同化部1415は、第1の実施形態のデータ同化部115と同様に、第1の確率密度関数導出部1412により導出された時刻tでの状態ベクトルxtの第1の確率密度関数p(xt|y0:t-Δt)と、尤度関数導出部1414により導出された時刻tでの観測データが得られたときの当該時刻tでの状態ベクトルxtの尤度関数L(xt|yt)とを、ベイズの定理を基礎としたベイズ統計のモデリングによるデータ同化を行うフィルタのアルゴリズムに与えることにより、時刻tまでの観測データが得られたときの時刻tでの状態ベクトルxtの条件付き確率密度関数p(xt|yt0:t)を第2の確率密度関数として導出する。本実施形態のデータ同化部1415と第1の実施形態のデータ同化部115とは、状態ベクトルxt、第1の確率密度関数p(xt|y0:t-Δt)および尤度関数L(xt|yt)の内容が異なるだけであり、これらの処理内容は同じである。
(Data Assimilation Unit 1415)
Similar to the
尚、データ同化部1415は、状態ベクトルxtの第2の確率密度関数p(xt|yt)を第1の係数導出部1411aに出力する。第1の係数導出部1411aは、状態ベクトルxtの第2の確率密度関数p(xt|yt)のアンサンブルメンバーに基づいて、時刻tでの係数au
1(t0)~au
Mu(t)、av
1(t)~av
Mv(t)、aT
1(t)~aT
MT(t)、並びに時刻tから時刻t+Δtまでの間での境界条件のパラメータを導出する。即ち、第1の係数導出部1411aは、第1の実施形態のように、時刻tでの各計算格子点lにおける溶鋼3の流速u(l,t)、v(l,t)および温度T(l,t)と、流速u、vの基底ベクトルφu
1(l)~φu
Mu(l)、φv
1(l)~φv
Mv(l)および温度Tの基底ベクトルφT
1(l)~φT
MT(l)と、に基づいて、(14)式~(16)式により、時刻tでの係数au
1(t)~au
Mu(t)、av
1(t)~av
Mv(t)、aT
1(t)~aT
MT(t)を導出する必要はなく、時刻tでの係数au
1(t0)~au
Mu(t0)、av
1(t0)~av
Mv(t0)、aT
1(t0)~aT
MT(t0)を、状態ベクトルxtの第2の確率密度関数p(xt|yt)のアンサンブルメンバーから直接的に導出する。
The
(推定物理量導出部1418)
推定物理量導出部1418は、データ同化部1415により状態ベクトルの推定値として導出された、時刻tでの状態ベクトルの推定値を構成する係数au
1(t)~au
Mu(t)と、基底ベクトル導出部103により導出された、流速uの基底ベクトルφu
1(l)~φu
Mu(l)とを(4)式に代入することにより、時刻t0~teの各時刻tでの各計算格子点lにおける溶鋼3の流速uを導出する。また、推定物理量導出部1418は、データ同化部1415により状態ベクトルの推定値として導出された、時刻tでの状態ベクトルの推定値を構成する係数av
1(t)~av
Mv(t)と、基底ベクトル導出部103により導出された、流速vの基底ベクトルφv
1(l)~φv
Mv(l)とを(5)式に代入することにより、時刻t0~teの各時刻tでの各計算格子点lにおける溶鋼3の流速vを導出する。また、推定物理量導出部1418は、データ同化部1415により状態ベクトルの推定値として導出された、時刻tでの状態ベクトルの推定値を構成する係数aT
1(t)~aT
MT(t)と、基底ベクトル導出部103により導出された温度Tの基底ベクトルφT
1(l)~φT
MT(l)とを(6)式に代入することにより、時刻t0~teの各時刻tでの各計算格子点lにおける溶鋼3の温度Tを導出する。時刻teは、予め設定された時刻であり、可視化データの作成期間の終期の時刻である。このように本実施形態では、推定物理量導出部1418は、データ同化部115によるデータ同化が終了した後に、各時刻tでの各計算格子点lにおける溶鋼3の流速u、v、温度Tを導出する。
(Estimated physical quantity deriving unit 1418)
The estimated physical quantity derivation unit 1418 derives the flow velocity u of
(可視化データ作成部1416)
可視化データ作成部1416は、推定物理量導出部1418により導出された時刻tでの各位置(各計算格子点l)における溶鋼3の流速u、v、温度Tの少なくとも何れか1つの表示データを作成する。第1の実施形態の可視化データ作成部116では、データ同化部115から、時刻tでの状態ベクトルの推定値を各位置における溶鋼3の温度Tと流速u、vとして取得する。これに対し、本実施形態の可視化データ作成部1416では、推定物理量導出部1418により導出された各時刻tでの各位置における溶鋼3の流速u、v、温度Tを取得する。第1の実施形態の可視化データ作成部116と本実施形態の可視化データ作成部1416とは、この点が異なるだけであり、これらの処理内容は同じである。
(Visualization Data Creation Unit 1416)
The visualization data creation unit 1416 creates display data of at least one of the flow velocities u, v, and temperature T of the
(出力部1417)
出力部1417は、第1の実施形態の出力部117と同様に、可視化データ作成部1416で作成された表示データをコンピュータディスプレイに表示する。また、出力部1417は、このような表示データに代えてまたは加えて、推定物理量導出部1418により導出された各時刻での各位置(各計算格子点l)における溶鋼3の流速と温度のデータを出力することができる。第1の実施形態の出力部117では、データ同化部115から、時刻tでの状態ベクトルの推定値を各位置における溶鋼3の温度Tと流速u、vとして取得する。これに対し、本実施形態の出力部117では、推定物理量導出部1418により導出された各時刻tでの各位置における溶鋼3の流速u、v、温度Tを取得する。第1の実施形態の出力部117と本実施形態の出力部1417とは、この点が異なるだけであり、これらの処理内容は同じである。
(Output unit 1417)
The
(動作フローチャート)
次に、図16のフローチャートを参照しながら、情報処理装置1400のオフラインにおける処理の一例を説明する。
図16のステップS301~S303、S304の処理は、図3のフローチャートのステップS301~S303、S304の処理と同じである。
(Operation Flowchart)
Next, an example of offline processing by the
The processes in steps S301 to S303 and S304 in FIG. 16 are the same as the processes in steps S301 to S303 and S304 in the flowchart in FIG.
即ち、ステップS301において、流速u(l,t)、v(l,t)および温度T(l,t)が導出される。次に、ステップS302において、流速uのスナップショット行列MU、流速vのスナップショット行列MV、温度Tのスナップショット行列MTが導出される。次に、ステップS303において、流速uの基底ベクトルφu j(l)、流速vの基底ベクトルφv j(l)、温度Tの基底ベクトルφT j(l)と、流速u、v、温度Tの基底ベクトルφu j(l)、φv j(l)、φT j(l)の数Mu、Mv、MTとが導出される。 That is, in step S301, flow velocities u(l,t), v(l,t) and temperature T(l,t) are derived. Next, in step S302, a snapshot matrix MU of flow velocity u, a snapshot matrix MV of flow velocity v, and a snapshot matrix MT of temperature T are derived. Next, in step S303, a basis vector φ u j (l) of flow velocity u, a basis vector φ v j (l) of flow velocity v, a basis vector φ T j (l) of temperature T, and numbers M u , M v , and M T of the basis vectors φ u j (l), φ v j (l), and φ T j (l) of flow velocities u , v , and temperature T are derived.
ステップS303の処理が終わると、処理は、ステップS1601に進む。ステップS1601において、第2の係数導出部1405は、ステップS302で導出された、流速u、v、温度Tのスナップショット行列MU、MV、MTを特異値分解することにより得られるモード番号jの特異値とモード番号jの特異値に対する右特異ベクトルとの積を、係数au j(t)、av j(t)、aT j(t)として導出する。尚、ステップS1601はステップS303の前に実行してもよい。 When the process of step S303 ends, the process proceeds to step S1601. In step S1601, the second coefficient derivation unit 1405 derives the product of the singular value of mode number j obtained by singular value decomposition of the snapshot matrices MU, MV, MT of the flow velocities u, v and temperature T derived in step S302 and the right singular vector for the singular value of mode number j as coefficients a u j (t), a v j (t), and a T j (t). Note that step S1601 may be executed before step S303.
ステップS1601の処理が終わると、処理は、ステップS304の処理に進み、縮約シミュレータが導出される。第1の実施形態で説明したように、縮約シミュレータは、(10)式~(12)式において、定数Ai
u、Buu
ij、Buv
ij、Cuu
ijk、Cuv
ijk、Ai
v、Bvu
ij、Bvv
ij、Cvu
ijk、Cvv
ijk、Ai
T、BTT
ij、BTu
ij、BTv
ij、CTu
ijk、CTv
ijk、Mu、Mv、MTが設定されたものである。ステップS304の処理が終了すると、処理は、ステップS1602に進む。
When the process of step S1601 is completed, the process proceeds to step S304, where a contraction simulator is derived. As described in the first embodiment, the contraction simulator is a simulator in which constants A i u , B uu ij , B uv ij , C uu ijk , C uv ijk , A i v , B vu ij , B vv ij , C vu ijk , C vv ijk , A i T , B TT ij , B Tu ij , B Tv ij , C Tu ijk , C Tv ijk , M u , M v , and M T are set in
ステップS1602において、観測行列導出部1406は、ステップS1601で導出された係数au
j(t)、av
j(t)、aT
j(t)についての相関を示す値を、観測行列Htの成分として導出し、観測行列Htを導出する((22)式~(28)式を参照)。ステップS1602の処理が終了すると、図16のフローチャートによる処理は終了する。尚、ステップS1602は、ステップS304の前に実行してもよい。この場合、ステップS304の処理が終了すると、図16のフローチャートによる処理は終了する。
In step S1602, the observation
次に、図17のフローチャートを参照しながら、情報処理装置1400のオンラインにおける処理の一例を説明する。
まず、ステップS1701において、第1の係数導出部1411aは、時刻tを初期値(t0)に設定する。
次に、ステップS1702において、第1の係数導出部1411aは、時刻tでの係数au
1(t)~au
Mu(t)、av
1(t)~av
Mv(t)、aT
1(t)~aT
MT(t)と、時刻tから時刻t+Δtまでの間での境界条件のパラメータとの組の初期値として、複数の組の初期値を設定する。
Next, an example of online processing by the
First, in step S1701, the first
Next, in step S1702, the first
次に、ステップS1703において、第1の係数導出部1411aは、時刻tを時刻t+Δtに更新して、処理は、ステップS1704に進む。
次に、ステップS1704において、第1の係数導出部1411aは、時刻tでの(10)式~(12)式の微分方程式の解(au
1(t)~au
Mu(t)、av
1(t)~av
Mv(t)、aT
1(t)~aT
MT(t))を導出する。係数au
1(t)~au
Mu(t)、av
1(t)~av
Mv(t)、aT
1(t)~aT
MT(t)は、時刻t-Δtでの各計算格子点lにおける溶鋼3の流速u(l,t-Δt)、v(l,t-Δt)および温度T(l,t-Δt)と、時刻t-Δtから時刻tまでの間での境界条件のパラメータとの組のそれぞれについて導出される。
Next, in step S1703, the first
Next, in step S1704, the first
次に、ステップS1705において、第1の確率密度関数導出部1412は、時刻tでの状態ベクトルxt(係数au
1(t)~au
Mu(t)、av
1(t)~av
Mv(t)、aT
1(t)~aT
MT(t))の第1の確率密度関数p(xt|yt-Δt)を導出する。
Next, in step S1705, the first probability density
次に、ステップS1706において、観測データ取得部1413は、温度計F、Lで測定される温度と、湯面レベル計9で測定される湯面レベルとを含む測定値を取得し、取得した測定値に基づいて、各観測位置αにおける溶鋼3の温度と鋳造方向(l2軸方向)の流速を導出する。そして、観測データ取得部1413は、可視化対象面Sの範囲内で、図16のステップS303で導出された温度Tの基底ベクトルφT
j(l)と、観測位置αでの温度T(l,t)との内積をとることを、モード番号jのそれぞれについて個別に実行することにより、係数a^T
j(t)を導出する((30)式を参照)。また、観測データ取得部1413は、融液自由表面Cの範囲内で、図16のステップS303で導出された流速vの基底ベクトルφv
j(l)と、観測位置αでの流速v(l,t)との内積をとることを、モード番号jのそれぞれについて個別に実行することにより、係数a^v
j(t)を導出する((31)式を参照)。そして、観測データ取得部1413は、係数a^T
j(t)、a^v
j(t)の全てで構成されるベクトルを、時刻tでの観測データyt
(0)として導出する。
Next, in step S1706, the observation data acquisition unit 1413 acquires measurement values including the temperatures measured by the thermometers F and L and the molten steel level measured by the molten steel level gauge 9, and derives the temperature of the
次に、ステップS1707において、尤度関数導出部1414は、観測データyt
(0)に基づいて、時刻tでの状態ベクトルxtの尤度関数L(xt|yt)を導出する。
次に、ステップS1708において、データ同化部1415は、時刻tでの状態ベクトルxtの第1の確率密度関数p(xt|y0:t-Δt)と、当該時刻tでの状態ベクトルxtの尤度関数L(xt|yt)とを、ベイズ統計のモデリングによるデータ同化を行うフィルタのアルゴリズムに与えることにより、時刻tでの状態ベクトルxtの第2の確率密度関数p(xt|yt0:t)を導出する。
Next, in step S1707, the likelihood
Next, in step S1708, the
次に、ステップS1709において、データ同化部1415は、時刻tでの状態ベクトルxtの第2の確率密度関数p(xt|yt)の最頻値を、時刻tでの状態ベクトルの推定値として導出する。図4のステップS410においては、溶鋼3の流速u、vと温度Tが状態ベクトルの推定値として導出される。これに対し、ステップS1709においては、係数au
1(t)~au
Mu(t)、av
1(t)~av
Mv(t)、aT
1(t)~aT
MT(t)が状態ベクトルの推定値として導出される。
Next, in step S1709, the
次に、ステップS1710において、データ同化部1415は、時刻tが、予め設定された時刻teになったか否かを判定する。この判定の結果、時刻tが、予め設定された時刻teになっていない場合、処理は、ステップS1711に進む。
Next, in step S1710, the
ステップS1711において、第1の係数導出部1411aは、(10)式~(12)式の微分方程式の解(au
1(t+Δt)~au
Mu(t+Δt)、av
1(t+Δt)~av
Mv(t+Δt)、aT
1(t+Δt)~aT
MT(t+Δt))を導出する際に使用するデータである時刻tでの係数au
1(t)~au
Mu(t)、av
1(t)~av
Mv(t)、aT
1(t)~aT
MT(t)と、時刻tから時刻t+Δtまでの間での境界条件のパラメータとの組を、ステップS1708でデータ同化部1415により導出された、時刻tでの状態ベクトルxtの第2の確率密度関数p(xt|yt0:t)に基づいて複数導出する。
In step S1711, the first
そして、処理は、ステップS1703に進み、時刻tが、予め設定された時刻teになるまで、ステップS1703~S1711の処理が繰り返し実行される。尚、ステップS1704においては、ステップS1704の直前のステップS1711で導出された係数と、境界条件のパラメータとの組が用いられる。
以上のようにして、ステップS1710において、時刻tが、予め設定された時刻teになると、処理は、ステップS1712に進む。ステップS1712において、推定物理量導出部1418は、流速u、v、温度Tに対する係数au
1(t)~au
Mu(t)、av
1(t)~av
Mv(t)、aT
1(t)~aT
MT(t)と、図16のステップS303で導出された、流速u、v、温度Tの基底ベクトルφu
1(l)~φu
Mu、φv
1(l)~φv
Mv、φT
1(l)~φT
MTとを、(4)式、(5)式、(6)式にそれぞれ代入することにより、時刻t0~teの各時刻tでの各計算格子点lにおける溶鋼3の流速u(l,t)、v(l,t)、温度T(l,t)を導出する。
Then, the process proceeds to step S1703, and steps S1703 to S1711 are repeatedly executed until time t reaches a preset time t e . In step S1704, the set of the coefficient and the boundary condition parameter derived in step S1711 immediately before step S1704 is used.
In this manner, when the time t reaches the preset time t e in step S1710, the process proceeds to step S1712. In step S1712, the estimated physical quantity derivation unit 1418 derives the flow velocity u(l , t), v( l , t) and temperature T( l , t) of the
尚、図4のステップS410では、時刻tでの各計算格子点lにおける溶鋼3の流速u(l,t)、v(l,t)、温度T(l,t)は、時刻t-Δtでの各計算格子点lにおける溶鋼3の流速u(l,t-Δt)、v(l,t-Δt)および温度T(l,t-Δt)と、時刻t-Δtから時刻tまでの間での境界条件のパラメータとの組のそれぞれについて導出される。これに対し、ステップS1712では、時刻t0~teの各時刻tでの各計算格子点lにおける溶鋼3の流速u(l,t)、v(l,t)、温度T(l,t)が導出される。
4, the flow velocity u(l,t), v(l,t) and temperature T(l,t) of the
次に、ステップS1713において、可視化データ作成部1416は、ステップS1713で導出された時刻t0~teの各時刻tでの各位置における溶鋼3の流速と温度の表示データを作成する。尚、図4のステップS411においては、1つの時刻での各位置における溶鋼3の流速と温度の表示データが作成される。これに対し、ステップS1714では、時刻t0~teの各時刻tでの各位置における溶鋼3の流速と温度の表示データが一括で作成される。
Next, in step S1713, the visualization data creation unit 1416 creates display data of the flow velocity and temperature of the
次に、ステップS1714において、出力部1417は、可視化データ(ステップS1713で導出された各時刻および各位置における溶鋼3の流速と温度のデータと、ステップS1714で導出された表示データ)を出力する。ステップS1714の処理が終了すると、図17のフローチャートによる処理は終了する。
Next, in step S1714, the
(実施例2)
次に、実施例2を説明するが、本発明は、以下の実施例2に限定されるものではない。実施例2では、第2の実施形態と第1の実施形態とを比較する。
Example 2
Next, Example 2 will be described, but the present invention is not limited to the following Example 2. In Example 2, the second embodiment and the first embodiment are compared.
実施例1の発明例1、2では計算格子点lの数をそれぞれ、2600点、3600点とした。これに対し、実施例2では計算格子点lの数を26000点とし、より高精度な数値シミュレーションが実行されるようにした。 In the first and second invention examples of the first embodiment, the number of computational grid points l was set to 2600 and 3600, respectively. In contrast, in the second embodiment, the number of computational grid points l was set to 26000, allowing for a more accurate numerical simulation to be performed.
また、実機のシステムでは、全ての観測位置での測定データ(可視化対象面の各位置での溶鋼3の温度T、可視化対象面の溶鋼3の湯面に対応する各位置での流速v)を取得し、データ同化の計算機システムである情報処理装置100、1300に転送するのに0.5秒程度かかる。このことから、1回のデータ同化のステップを終えた後、次の測定データが取得されるまでの時間を0.5秒とした。第1の実施形態の図4のフローチャートでは、この時間は、ステップS409の処理が終了してから次のステップS407の処理が終了するまでの時間である。第2の実施形態の図17のフローチャートでは、この時間は、ステップS1709の処理が終了してから次のステップS1706において測定値(温度計F、Lで測定される温度と、湯面レベル計9で測定される湯面レベル)が取得されるまでの時間である。また、実施例2では、流速u、v、温度Tの基底ベクトルφu
j(l)、φv
j(l)、φT
j(l)の数Mu、Mv、MTを、それぞれ、8個、8個、10個とした。
In addition, in the actual system, it takes about 0.5 seconds to acquire measurement data (temperature T of the
この他の実施例2の計算条件は、実施例1の発明例1の計算条件と同じである。実施例1で使用した測定データと同じ測定データを使用して、第1の実施形態で説明した手法および第2の実施形態で説明した手法のそれぞれにおいて可視化データを作成した。以下では、第2の実施形態で説明した手法を発明例3と称し、第1の実施形態で説明した手法を参考例1と称する。 The other calculation conditions of Example 2 are the same as those of Example 1. The same measurement data as that used in Example 1 was used to create visualization data in the method described in the first embodiment and the method described in the second embodiment. Hereinafter, the method described in the second embodiment will be referred to as Example 3, and the method described in the first embodiment will be referred to as Reference Example 1.
その結果、発明例3では、各位置における溶鋼3の流速と温度の計算精度は、実施例1の発明例1および発明例2と同等になった。
次に、発明例3と参考例1の計算時間を比較する。表3は、発明例3の計算時間を示す。表4は、参考例1の計算時間を示す。発明例3の計算時間は、図17のフローチャートの処理を実行するのに要する時間であり、参考例1の計算時間は、図4のフローチャートの処理を実行するのに要する時間である。尚、表3および表4に示す数値の単位は秒である。
As a result, in Example 3, the calculation accuracy of the flow velocity and temperature of the
Next, the calculation times of Example 3 and Reference Example 1 are compared. Table 3 shows the calculation times of Example 3. Table 4 shows the calculation times of Reference Example 1. The calculation times of Example 3 are the time required to execute the process of the flowchart in Fig. 17, and the calculation times of Reference Example 1 are the time required to execute the process of the flowchart in Fig. 4. The numerical values shown in Tables 3 and 4 are in seconds.
表3および表4において、想定時間は、実施例1の表1に示した想定時間と同じ意味であり、シミュレーションで想定している時間である。実操業において想定時間の欄に示される時間分のシミュレーションを実行するために必要な時間が、表3の発明例3の計算時間の合計の欄および表4の参考例1の計算時間の合計の欄に記載された時間である。 In Tables 3 and 4, the estimated time has the same meaning as the estimated time shown in Table 1 of Example 1, and is the time assumed in the simulation. The time required to execute the simulation for the time shown in the estimated time column in actual operation is the time shown in the column for total calculation time for invention example 3 in Table 3 and the column for total calculation time for reference example 1 in Table 4.
表3および表4において、測定値収集は、図17、図4のフローチャートの処理を開始してから終了するまでの期間において測定値(可視化対象面の各位置での溶鋼3の温度T、可視化対象面の溶鋼3の湯面に対応する各位置での流速v)を取得するのに要する時間である。前述したように、本実施例では、1回のデータ同化のステップを終えた後、次の測定データが取得されるまでの時間を0.5秒としている。従って、この時間(=0.5秒)に想定時間を乗算した値が測定値収集の欄の値になる。
In Tables 3 and 4, measurement value collection is the time required to obtain measurement values (temperature T of
PODSimは、オンラインにおいて、係数au
1(t)~au
Mu(t)、av
1(t)~av
Mv(t)、aT
1(t)~aT
MT(t)を導出するのに要する時間(第1の係数導出部111a、1411aにおける処理時間)であり、実施例1の表1に示した発明例1および発明例2の欄に示す時間に対応する。
その他は、オンラインにおけるその他の処理時間である。
PODSim is the time required to derive the coefficients a u 1 (t) to a u Mu (t), a v 1 (t) to a v Mv (t), and a T 1 (t) to a T MT (t) online (the processing time in the first
The other is other online processing time.
表3および表4のPODSimの欄に示す時間と表1の発明例1および発明例2の欄に示す時間とを比較すると、計算格子点lの数に関わらず(発明例1では2600点、発明例2では3600点、発明例3および参考例1では26000点)、縮約シミュレータの稼働に要する時間(係数導出部111a、1411aにおける処理時間)は略同じである。
When comparing the time shown in the PODSim column of Tables 3 and 4 with the time shown in the columns of Example 1 and Example 2 of Table 1, the time required to operate the contraction simulator (processing time in
表3の発明例3の計算時間の合計の欄に示す時間は、何れも想定時間を下回っており、発明例3(第2の実施形態の手法)では、想定時間内に可視化データを作成して出力することができることが分かる。一方、表4の参考例1の計算時間の合計の欄に示す時間は、何れも想定時間を上回っており、参考例1(第1の実施形態の手法)では、計算格子点lの数が多くなると、想定時間内に可視化データを作成して出力することができなくなる虞があることが分かる。 The times shown in the column for total calculation time for invention example 3 in Table 3 are all shorter than the estimated time, and it is clear that in invention example 3 (method of the second embodiment), visualization data can be created and output within the estimated time. On the other hand, the times shown in the column for total calculation time for reference example 1 in Table 4 are all longer than the estimated time, and it is clear that in reference example 1 (method of the first embodiment), if the number of computational grid points l becomes large, there is a risk that visualization data cannot be created and output within the estimated time.
表4の結果は、参考例1では、ステップS404の処理の一部(時刻tでの各計算格子点lにおける溶鋼3の流速u(l,t-Δt)、v(l,t-Δt)、温度T(l,t-Δt)に基づいて、流速u、v、温度Tに対する係数au
1(t-Δt)~au
Mu(t-Δt)、av
1(t-Δt)~av
Mv(t-Δt)、aT
1(t-Δt)~aT
MT(t-Δt)を導出すること)とステップS405の処理(流速u、v、温度Tに対する係数au
1(t)~au
Mu(t)、av
1(t)~av
Mv(t)、aT
1(t)~aT
MT(t)に基づいて、時刻tでの各計算格子点lにおける溶鋼3の流速u(l,t)、v(l,t)、温度T(l,t)を導出すること)を、全てのアンサンブルメンバーに対してステップS409においてデータ同化を行う度(時刻tを更新する度)に繰り返し実行するために、計算格子点lの数が多くなると、計算時間が長くなることに対応する。
The results in Table 4 are obtained in Reference Example 1 by deriving coefficients a u 1 (t-Δt) to a u Mu (t-Δt), a v 1 (t-Δt) to a v Mv (t-Δt), and a T 1 (t-Δt) to a T MT (t-Δt) for the flow velocities u, v, and temperature T based on the flow velocities u ( l ,t-Δt), v ( l ,t-Δt), and temperature T ( l ,t-Δt) of the molten steel 3 at each computational grid point l at time t) in part of the processing in step S404 and processing in step S405 (deriving coefficients a u 1 (t) to a u Mu (t), a v 1 (t) to a v Mv (t), and a T 1 (t) to a T MT Since the process of deriving the flow velocity u(l, t), v(l, t), and temperature T(l, t) of the
表3の結果は、発明例3では、ステップS1706の処理の一部(時刻tでの各観測位置αにおける溶鋼3の温度Tと流速vとに基づいて、温度T、流速vに対する係数a^v 1(t)~a^v Mv(t)、a^T 1(t)~a^T MT(t)を導出すること)を流速vおよび温度Tの測定値のみに対して実行し、ステップS1713の処理(流速u、v、温度Tに対する係数au 1(t)~au Mu(t)、av 1(t)~av Mv(t)、aT 1(t)~aT MT(t)に基づいて、時刻tでの各計算格子点lにおける溶鋼3の流速u(l,t)、v(l,t)、温度T(l,t)を導出すること)を、ステップS1709のデータ同化が全て終了した後(時刻teの後)に状態ベクトルの推定値のみに対して実行するために、計算格子点lの数が多くなっても、参考例1に比べ、計算時間が短くなる。 The results in Table 3 show that, in invention example 3, a part of the process of step S1706 (deriving coefficients a^ v1 (t) to a^ vMv (t) and a^ T1 (t) to a^ TMT (t) for temperature T and flow velocity v based on the temperature T and flow velocity v of the molten steel 3 at each observation position α at time t) is executed only for the measured values of the flow velocity v and temperature T, and the process of step S1713 (deriving the flow velocity u ( l ,t), v ( l ,t) and temperature T( l,t) of the molten steel 3 at each computational grid point l at time t based on the coefficients au1(t) to auMu (t), av1(t) to avMv ( t ) , and aT1 (t) to aTMT(t) for the flow velocities u, v, and temperature T) is executed after all data assimilation in step S1709 is completed (at time t Since the calculation is performed only on the estimated value of the state vector (after e ), the calculation time is shorter than that of the first embodiment even if the number of calculation grid points l is increased.
(まとめ)
以上のように本実施形態では、情報処理装置1400は、係数au
1(t)~au
Mu(t)、av
1(t)~av
Mv(t)、aT
1(t)~aT
MT(t)の全てで構成される状態ベクトルを、時刻tでの数値解析データ{xt|t-Δt
(k)}kとしてオンラインで導出する。また、情報処理装置1400は、係数a^T
j(t)、a^v
j(t)の全てで構成される観測ベクトルを、時刻tでの観測データyt
(0)としてオンラインで導出する。情報処理装置1400は、これら数値解析データ{xt|t-Δt
(k)}kおよび観測データyt
(0)を用いてデータ同化を実行する。そして、情報処理装置1400は、データ同化により導出した状態ベクトルの推定値を構成する係数au
1(t)~au
Mu(t)、av
1(t)~av
Mv(t)aT
1(t)~aT
MT(t)を用いて、基底ベクトルφu
1(l)~φu
Mu(l)、φv
1(l)~φv
Mv(l)、φT
1(l)~φT
MT(l)を線形結合することにより、各時刻tでの各計算格子点lにおける溶鋼3の流速u、vおよび温度Tを導出する。従って、全てのアンサンブルメンバーに対してデータ同化を実行する度に、係数au
1(t)~au
Mu(t)、av
1(t)~av
Mv(t)aT
1(t)~aT
MT(t)を、各計算格子点lにおける溶鋼3の流速u、vおよび温度Tに変換する必要がなくなる。従って、第1の実施形態で説明した効果に加えて、オンラインにおける計算時間をより短くすることができるという効果が得られる。
(summary)
As described above, in this embodiment, the
また、本実施形態では、情報処理装置1400は、スナップショット行列MU、MVに対して特異値分解を実行することにより、基底ベクトルφu
1(l)~φu
Mu(l)、φv
1(l)~φv
Mv(l)、φT
1(l)~φT
MT(l)を線形結合する際に当該基底ベクトルφu
1(l)~φu
Mu(l)、φv
1(l)~φv
Mv(l)、φT
1(l)~φT
MT(l)に乗算される係数au
1(t)~au
Mu(t)、av
1(t)~av
Mv(t)aT
1(t)~aT
MT(t)を導出し、これらの相関に基づく成分を有する観測行列Htを導出する。従って、係数a^v
j(t)、a^T
j(t)のそれぞれ異なる物理量に対する係数au
1(t)~au
Mu(t)、av
1(t)~av
Mv(t)aT
1(t)~aT
MT(t)への依存度を考慮して、観測行列Htの成分を導出することができる。よって、データ同化の計算精度をより向上させることができる。
尚、本実施形態においても第1の実施形態で説明した種々の変形例を採用してもよい。
Furthermore, in this embodiment, the
Incidentally, in this embodiment as well, the various modified examples described in the first embodiment may be adopted.
尚、以上説明した本発明の実施形態は、コンピュータがプログラムを実行することによって実現することができる。また、前記プログラムを記録したコンピュータ読み取り可能な記録媒体及び前記プログラム等のコンピュータプログラムプロダクトも本発明の実施形態として適用することができる。記録媒体としては、例えば、フレキシブルディスク、ハードディスク、光ディスク、光磁気ディスク、CD-ROM、磁気テープ、不揮発性のメモリカード、ROM等を用いることができる。
また、以上説明した本発明の実施形態は、何れも本発明を実施するにあたっての具体化の例を示したものに過ぎず、これらによって本発明の技術的範囲が限定的に解釈されてはならないものである。すなわち、本発明はその技術思想、またはその主要な特徴から逸脱することなく、様々な形で実施することができる。
The above-described embodiment of the present invention can be realized by a computer executing a program. A computer-readable recording medium on which the program is recorded and a computer program product such as the program can also be applied as an embodiment of the present invention. Examples of the recording medium that can be used include a flexible disk, a hard disk, an optical disk, a magneto-optical disk, a CD-ROM, a magnetic tape, a non-volatile memory card, and a ROM.
Furthermore, the above-described embodiments of the present invention are merely examples of the implementation of the present invention, and the technical scope of the present invention should not be interpreted as being limited by these. In other words, the present invention can be implemented in various forms without departing from its technical concept or main features.
1:取鍋、2:タンディッシュ、3:溶鋼、4:鋳型、5:スライディングノズル、6:浸漬ノズル、7:吐出口、8:凝固シェル、9:湯面レベル計、100、1400:情報処理装置、101:数値シミュレーション部、102:行列導出部、103:基底ベクトル導出部、104:縮約シミュレータ導出部、111:状態導出部、111a:第1の係数導出部、111b:物理量導出部、112:第1の確率密度関数導出部、113:観測データ取得部、114:尤度関数導出部、115:データ同化部、116:可視化データ作成部、117:出力部、1405:第2の係数導出部、1406:観測行列導出部、1411:状態導出部、1411a:第1の係数導出部、1412:第1の確率密度関数導出部、1413:観測データ取得部、1414:尤度関数導出部、1415:データ同化部、1416:可視化データ作成部、1417:出力部、1418:推定物理量導出部1418、F1~F12、L1~L12:温度計 1: Ladle, 2: Tundish, 3: Molten steel, 4: Mold, 5: Sliding nozzle, 6: Submerged nozzle, 7: Discharge outlet, 8: Solidified shell, 9: Molten metal level gauge, 100, 1400: Information processing device, 101: Numerical simulation unit, 102: Matrix derivation unit, 103: Basis vector derivation unit, 104: Reduced simulator derivation unit, 111: State derivation unit, 111a: First coefficient derivation unit, 111b: Physical quantity derivation unit, 112: First probability density function derivation unit, 113: Observation data acquisition unit, 114: Likelihood function derivation unit, 115: data assimilation unit, 116: visualization data creation unit, 117: output unit, 1405: second coefficient derivation unit, 1406: observation matrix derivation unit, 1411: state derivation unit, 1411a: first coefficient derivation unit, 1412: first probability density function derivation unit, 1413: observation data acquisition unit, 1414: likelihood function derivation unit, 1415: data assimilation unit, 1416: visualization data creation unit, 1417: output unit, 1418: estimated physical quantity derivation unit 1418, F1 to F12, L1 to L12: thermometer
Claims (12)
各計算位置および各時刻における前記溶融金属の状態を示す物理量である第1の物理量の値を格納する行列を導出する行列導出手段と、
前記行列導出手段により導出された、前記第1の物理量の値を格納する行列に対して特異値分解を実行することにより、1次~M次のモードの基底ベクトルを導出することで、前記第1の物理量の1次~M次のモードの基底ベクトルを導出する基底ベクトル導出手段と、
前記溶融金属の状態を示す量である状態量を導出する状態導出手段と、
を有し、
前記第1の物理量は、前記溶融金属の流速と、前記溶融金属の温度と、を含み、
前記行列導出手段は、前記溶融金属の流速の値を格納する行列を、前記流速の軸成分ごとに導出することと、前記溶融金属の温度の値を格納する行列を導出することと、を行い、
前記基底ベクトル導出手段は、前記行列導出手段により導出された、前記溶融金属の流速の値を格納する行列に対して特異値分解を実行することにより、1次~M次のモードの基底ベクトルを導出することを、前記行列導出手段により導出された、前記溶融金属の流速の値を格納する行列のそれぞれについて実行することで、前記溶融金属の流速の1次~M次のモードの基底ベクトルを、前記流速の軸成分ごとに導出することと、前記行列導出手段により導出された、前記溶融金属の温度の値を格納する行列に対して特異値分解を実行することにより、1次~M次のモードの基底ベクトルを導出することで、前記溶融金属の温度の1次~M次のモードの基底ベクトルを導出することと、を行い、
前記状態導出手段は、前記基底ベクトル導出手段により前記第1の物理量の1次~M次のモードの基底ベクトルが導出された後に、前記第1の物理量の1次~M次のモードの基底ベクトルを線形結合する際に用いる係数を導出する第1の係数導出手段を有し、
前記状態量は、前記第1の係数導出手段により導出された前記係数に基づいて定められ、
前記第1の係数導出手段は、前記第1の物理量の1次~M次のモードの基底ベクトルを線形結合する際に用いる係数として、少なくとも前記溶融金属の流速の1次~M次のモードの基底ベクトルを線形結合する際に用いる係数を、軸成分ごとに導出することと、前記基底ベクトル導出手段により前記溶融金属の温度の1次~M次のモードの基底ベクトルが導出された後に、前記溶融金属の温度の1次~M次のモードの基底ベクトルを線形結合する際に用いる係数を導出することと、を行い、
前記Mの値は、2以上の整数であり、且つ、前記流速の軸成分ごと、前記温度ごとに定まり、
前記溶融金属の流速は、前記基底ベクトル導出手段により導出された、前記溶融金属の流速の1次~M次のモードの基底ベクトルを、前記第1の係数導出手段により導出された係数を用いて線形結合することにより表され、
前記溶融金属の温度は、前記基底ベクトル導出手段により導出された、前記溶融金属の温度の1次~M次のモードの基底ベクトルを、前記第1の係数導出手段により導出された係数を用いて線形結合することにより表されることを特徴とする情報処理装置。 An information processing device that executes a process including deriving a state of molten metal poured into a mold of a continuous casting facility,
a matrix derivation means for deriving a matrix storing values of a first physical quantity, which is a physical quantity indicating a state of the molten metal at each calculation position and each time;
a basis vector derivation means for deriving basis vectors of first to M-th modes of the first physical quantity by performing singular value decomposition on a matrix storing the values of the first physical quantity derived by the matrix derivation means, thereby deriving basis vectors of first to M-th modes of the first physical quantity;
A state derivation means for deriving a state quantity which is a quantity indicating a state of the molten metal;
having
the first physical quantity includes a flow velocity of the molten metal and a temperature of the molten metal;
the matrix deriving means derives a matrix storing values of a flow velocity of the molten metal for each axial component of the flow velocity, and derives a matrix storing values of a temperature of the molten metal;
the basis vector derivation means derives basis vectors of 1st to Mth modes by performing singular value decomposition on a matrix storing values of the molten metal flow velocity derived by the matrix derivation means, for each of the matrices storing values of the molten metal flow velocity derived by the matrix derivation means, thereby deriving basis vectors of 1st to Mth modes of the molten metal flow velocity for each axial component of the flow velocity; and derives basis vectors of 1st to Mth modes of the molten metal temperature by performing singular value decomposition on a matrix storing values of the molten metal temperature derived by the matrix derivation means, thereby deriving basis vectors of 1st to Mth modes of the molten metal temperature;
the state derivation means includes a first coefficient derivation means for deriving a coefficient used in linearly combining the basis vectors of the first to M-th modes of the first physical quantity after the basis vectors of the first to M-th modes of the first physical quantity are derived by the basis vector derivation means;
the state quantity is determined based on the coefficient derived by the first coefficient derivation means,
the first coefficient derivation means derives, for each axial component, a coefficient used when linearly combining basis vectors of a first to Mth mode of at least the molten metal flow velocity as a coefficient used when linearly combining basis vectors of a first to Mth mode of the first physical quantity, and derives a coefficient used when linearly combining the basis vectors of a first to Mth mode of the molten metal temperature after the basis vectors of a first to Mth mode of the molten metal temperature have been derived by the basis vector derivation means;
The value of M is an integer of 2 or more, and is determined for each axial component of the flow velocity and for each temperature ,
the flow velocity of the molten metal is expressed by linearly combining basis vectors of first to Mth modes of the flow velocity of the molten metal derived by the basis vector derivation means, using coefficients derived by the first coefficient derivation means ;
The information processing device is characterized in that the temperature of the molten metal is represented by linearly combining basis vectors of 1st to Mth modes of the temperature of the molten metal derived by the basis vector derivation means, using coefficients derived by the first coefficient derivation means .
前記第1の微分方程式のそれぞれは、前記溶融金属の流動を記述する運動方程式であって、前記溶融金属の流速を、1次~M次のモードの基底ベクトルの線形結合で表した運動方程式の両辺に対し、1次~M次のモードの基底ベクトルとの内積をとることにより構成される微分方程式であり、
前記第2の微分方程式のそれぞれは、前記溶融金属の熱伝導を記述する熱伝導方程式であって、前記溶融金属の温度を、1次~M次のモードの基底ベクトルの線形結合で表した熱伝導方程式の両辺に対し、1次~M次のモードの基底ベクトルとの内積をとることにより構成される微分方程式であることを特徴とする請求項1に記載の情報処理装置。 the first coefficient derivation means derives coefficients used when linearly combining the basis vectors of the 1st to Mth modes of the molten metal flow velocity for each axial component of the flow velocity by simultaneously solving first differential equations, each of which is a differential equation expressing a change over time of the coefficients to be multiplied to the basis vectors of the 1st to Mth modes of the molten metal flow velocity, and configured for each axial component of the flow velocity ; and derives coefficients used when linearly combining the basis vectors of the 1st to Mth modes of the molten metal temperature by simultaneously solving second differential equations, each of which is a differential equation expressing a change over time of the coefficients to be multiplied to the basis vectors of the 1st to Mth modes of the molten metal temperature;
each of the first differential equations is a motion equation describing a flow of the molten metal, and is a differential equation obtained by taking an inner product of both sides of an equation of motion expressing a flow velocity of the molten metal as a linear combination of basis vectors of first to Mth modes and the basis vectors of first to Mth modes ;
The information processing device according to claim 1, characterized in that each of the second differential equations is a heat conduction equation describing the heat conduction of the molten metal, and is a differential equation constructed by taking the inner product of both sides of a heat conduction equation in which the temperature of the molten metal is expressed as a linear combination of basis vectors of first to Mth modes and the basis vectors of first to Mth modes.
前記行列導出手段は、前記行列として、前記数値シミュレーション手段により導出された前記第1の物理量であって、前記時間隔Δt1よりの長い時間隔Δt2で定まる各時刻での前記第1の物理量の値を格納するスナップショット行列を導出することを特徴とする請求項1または2に記載の情報処理装置。 a numerical simulation means for deriving the first physical quantity at each time interval Δt1 by a numerical simulation;
3. The information processing device according to claim 1, wherein the matrix derivation means derives, as the matrix, a snapshot matrix storing values of the first physical quantity derived by the numerical simulation means at each time point determined by a time interval Δt2 longer than the time interval Δt1 .
前記状態導出手段は、前記基底ベクトル導出手段により導出された、前記溶融金属の流速の1次~M次のモードの基底ベクトルを、前記第1の係数導出手段により導出された係数を用いて線形結合することにより、前記第1の物理量として前記溶融金属の流速を軸成分ごとに導出することと、前記基底ベクトル導出手段により導出された、前記溶融金属の温度の1次~M次のモードの基底ベクトルを、前記第1の係数導出手段により導出された係数を用いて線形結合することにより、前記第1の物理量として前記溶融金属の温度を更に導出することと、を行う物理量導出手段を更に有することを特徴とする請求項1~4の何れか1項に記載の情報処理装置。 the state quantity includes the first physical quantity,
The information processing device according to any one of claims 1 to 4, characterized in that the state derivation means further includes physical quantity derivation means for linearly combining basis vectors of 1st to Mth modes of the flow velocity of the molten metal derived by the basis vector derivation means, using coefficients derived by the first coefficient derivation means, to derive the flow velocity of the molten metal for each axial component as the first physical quantity, and for further deriving the temperature of the molten metal as the first physical quantity by linearly combining basis vectors of 1st to Mth modes of the temperature of the molten metal derived by the basis vector derivation means, using coefficients derived by the first coefficient derivation means.
前記情報処理装置は、前記状態導出手段により導出された前記複数のケースの時刻tでの前記数値解析データから当該時刻tでの前記状態ベクトルの第1の確率密度関数を導出する第1の確率密度関数導出手段と、
センサより測定される前記連続鋳造設備の各観測位置における物理量である第2の物理量に基づく量である観測量からなるベクトルである観測ベクトルについて、時刻tでの前記観測ベクトルのデータである観測データを前記時間隔Δtごとに取得する観測データ取得手段と、
前記観測データ取得手段により取得された時刻tでの観測データから当該時刻tでの前記状態ベクトルの尤度関数を導出する尤度関数導出手段と、
前記第1の確率密度関数導出手段により導出された時刻tでの前記状態ベクトルの前記第1の確率密度関数と、前記尤度関数導出手段により導出された当該時刻tでの前記状態ベクトルの前記尤度関数とに基づいて、ベイズ統計のモデリングによるデータ同化を行うフィルタにより、当該時刻tでの前記状態ベクトルの第2の確率密度関数を導出し、当該状態ベクトルの前記第2の確率密度関数に基づいて、当該時刻tでの前記状態ベクトルの推定値を導出するデータ同化手段と、
を更に有し、
前記第1の係数導出手段は、前記状態ベクトルの前記第2の確率密度関数に基づいて、前記複数のケースのそれぞれについて、時刻tでの前記係数をオンラインで導出することを特徴とする請求項1~5の何れか1項に記載の情報処理装置。 the state derivation means derives, for a state vector, which is a vector consisting of the state quantities at each calculation position in a visualization target area of the molten metal poured into a mold of the continuous casting equipment, numerical analysis data, which is data of the state vector at time t, for each time interval Δt, for each of a plurality of cases in which at least one of at least one of the first physical quantities at time t-Δt, which is the time interval Δt before time t, and at least one of parameters of a boundary condition between time t-Δt and time t, is different, online;
The information processing device includes a first probability density function derivation means for deriving a first probability density function of the state vector at time t from the numerical analysis data at time t of the plurality of cases derived by the state derivation means;
an observation data acquiring means for acquiring observation data, which is data of an observation vector at time t, for each time interval Δt, the observation vector being a vector of observation quantities that are quantities based on a second physical quantity that is a physical quantity measured by a sensor at each observation position of the continuous casting equipment;
a likelihood function derivation means for deriving a likelihood function of the state vector at time t from the observation data at time t acquired by the observation data acquisition means;
data assimilation means for deriving a second probability density function of the state vector at the time t by a filter that performs data assimilation through modeling of Bayesian statistics based on the first probability density function of the state vector at the time t derived by the first probability density function derivation means and the likelihood function of the state vector at the time t derived by the likelihood function derivation means, and for deriving an estimate of the state vector at the time t based on the second probability density function of the state vector;
Further comprising:
The information processing device according to any one of claims 1 to 5, characterized in that the first coefficient derivation means derives the coefficient at time t online for each of the plurality of cases based on the second probability density function of the state vector.
前記観測量は、前記第2の物理量であり、
前記状態導出手段は、前記基底ベクトル導出手段により導出された、前記第1の物理量の1次~M次のモードの基底ベクトルを、前記第1の係数導出手段により導出された係数を用いて線形結合することにより、前記第1の物理量を導出する物理量導出手段を更に有し、
前記物理量導出手段は、前記複数のケースのそれぞれにおける時刻tでの前記係数を用いて、当該時刻tでの前記数値解析データを前記複数のケースのそれぞれについてオンラインで導出することを特徴とする請求項6に記載の情報処理装置。 the state quantity includes the first physical quantity,
the observable is the second physical quantity,
the state derivation means further includes physical quantity derivation means for linearly combining basis vectors of 1st to Mth modes of the first physical quantity derived by the basis vector derivation means using a coefficient derived by the first coefficient derivation means, to derive the first physical quantity;
7. The information processing apparatus according to claim 6, wherein the physical quantity deriving means derives the numerical analysis data at time t online for each of the plurality of cases by using the coefficient at time t in each of the plurality of cases.
前記状態量は、前記第1の物理量を1次~M次のモードの基底ベクトルの線形結合で表現する際に用いる係数を含み、
前記観測量は、前記第2の物理量を1次~M次のモードの基底ベクトルの線形結合で表現する際に用いる係数であり、
前記基底ベクトル導出手段により導出された、前記第1の物理量の1次~M次のモードの基底ベクトルを、前記データ同化手段により導出された前記状態ベクトルの推定値に含まれる前記係数を用いて線形結合することにより、前記第1の物理量を導出する推定物理量導出手段を更に有し、
前記観測データ取得手段は、センサにより測定された前記各観測位置おける前記第2の物理量の測定値と、前記基底ベクトル導出手段により導出された当該第2の物理量の1次~M次のモードの基底ベクトルとに基づいて前記観測データを導出することを特徴とする請求項6に記載の情報処理装置。 the first physical quantity includes the second physical quantity,
the state quantity includes a coefficient used when expressing the first physical quantity as a linear combination of basis vectors of first to Mth modes,
the observation quantity is a coefficient used when expressing the second physical quantity as a linear combination of basis vectors of first to Mth modes,
the estimation means for deriving the first physical quantity by linearly combining basis vectors of 1st to Mth modes of the first physical quantity derived by the basis vector derivation means using the coefficient included in the estimated value of the state vector derived by the data assimilation means,
The information processing device according to claim 6, characterized in that the observation data acquisition means derives the observation data based on measurement values of the second physical quantity at each of the observation positions measured by a sensor and basis vectors of 1st to Mth modes of the second physical quantity derived by the basis vector derivation means.
前記状態量である、前記第1の物理量の1次~M次のモードの基底ベクトルを線形結合する際に用いる係数と、前記観測量である、前記第2の物理量の1次~M次のモードの基底ベクトルを線形結合する際に用いる係数との関係を示す観測行列を導出する観測行列導出手段と、を更に有し、
前記観測行列導出手段は、前記第2の係数導出手段により導出された係数に基づいて、前記状態量である、前記第1の物理量の1次~M次のモードの基底ベクトルを線形結合する際に用いる係数と、前記観測量である、前記第2の物理量の1次~M次のモードの基底ベクトルを線形結合する際に用いる係数との相関を示す値を、前記観測行列の成分として導出することを特徴とする請求項8に記載の情報処理装置。 a second coefficient derivation means for deriving coefficients used in linearly combining basis vectors of first to M-th modes of the first physical quantity by performing the singular value decomposition on a matrix storing values of the first physical quantity;
and an observation matrix derivation means for deriving an observation matrix indicating a relationship between a coefficient used in linearly combining basis vectors of a first to an Mth mode of the first physical quantity, which is the state quantity, and a coefficient used in linearly combining basis vectors of a first to an Mth mode of the second physical quantity, which is the observation quantity,
The information processing device according to claim 8, characterized in that the observation matrix derivation means derives, based on the coefficients derived by the second coefficient derivation means, values indicating a correlation between a coefficient used when linearly combining basis vectors of a first to an Mth mode of the first physical quantity, which is the state quantity, and a coefficient used when linearly combining basis vectors of a first to an Mth mode of the second physical quantity, which is the observation quantity, as components of the observation matrix.
各計算位置および各時刻における前記溶融金属の状態を示す物理量である第1の物理量の値を格納する行列を導出する行列導出工程と、
前記行列導出工程により導出された、前記第1の物理量の値を格納する行列に対して特異値分解を実行することにより、1次~M次のモードの基底ベクトルを導出することで、前記第1の物理量の1次~M次のモードの基底ベクトルを導出する基底ベクトル導出工程と、
前記溶融金属の状態を示す量である状態量を導出する状態導出工程と、
を有し、
前記第1の物理量は、前記溶融金属の流速と、前記溶融金属の温度と、を含み、
前記行列導出工程は、前記溶融金属の流速の値を格納する行列を、前記流速の軸成分ごとに導出することと、前記溶融金属の温度の値を格納する行列を導出することと、を行い、
前記基底ベクトル導出工程は、前記行列導出工程により導出された、前記溶融金属の流速の値を格納する行列に対して特異値分解を実行することにより、1次~M次のモードの基底ベクトルを導出することを、前記行列導出工程により導出された、前記溶融金属の流速の値を格納する行列のそれぞれについて実行することで、前記溶融金属の流速の1次~M次のモードの基底ベクトルを、前記流速の軸成分ごとに導出することと、前記行列導出工程により導出された、前記溶融金属の温度の値を格納する行列に対して特異値分解を実行することにより、1次~M次のモードの基底ベクトルを導出することで、前記溶融金属の温度の1次~M次のモードの基底ベクトルを導出することと、を行い、
前記状態導出工程は、前記基底ベクトル導出工程により前記第1の物理量の1次~M次のモードの基底ベクトルが導出された後に、前記第1の物理量の1次~M次のモードの基底ベクトルを線形結合する際に用いる係数を導出する第1の係数導出工程を有し、
前記状態量は、前記第1の係数導出工程により導出された前記係数に基づいて定められ、
前記第1の係数導出工程は、前記第1の物理量の1次~M次のモードの基底ベクトルを線形結合する際に用いる係数として、少なくとも前記溶融金属の流速の1次~M次のモードの基底ベクトルを線形結合する際に用いる係数を、軸成分ごとに導出することと、前記基底ベクトル導出工程により前記溶融金属の温度の1次~M次のモードの基底ベクトルが導出された後に、前記溶融金属の温度の1次~M次のモードの基底ベクトルを線形結合する際に用いる係数を導出することと、を行い、
前記Mの値は、2以上の整数であり、且つ、前記流速の軸成分ごと、前記温度ごとに定まり、
前記溶融金属の流速は、前記基底ベクトル導出工程により導出された、前記溶融金属の流速の1次~M次のモードの基底ベクトルを、前記第1の係数導出工程により導出された係数を用いて線形結合することにより表され、
前記溶融金属の温度は、前記基底ベクトル導出工程により導出された、前記溶融金属の温度の1次~M次のモードの基底ベクトルを、前記第1の係数導出工程により導出された係数を用いて線形結合することにより表されることを特徴とする情報処理方法。 An information processing method for performing a process including deriving a state of molten metal poured into a mold of a continuous casting facility, comprising:
a matrix derivation step of deriving a matrix storing values of a first physical quantity, which is a physical quantity indicating a state of the molten metal at each calculation position and each time;
a basis vector derivation step of deriving basis vectors of first to M-th modes of the first physical quantity by performing singular value decomposition on a matrix storing values of the first physical quantity derived in the matrix derivation step, thereby deriving basis vectors of first to M-th modes of the first physical quantity;
A state derivation step of deriving a state quantity which is a quantity indicating a state of the molten metal;
having
the first physical quantity includes a flow velocity of the molten metal and a temperature of the molten metal;
The matrix deriving step includes deriving a matrix storing values of a flow velocity of the molten metal for each axial component of the flow velocity, and deriving a matrix storing values of a temperature of the molten metal;
the basis vector derivation step performs singular value decomposition on a matrix storing values of the molten metal flow velocity derived in the matrix derivation step to derive basis vectors of 1st to Mth modes, for each of the matrices storing values of the molten metal flow velocity derived in the matrix derivation step, thereby deriving basis vectors of 1st to Mth modes of the molten metal flow velocity for each axial component of the flow velocity ; and performs singular value decomposition on a matrix storing values of the molten metal temperature derived in the matrix derivation step to derive basis vectors of 1st to Mth modes of the molten metal temperature, thereby deriving basis vectors of 1st to Mth modes of the molten metal temperature .
the state derivation step includes a first coefficient derivation step of deriving coefficients used in linearly combining basis vectors of first to M-th modes of the first physical quantity after the basis vectors of first to M-th modes of the first physical quantity are derived by the basis vector derivation step,
the state quantity is determined based on the coefficient derived in the first coefficient derivation step,
the first coefficient derivation step comprises deriving, for each axial component, a coefficient used in linearly combining basis vectors of a first to an Mth mode of at least the molten metal flow velocity as a coefficient used in linearly combining basis vectors of a first to an Mth mode of the first physical quantity; and deriving a coefficient used in linearly combining the basis vectors of a first to an Mth mode of the molten metal temperature after the basis vectors of a first to an Mth mode of the molten metal temperature have been derived in the basis vector derivation step;
The value of M is an integer of 2 or more, and is determined for each axial component of the flow velocity and for each temperature ,
the flow velocity of the molten metal is expressed by linearly combining basis vectors of first to Mth modes of the flow velocity of the molten metal derived in the basis vector derivation step, using coefficients derived in the first coefficient derivation step ;
The information processing method is characterized in that the temperature of the molten metal is represented by linearly combining basis vectors of 1st to Mth modes of the temperature of the molten metal derived in the basis vector derivation process using coefficients derived in the first coefficient derivation process .
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2021130708A JP7698196B2 (en) | 2021-08-10 | 2021-08-10 | Information processing device, information processing method, and program |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2021130708A JP7698196B2 (en) | 2021-08-10 | 2021-08-10 | Information processing device, information processing method, and program |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2023025449A JP2023025449A (en) | 2023-02-22 |
| JP7698196B2 true JP7698196B2 (en) | 2025-06-25 |
Family
ID=85251417
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2021130708A Active JP7698196B2 (en) | 2021-08-10 | 2021-08-10 | Information processing device, information processing method, and program |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP7698196B2 (en) |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2016016414A (en) | 2014-07-07 | 2016-02-01 | Jfeスチール株式会社 | Method for estimating flow state of molten steel, and flow state estimation device |
| JP2019217510A (en) | 2018-06-15 | 2019-12-26 | 日本製鉄株式会社 | Device for visualizing the inside of continuous casting mold, method, and program |
| JP2021102221A (en) | 2019-12-25 | 2021-07-15 | 日本製鉄株式会社 | Device for visualizing inside of continuous casting mold, method therefor and program therefor |
-
2021
- 2021-08-10 JP JP2021130708A patent/JP7698196B2/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2016016414A (en) | 2014-07-07 | 2016-02-01 | Jfeスチール株式会社 | Method for estimating flow state of molten steel, and flow state estimation device |
| JP2019217510A (en) | 2018-06-15 | 2019-12-26 | 日本製鉄株式会社 | Device for visualizing the inside of continuous casting mold, method, and program |
| JP2021102221A (en) | 2019-12-25 | 2021-07-15 | 日本製鉄株式会社 | Device for visualizing inside of continuous casting mold, method therefor and program therefor |
Also Published As
| Publication number | Publication date |
|---|---|
| JP2023025449A (en) | 2023-02-22 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Aly et al. | Modelling of surface tension force for free surface flows in ISPH method | |
| JP6816794B2 (en) | Flow state estimation method of molten steel, flow state estimation device, online display device of flow state of molten steel and continuous casting method of steel | |
| JP6443165B2 (en) | State estimation method and state estimation device | |
| CN119180176B (en) | Method, device and storage medium for predicting mold life in continuous casting process | |
| CN118536377A (en) | Method and device for predicting flow state of molten steel in continuous casting crystallizer, and storage medium | |
| JP6536384B2 (en) | State estimation method, level control method, program and state estimation device | |
| US10583477B2 (en) | Apparatus, method, and program for detecting molten metal surface level in continuous casting mold | |
| JP2019217510A (en) | Device for visualizing the inside of continuous casting mold, method, and program | |
| CN104769593B (en) | Simulator, emulation mode and program | |
| JP2014146302A (en) | Simulation program, simulation method, and simulation device | |
| JP7698196B2 (en) | Information processing device, information processing method, and program | |
| JP7367733B2 (en) | Method and device for estimating surface defects in steel materials, manufacturing specification determination support method and device for steel materials, and manufacturing method for steel materials | |
| Rozos et al. | Stochastic Analysis and Modeling of Velocity Observations in Turbulent Flows | |
| JP4032848B2 (en) | Molding simulation method, molding simulation apparatus, molding simulation program, and computer-readable recording medium recording the molding simulation program | |
| JP4797157B2 (en) | Method and program for estimating fluid and thermal characteristics of turbulent flow with buoyancy | |
| JP5981901B2 (en) | Airflow estimation method and apparatus | |
| JP6372217B2 (en) | Method and apparatus for estimating state of molten metal level in continuous casting mold | |
| JP7332875B2 (en) | Continuous casting mold visualization device, method, and program | |
| JP2021102221A (en) | Device for visualizing inside of continuous casting mold, method therefor and program therefor | |
| JP5442242B2 (en) | Cast hole analysis method and cast hole analysis program for die casting | |
| JP6464949B2 (en) | State estimation method, hot water level control method, and state estimation device | |
| JP6471632B2 (en) | Mold level estimation method, level level control method and apparatus in mold | |
| JP7335499B2 (en) | Continuous casting mold visualization device, method, and program | |
| Fu et al. | Study on the inverse problem of thin slab continuous casting mold based on fluid-solid coupled heat transfer | |
| JP6097200B2 (en) | Data estimation method and apparatus |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20240423 |
|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20250219 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20250311 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20250410 |
|
| 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: 20250513 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20250526 |
|
| R150 | Certificate of patent or registration of utility model |
Ref document number: 7698196 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |