[发明专利]一种测量结晶器热面热流密度、温度的方法

权利要求书

1.一种测量结晶器热面热流密度、温度的方法，其特征在于包括下述步骤：

步骤一

沿结晶器拉坯方向方向，在结晶器壁内纵剖面内，选取垂直结晶器热面的、高度为H、宽度为d₂的矩形区域ABCD，所述矩形区域的竖直边分别为AB边、CD边，且AB边位于结晶器热面上，CD边位于结晶器壁内；选取矩形区域ABCD后；在CD边上设置一组热电偶，并将该组热电偶计为第一组热电偶，且第一组热电偶位于同一条竖直线上；在在第一组热电偶与其所对应的结晶器热面间设有第二组热电偶；所述H≤结晶器的高度；所述d₂≤结晶器的壁厚；

步骤二

连铸时，以一定的采集频率f测量、存储结晶器壁内，时间[t₁,t₂]热电偶温度；

步骤三

采用二维传热反问题把[t₁,t₂]时间段，结晶器壁测量的温度转换为结晶器热面热流密度、温度与结晶器壁内的温度；其过程如下：

a确定计算域Ω

选取矩形区域ABCD为数学计算域，并记为Ω，同时令B为原点O；

所述Ω有四个边界，其上、下边界AD、BC分别记为和其左、右边界AB、CD分别记为和

b反算法求解传热数学模型

反问题目的是：为Ω的三个边界寻找边界热流密度函数和使得正问题方程封闭、而且使得正问题中计算出热电偶所在位置处(x_m,y_m)的温度值等于热电偶的测量值Y_m；

于是传热反问题简化为目标函数的最小化过程；所述目标函数的表达式为式(1)：

s[Q(∂Ω1,2,3,t)]=Σm=1M∫t1t2{Ym-Tm[t;Q(∂Ω1,2,3,t)]}2dt---(1)]]>

式(1)中：

M为矩形区域ABCD不包括CD边在内，所设置热电偶的数目，

Y_m和分别为热电偶所在位置处(x_m,y_m)测量的温度值和通过正问题计算的温度值；

∂Ω1={(x,y)|0≤x≤d2,y=H},∂Ω2={(x,y)|x=0,0≤y≤H},∂Ω3={(x,y)|0≤x≤d2,y=0},]]>

Q(∂Ω1,2,3)=Q(∂Ω1)∪Q(∂Ω2)∪Q(∂Ω3);]]>

然后采用共轭梯度法求解目标函数的最小值，其过程如下：

第1步

令迭代步数i＝0,在时间段[t₁,t₂]内，假设边界的热流密度函数和为常数函数，其值都为常数，所述常数选自0-2×10⁶中任意一个数值；

第2步

求解计算域Ω内的传热过程，将待求解的问题转为求解传热偏微分的初边界(正)问题；把假设的j＝1 to 3带入下列正问题T(x,y,t)偏微分方程：

cρ∂T∂t=∂∂x(k∂T∂x)+∂∂y(k∂T∂y),(x,y)∈Ω,t∈[t1,t2]---(2a)]]>

k∂T∂y=Q(∂Ω1,t)at∂Ω1---(2b)]]>

-k∂T∂x=Q(∂Ω2,t)at∂Ω2---(2c)]]>

-k∂T∂y=Q(∂Ω3,t)at∂Ω3---(2d)]]>

T(x,y,t)=f(∂Ω4,t)at∂Ω4---(2e)]]>

T(x,y,t)＝T_ini                  for t＝t₁.(2f)

式(2a)-(2f)中：

c为结晶器的热容，其单位为J/kg；

ρ结晶器的密度，其单位为kg/m³；

t为时间，其单位为s；

k为结晶器的导热系数，其单位为J/(m·s·K)；

T_ini为反应计算域Ω内t₁时刻温度分布的函数(函数自变量为x,y)；

为t时刻边界的温度，其值由边界上热电偶测量；或者对相邻的两个热电偶温度进行空间线性插值计算他们之间没有热电偶地方的温度，即：

f(Ω4,t)=Y(ym,t),(2g)Y(ym,t)-Y(ym-1,t)ym-ym-1(y-ym-1),(2h)]]>

式(2g)是指有热电偶时，边界Ω₄的温度，其值为热电偶位置(d₂,y_m)处t时刻温度测量值Y(y_m,t)，且d₂为第一组热电偶到结晶器热面的距离；

式(2h)是指无热电偶时，边界Ω₄的温度，其值为由边界Ω₄上相邻的两个热电偶Y(y_m,t)和Y(y_m-1,t)的温度对空间进行线性插值计算得到；

求解正问题，得到计算域Ω在时间段[t₁,t₂]内的结晶器壁内的温度变化；同时也计算出结晶器热面温度和结晶器内热电偶所在位置处的温度值Tm[t;Qi(∂Ω1,2,3,t)];]]>

第3步.把第2步计算的热电偶所在位置的温度值代入方程(1)求解得目标函数值并判断下面收敛标准是否成立，

s[Q(∂Ω1,2,3,t)]≤ϵ---(3)]]>

式(3)中：

由于热电偶测量时含有误差；所述热电偶本身测量误差的标准差为σ，故可以依据Discrepancy Principle计算收敛容差，得ε＝Mσ²(t₂-t₁)；

如果满足，则认为时间段[t₁,t₂]内结晶器热面热流密度为第2步计算得到的结晶器壁温度变化和结晶器热面温度变化为真实值；否则，进入第4步；

