[論文レビュー] Identification of Convection Heat Transfer Coefficient of Secondary Cooling Zone of CCM based on Least Squares Method and Stochastic Approximation Method

А. А. Ivanova|arXiv (Cornell University)|Mar 24, 2010

Metallurgy and Material Forming被引用数 4

ひとこと要約

本稿では、曲がり傾きのある連続鋳造機（CCM）の二次冷却ゾーンにおける空間的に分布する対流熱伝達係数（CHTC）を同定する二段階手法を提案する。最初の段階では最小二乗法を用いたパラメータ調整を行い、二番目の段階では確率的近似法を用いた運用手順を実施する。本手法は、限られた温度測定データからCHTCを高精度に推定でき、最適なステップサイズ系列を用いることで200回未満の反復で収束を達成し、産業的鋼鉄鋳造プロセスにおける熱モデル精度を顕著に向上させる。

ABSTRACT

The detailed mathematical model of heat and mass transfer of steel ingot of curvilinear continuous casting machine is proposed. The process of heat and mass transfer is described by nonlinear partial differential equations of parabolic type. Position of phase boundary is determined by Stefan conditions. The temperature of cooling water in mould channel is described by a special balance equation. Boundary conditions of secondary cooling zone include radiant and convective components of heat exchange and account for the complex mechanism of heat-conducting due to airmist cooling using compressed air and water. Convective heat-transfer coefficient of secondary cooling zone is unknown and considered as distributed parameter. To solve this problem the algorithm of initial adjustment of parameter and the algorithm of operative adjustment are developed.

研究の動機と目的

• 曲がり傾きのある連続鋳造機（CCM）の二次冷却ゾーンにおける、空間的・時間的に変化する対流熱伝達係数（CHTC）を正確に同定する課題に対処すること。これは、複雑なミスト冷却による影響を受けるためである。
• 限られた温度測定データを用いてCHTCを推定する強固な逆問題手法を開発すること。これにより、経験的または固定パラメータに基づく手法の限界を克服する。
• 二つの異なるアルゴリズムを提案する：複数の測定点における温度データを用いた最小二乗法による初期調整、および単一測定点における温度データを用いた確率的近似法による運用手順。
• 確率的近似法の収束条件を厳密に定義したステップサイズ系列を用いることで、パrameter同定プロセスの収束性と安定性を保証すること。
• 数値シミュレーションを用いて手法を検証し、異なるステップサイズ系列における収束挙動を比較することで、最適なチューニングパラメータを同定すること。

提案手法

• 熱伝達プロセスは、相変化をステファン条件で扱い、非一様な材料特性を考慮した非線形かつ時間依存の放物型偏微分方程式でモデル化される。
• 境界条件には放射と対流の併存熱伝達が含まれ、対流成分は未知の分布関数としてのCHTCパラメータでモデル化される。
• 初期調整では、鋳錠表面の複数の空間的点で測定された温度データを用いて、最小二乗法によりCHTCを推定する。
• 運用手順では、確率的近似アルゴリズムが用いられ、単一測定点の温度データを用いてCHTC推定値を反復的に更新する：α_{j+1} = α_j - k_j(T_j^* - T_j)。
• ステップサイズ系列k_jは収束条件を満たすように定義される：lim k_j = 0、Σk_j = ∞、およびΣ(k_j)^2 < ∞。例としてk_j = a/(b+j)や、残差誤差の符号変化に基づく適応型形式が挙げられる。
• 収束は、αの最後のm個の近似値が互いに小さなε近傍内に収束しているかどうかをチェックすることで判定され、安定化を示す。

実験結果

リサーチクエスチョン

• RQ1 限られた温度測定データしか入手できない状況下で、曲がり傾きのある連続鋳造機（CCM）の二次冷却ゾーンにおける空間的分布する対流熱伝達係数（CHTC）をどのように正確に同定できるか？
• RQ2 CHTC同定において、高速かつ安定な収束を保証する最適なステップサイズ系列は何か？
• RQ3 ステップサイズ系列k_j = a/(b+j)におけるパラメータaとbの選択が、CHTC推定の収束速度および精度にどのように影響するか？
• RQ4 確率的近似法は単一測定点からの温度データからCHTCを信頼性高く同定できるか？また、複数測定点を用いた最小二乗法と比較して、その性能はいかがであるか？
• RQ5 確率的近似アルゴリズムが真のCHTC値に収束する条件は何か？また、初期パラメータ推定値にどの程度感度を示すか？

主な発見

• 確率的近似アルゴリズムは、k_j = a/j（a ∈ [1,2]）のステップサイズ系列を用いることで、200回未満の反復でCHTCの空間的分布を正確に同定し、収束を達成した。収束は安定的かつ高精度であった。
• k_j = a/jの系列において、aが1 ≤ a ≤ 2の範囲にある場合に最適な収束が観察された。a < 1の場合は収束が遅く、a > 2では増幅する振動が生じ、収束が遅くなる傾向を示した。
• k_j = a/(b+j)の系列におけるパラメータbは非負でなければならない。b < 0の場合は誤った調整方向が生じ、bを大きくすると収束速度が低下するため、最適な性能を得るためにはb = 0が推奨される。
• 残差誤差（T_j^* - T_j）の符号変化に基づく適応型ステップサイズ系列は、a ∈ [1,3]の範囲で収束を達成し、固定系列よりもロバストで、より速やかに安定化することが示された。
• 残差の符号変化時にのみn_jを更新する修正適応型系列では、収束は遅くなるが、a ∈ [1.2,1.5]の範囲で20回の反復後には真値の6％以内の精度を達成でき、許容可能な精度を満たしていた。
• 本手法は強力なロバスト性を示しており、広範な初期CHTC推定値の範囲で信頼性高く収束した。これは、産業的逆熱伝達問題への実用的妥当性を裏付けている。