1、放射性废料的处理问题美国原子能委员会以往处理浓缩的放射性废料的方法,一直是把它们装入密封的圆桶里,然后扔到水深为 90 多米的海底。生态学家和科学家们表示担心,怕圆桶下沉到海底时与海底碰撞而发生破裂,从而造成核污染。原子能委员会分辨说这是不可能的。为此工程师们进行了碰撞实验。发现当圆桶下沉速度超过 12.2 m/s 与海底相撞时,圆桶就可能发生碰裂。这样为避免圆桶碰裂,需要计算一下圆桶沉到海底时速度是多少? 这时已知圆桶重量为 239.46 kg,体积为0.2058m3,海水密度为 1035.71kg/m3,如果圆桶速度小于 12.2 m/s 就说明这种方法是安全可靠的,否则就要禁止使用这种方
2、法来处理放射性废料。假设水的阻力与速度大小成正比例,其正比例常数 。现要求建立合理的数学模6.0k型,解决如下实际问题:1. 判断这种处理废料的方法是否合理?2. 一般情况下, 大, 也大; 小, 也小。当 很大时,常用 来代vkvkvkv替 ,那么这时速度与时间关系如何? 并求出当速度不超过 12.2 m/s,圆桶的运k动时间和位移应不超过多少? ( 的值仍设为 0.6)3.1 问题分析与建立模型的建立(1) 首先要找出圆桶的运动规律,由于圆桶在运动过程中受到本身的重力以及水的浮力 F 和水的阻力 f 的作用,所以根据牛顿运动定律得到以下的方程 :(1)G合又因为 以及 ,可得到gVFmdt
3、stvma,2合 dtskvf圆桶的位移和速度分别满足下面的微分方程: )3(22kvgVmdtvdtss(2) 由题设这时圆桶受到的阻力应改为 , 类似上面,可得22dtskvf到这时圆桶的速度应满足如下的微分方程:(4)2kvgVmdtv3.2 模型求解(1) 根据方程(2)和(3),加上初始条件 以及题设的0,0ttt svds初始数据。通过 Matlab 就可以求出圆桶的位移合速度的方程。Matlab 源程序如下:m=239.46; w=0.2058;g=9.8;p=1035.71;k=0.6;dsolve(m*D2s=m*g-p*g*w-k*Ds,t)C3/exp(k*t)/m) -
4、 (C2 - g*m*t + g*p*t*w)/k - (g*m*(m - p*w)/k2 ,0,0,* ;6.71.035;8.9205.;46.239 tsstskwgpmtsDsolveChpgwm ans =(2362544856843450433*t)/5497558138880000 - (5*C2)/3 + C3/exp(10*t)/3991) - 9428916523662210678103/54975581388800000 ,0,* tvtvkwgpmtvDSole 得出位移方程为:(5)tetts 02564.174.29175)( 速度方程为 (6)ttv02564.4
5、)(通过方程(5)及 ,利用下面 Mathematica 程序:ms9013,7502564./175.2 ttExptFindRot 求出时间 ,再把它带入方程(6)求出圆桶的速度为 .s.13 smv/29.显然此圆桶的速度已超过 ,可以得出这种处理废料的方法不合理。s/.因此,美国原子能委员会已经禁止用这种方法来处理放射性废料。(2) 根据方程 (4),加上初始条件 ,利用下面的 Mathematica 程序:0tv,0,2)(* tvtkVgpmtvDSole 求出圆桶的速度: ).()(tFTanhkt若把题设中的 (仍为 0.6 的话)和 m 的值带入上式,可得到圆桶的速度,为)0
6、519436.(73.20)( tTanhtv这时若该速度要小于 ,那么利用 Mathematica 可得圆桶的运动时间就不s/能超过 , 位移不能超过 84.8m。对应的 Mathematica 源程序:s1(求时间)12,.0519426.*73.20 ttTanhFindRot (求位移)3,Itegra从这个模型, 我们也可以得到原来处理核废料的方法是不合理的.3.3 结果分析(1) 由于在实际中 k 与 v 的关系很难确定, 所以上面的模型有它的局限性,而且对不同的介质比如在水中与在空气中 k 与 v 的关系也不同 .如果按照(1)的假设 k 为常数的话,那么水中的这个 k 就比在空气中对应的 k 要大一些.在一般情况下,k 应是 v 的函数,即 k=k(v),至于是什么样的函数,这个问题至今还没有能解决.(2) 这个模型还可以推广到其他方面,比如说一个物体从高空落向地面的道理也是一样的.尽管物体越高,落到地面的速度越大,但决不会无限大.3.4 相关的思考题(1) 若在其他介质如在空气中物体的运行速度与它所受的空气阻力仍成正比, 那么这个正比例系数与在水中的正比例系数是否一样?(2) 把所建立的数学模型应用到求火车的速度问题上, 并说明为何火车的速度不能无限增大?