第4步.把第二步计算的热电偶所在位置的温度值代入下列伴随问题偏微分方程，计算得伴随问题λ(x,y,t)：

[∂∂x(k∂λ∂x)+∂∂y(k∂λ∂y)]+Σm=1Mi2{Tm[t;Q(∂Ω1,2,3,t)]-Ym}δ(x-xm)δ(y-ym)=-ρc∂λ∂t(x,y)∈Ω,t∈(t1,t2]---(4a)]]>

k∂λ∂y=Σm=1M12{Tm[t;Q(∂Ω1,2,3,t)]-Ym}δ(x-xm)δ(y-ym)at∂Ω1---(4b)]]>

-k∂λ∂x=0at∂Ω2---(4c)]]>

-k∂λ∂y=Σm=1M32[Tm[t;Q(∂Ω1,2,3,t)]-Ym]δ(x-xm)δ(y-ym)at∂Ω3---(4d)]]>

λ(x,y,t)=0at∂Ω4---(4e)]]>

λ(x,y,t)＝0              for t＝t₂.(4f)

式(4a)-(4d)中：δ(.)为Dirac delta算子；

把伴随问题λ(x,y,t)带入梯度公式(5)中，计算得目标函数的梯度▿s[Q(∂Ωj,t)],]]>

▿s[Q(∂Ωj,t)]=λ(∂Ωj,t).j=1to3---(5)]]>

第5步.把第4步的计算的j＝1 to 3带入共轭系数公式(6)中，计算得共轭系数(j＝1 to 3)：

γji=max(0,γj,DYi,γj,HSi)]]>当i＝1时，γji=0.]]>       (6)

式(6)中：

γj,DYi=∫t1t2∫Ωj{▿s[Qi(∂Ωj,t)]}2dΩjdt∫t1t2∫Ωjdi-1(∂Ωj,t){▿s[Qi(∂Ωj,t)]-▿s[Qi-1(∂Ωj,t)]}dΩjdt]]>

γj,HSi=∫t1t2∫Ωj▿s[Qi(∂Ωj,t)]{▿s[Qi(∂Ωj,t)]-▿s[Qi-1(∂Ωj,t)]}dΩjdt∫t1t2∫Ωjdi-1(∂Ωj,t){▿s[Qi(∂Ωj,t)]-▿s[Qi-1(∂Ωj,t)]}dΩjdt]]>

第6步.把第4步的计算梯度和第5步计算的共轭系数代入搜索方向公式(7)中，计算搜索方向

di(∂Ωj,t)=▿s[Qi(∂Ωj,t)]+γjidi-1(∂Ωj,t)---(7)]]>

第7步.解灵敏度问题偏微分方程

ρc∂ΔT∂t=∂∂x(k∂ΔT∂x)+∂∂y(k∂ΔT∂y),(x,y)∈Ω,t∈[t1,t2]---(8a)]]>

k∂ΔT∂y=ΔQ(∂Ω1,t)at∂Ω1---(8b)]]>

-k∂ΔT∂x=ΔQ(∂Ω2,t)at∂Ω2---(8c)]]>

-k∂ΔT∂y=ΔQ(∂Ω3,t)at∂Ω3---(8d)]]>

ΔT(x,y,t)=0at∂Ω4---(8e)]]>

ΔT(x,y,t)＝0        for t＝t₁.                  (8f)

以第6步计算的搜索方向为已知条件，联立式(8a)-(8f)，并

令ΔQ(∂Ω1,t)=di(∂Ω1,t)]]>和ΔQ(∂Ω2,t)=ΔQ(∂Ω3,t)=0,]]>计算

以第6步计算的搜索方向为已知条件，联立式(8a)-(8f)，并

令ΔQ(∂Ω2,t)=di(∂Ω2,t)]]>和ΔQ(∂Ω1,t)=ΔQ(∂Ω3,t)=0,]]>计算

以第6步计算的搜索方向为已知条件，联立式(8a)-(8f)，并

令ΔQ(∂Ω3,t)=di(∂Ω3,t)]]>和ΔQ(∂Ω1,t)=ΔQ(∂Ω2,t)=0]]>计算

第8步.把第7步的计算的热电偶所在位置处(x_m,y_m)的灵敏度(j＝1 to 3)值代入下面搜索步长方程，计算搜索步长(j＝1 to 3)；

Σm=1M∫t1t2(ΔT1,mi)2ΔT1,miΔT2,miΔT1,miΔT3,miΔT2,miΔT1,mi(ΔT2,mi)2ΔT2,miΔT3,miΔT3,miΔT1,miΔT3,miΔT2,mi(ΔT3,mi)2β1iβ2iβ3idt=Σm=1M∫t1t2(Tmi-Ymi)ΔT1,mi(Tmi-Ymi)ΔT2,mi(Tmi-Ymi)ΔT3,midt---(9)]]>

第9步.把第6步的计算搜索方向和第8步的计算搜索步长代入下面热流密度更新公式，计算新的j＝1 to 3；

qi+1(∂Ωj,t)=Qi(∂Ωj,t)-βjidi(∂Ωj,t)for j=1to3---(10)]]>

第10步.令i＝i+l返回第2步，依次循环；直至满足第3步条件，即认为第9步计算得到的热流密度为真实的结晶器热面热流密度，以及第2步计算得到的结晶器壁温度变化和结晶器热面温度变化为真实值。