2020国赛A题偏微分方程的数值解法


国赛A题根据Fourier热传导定律和牛顿冷却定律建立了热传导模型,模型的微分方程如下:
{ ∂ u ∂ t = k c ρ ∂ 2 u ∂ x 2 ∂ u ∂ x ∣ x = 0 = − h ( T 炉 − u ( 0 , t ) ) k , ∂ u ∂ x ∣ x = l = − h ( T 炉 − u ( l , t ) ) k , u ( x , 0 ) = T 0 . \begin{cases} \frac{\partial u}{\partial t}=\frac{k}{c\rho}\frac{\partial ^2u}{\partial x^2}\\ \frac{\partial u}{\partial x}\vert _{x=0}=-\frac{h(T_炉-u(0,t))}{k},\\ \frac{\partial u}{\partial x}\vert _{x=l}=-\frac{h(T_炉-u(l,t))}{k},\\ u(x,0)=T_0. \end{cases} tu=cρkx22uxux=0=kh(Tu(0,t)),xux=l=kh(Tu(l,t)),u(x,0)=T0.
其中, T 炉 ( x ) = T_炉(x)= T(x)=
{ T 1 − 5 , T 6 − T 1 − 5 5 ( x − 172.5 ) + T 1 − 5 , T 6 , T 7 − T 6 5 ( x − 208 ) + T 6 , T 7 , T 8 − 9 − T 7 5 ( x − 243.5 ) + T 7 , T 8 − 9 , T 10 − 11 − T 8 − 9 5 ( x − 314.5 ) + T 8 − 9 , T 10 − 11 \begin{cases} T_{1-5},\\ \frac{T_6-T_{1-5}}{5}(x-172.5)+T_{1-5},\\ T_6,\\ \frac{T_7-T_6}{5}(x-208)+T_{6},\\ T_7,\\ \frac{T_{8-9}-T_7}{5}(x-243.5)+T_{7},\\ T_{8-9},\\ \frac{T_{10-11}-T_{8-9}}{5}(x-314.5)+T_{8-9},\\ T_{10-11} \end{cases} T15,5T6T15(x172.5)+T15,T6,5T7T6(x208)+T6,T7,5T89T7(x243.5)+T7,T89,5T1011T89(x314.5)+T89,T1011
下面来考虑这个偏微分方程的求解问题,设参数 k c ρ \frac{k}{c\rho} cρk k 1 k_1 k1,参数 h k \frac{h}{k} k
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值