2011年集训微分方程建模问题.ppt

微分方程建模的若干问题 这样，问题定解的初边值条件为： 根据假设（4），初始时刻的鱼群密度分布函数 u ( x , 0 ) 为已知函数 u 0 ( x ) ： 根据假设（2），休渔边界 x = l 外，鱼群将全部被外 海渔船打尽： 显然，鱼群数向岸上游动的扩散流量 Φ( 0 , t ) 为零 ， 根据假设（6）： 故 0 l x t 这个偏微分方程的初、边值问题也是 适定的 ， 即问题的解是 存在、唯一 的，且 连续依赖 于初边值数据 自由边界问题 自由边界问题是一类较为复杂的偏微分方程问题，这种 类型的问题在各种各样的应用中非常频繁地出现，例如 它可出现在物相变化过程、化学反应过程、生物扩散 过程、土壤冻过程等等的物理、化学现象之中，甚至 还出现在金融衍生物价格计算、抵押贷款评估研究等等 的经济现象之中。 （1）一相 Stefan 问题 考虑一根套在与四周完全绝缘隔热的管子中而正在融 化的细冰棍；其右端为冰，左端为融化而成的水。拟 建立一个融化水区域上任意点处温度随时间演变的模型。 建模假设： （1）假定冰区域温度恒等于零度； （2）假定水区域中热量传导服从 Fourier 定律 ，即 单位时间中高温点到低温点通过单位面积的热 流量大小与与温度关于位置量 x 的下降率成正比 ； 由此可推出以下等式 ： （3）假定水的密度 ρ、比热 c 、热传导系数 k 和 为了融化冰为水的潜热 L 均为常数 。 取细棍的一小段 [ x , x +Δx ] , 设细棍的截面积为 s 0 厘米2 ； 记 q ( x , t ) 为热流密度（卡 / 秒 · 厘米2 , 单位时间内通过 单位面积 的热量）， 则在 Δt 时间内，沿 x 方向流入小段 [ x , x +Δx ] 的总热量数近似为：q ( x , t ) · s 0 ·Δt （卡） , 流出小段 [ x , x +Δx ] 的总热量数近似为： q ( x +Δx , t ) · s 0 ·Δt （卡） , 流入小段与流出小段的热量差使得小段中水的温度升 高，这个热量差可以根据下式计算： （ρ· Δx · s 0 ）·c · [ u ( x , t +Δt ) – u ( x , t )] （卡） ， 这样便可得： 根据 Fourier 定律，有： 在融化而成的水域里 ，水的温度 u ( x , t ) 服从 热传导 方程 ： u t = a2 uxx , x ? ( 0 , s0 ) , t ? ( 0 , +? ) . 为求解这个偏微分方程，还需知道左、右边界值和初值。 在 左边界上 水温为已知函数： u ( 0 , t ) = u 1 ( t ) 0 ; 假定水温的 初值 为已知函数 ： u ( x , 0 ) = u 0 ( x ) ； 由于右边界端处的 热传导，冰在不断融化，故水域的 右边界是一条 移动边界 ，或称为 自由边界 。 这条 自由边界 本身也是需要求解的 未知一元函数！ 0 L 冰 水 x t s0 x = s ( t ) 易知，在移动的右边界 s ( t ) 上水温函数应满足： u ( s ( t ) , t ) = 0 ； 为了决定 自由边界 的位置，还需导出边界上另一个条件 。 t1 t2 t3 t4 每次用药时间间隔尽可能长， 但必须应保证 即 （2）最佳用药时间间隔 T 的计算公式为 进一步可思考的问题：如给药方式为 恒速静脉点滴或口服药片时， 模型将如何变化？又如何求解？（ref. 数学模型（第三版），姜启源 等著 , pp. 153 , 第五章 §5 . 4 ） 野猪的